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Abstract 

This paper develops the theory of affine Euler-Poincare and affine Lie-Poisson re- 
ductions and apphes these processes to various examples of complex fluids, including 
Yang- Mills and Hall magnetohydrodynamics for fluids and superfiuids, spin glasses, 
microfluids, and liquid crystals. As a consequence of the Lagrangian approach, the 
variational formulation of the equations is determined. On the Hamiltonian side, the 
associated Poisson brackets are obtained by reduction of a canonical cotangent bundle. 
A Kelvin-Noether circulation theorem is presented and is applied to these examples. 
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1 Introduction 



The equations of motion of an adiabatic compressible fluid are given by 



— + VuU = -gradp, 
at p 



dp 
di 



+ div(pu) = 0, 



as 

'dt 



'1.1] 



div(5u) = 0, 



where p is the mass density, 5* is the entropy density, and p is the pressure. It was shown in 
"iol that this system, as well as its magnetohydrodynamic extension, admit a non-canonical 



Poisson formulation, that is, equation f 1 1.11) can be written as 

/ = {/,M, 

relative to a Hamiltonian function h. It is of great (mathematical and physical) interest 
to obtain these Poisson brackets by a reduction procedure from a canonical Hamiltonian 
formulation on a cotangent bundle. 



In 28 



the non-canonical Poisson bracket associated 
to fll.ip is obtained via Lie-Poisson reduction for a semidirect product group involving the 
diffeomorphisms group of the fluid container and the space of the advected quantities p and 



S. The Lagrangian formulation of these equations is given in [21 . 

In the same spirit, the non-canonical Hamiltonian structure for adiabatic Yang-Mills 



charged fluids discovered in 13| is obtained by reduction from a canonical formulation in 
12l | , by using a Kaluza-Klein point of view involving the automorphism group of the principal 
bundle of the theory. This paper gives also the Euler-Poincare formulation of these equations. 

Non-canonical Hamiltonian structures for a wide class of non-dissipative fluid models 
were derived in 18|, IJ], 19|], 20|, and 15|]. These examples include Yang-Mills magne- 



tohydrodynamics, spin glasses, and various models of superfluids, and involve Lie-Poisson 
brackets with cocycles. Remarkably, from a mathematical point of view, the Hamiltonian 
structures of many of these models are identical. This Hamiltonian structure together with 
the corresponding variational principles are studied in more detail, with an application to 



liquid crystals, in [16| . We will refer to all these models as complex fluids. 
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In this paper we show the remarkable property that these Lie-Poisson brackets with 
cocycles can also be obtained by Poisson reduction from a canonical Hamiltonian structure. 
The cocycle in the Hamiltonian structure appears only after reduction and it is due to the 
presence of an affine term added to the cotangent lifted action. The associated reduction 
process is naturally called affine Lie-Poisson reduction. 

An important example of such an affine action is given by the usual action of the au- 
tomorphism group of a principal bundle on the connection forms. As a result, we obtain, 
in a natural way, covariant differentials and covariant divergences in the expression of the 
Poisson brackets and of the reduced equations. These gauge theory aspects in the case of 
complex fluids are mathematically and physically interesting since they represent a bridge 
to other possible gauge theories in physics. 

The new affine Lie-Poisson reduction principle presented in this paper unifies all previous 
approaches in the problem of complex fluid dynamics. Eringen and Holm have extensively 
studied the roles of auxiliary variables that were Lie algebra valued 1-forms frozen into the 
fluid flow. In particular. Holm observed that these variables contribute generalized 2-cocycles 
in the Lie-Poisson structures of ideal complex fluids. We implement Holm's observation as 
a fundamental principle which, for physical interactions in complex fluids, is similar to the 
principle of spontaneous symmetry breaking in the coupling of charges and flelds in electro- 
magnetism. This new geometrically formulated principle applies also to Lie algebra valued 
charges (e.g., spins, Yang-Mills charges, etc). The present paper reformulates the theory of 
complex fluids in the framework of the fundamental principles of condensed-matter physics. 
It is clear that other geometric applications in condensed-matter physics are forthcoming, 
such as micropolar elasticity, for example. 

The paper also formulates a parallel theory on the Lagrangian side by extending the 



Euler-Poincare framework for semidirect products developed in [21[ to the case of an affine 
representation of the configuration space Lie group on the vector space whose dual are the 
convected variables. The resulting affine Euler-Poincare reduction is natural in two senses. 
Firstly, in the case of complex fluids and at the reduced level, it coincides with that given 
in 161]. Secondly, in the hyperregular case and through the Legendre transformation, it 
is compatible with the affine Lie-Poisson reduction for semidirect products discussed pre- 
viously. In addition, the Euler-Poincare formulation immediately leads to analogs of the 
Kelvin circulation theorem for complex fluids. 

Organization of the paper. The paper is organized as follows. The flrst third, comprising 
presents only theoretical results. These are developed and applied to a wide range of 
examples in the last section ^ forming the bulk of the paper. 

The theoretical part begins by recalling at the end of this introduction some needed 
facts about the Lagrangian and Hamiltonian reductions for semidirect products. It is well- 
known that these processes form the basic framework for the geometric formulation of various 
models of simple fluids. In ^ we present the theory of affine Euler-Poincare reduction for 
a general Lie group G acting by affine representation on a dual vector space V*. This 
theory is specialized, in §31 to the case of general complex fluids. More precisely, we describe 
concretely the group G and the cocycle needed in the affine representation, in order to obtain 
by reduction the general equations for complex fluids. In order to carry out the Hamiltonian 
side of the theory, we need to recall and state some results concerning the reduction of 
a canonical symplectic structure with a magnetic term. This is the subject of ^ whose 
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principal result states that, under some conditions, reducing a canonical symplectic form 
relative to a cotangent lift with an affine term is equivalent to reducing a magnetic symplectic 
form relative to the cotangent lift. This observation is used in ^ to obtain the theory of 
affine Lie-Poisson reduction for a general Lie group G acting by affine representation on a 
dual vector space V* . In particular, we compute the associated momentum map and Poisson 
bracket as well as the symplectic reduced spaces. This theory is also shown to be a particular 
case of the process of reduction by stages for nonequivariant momentum maps. In ^ we 
specialize these results to the case of the group and the cocycle involved in the description 
of complex fluid dynamics. In particular, we compute the associated Poisson bracket and 
momentum map. This part constitutes the Hamiltonian side of the theory developed in ^ 
The Kelvin-Noether theorem is a version of the classical Noether conservation law that holds 
for solutions of the Euler-Poincare equations. For example, an application of this theorem 
to the compressible fluid gives the Kelvin circulation theorem. The generalization of this 
result to the case of affine Euler-Poincare equations is the subject of §71 

The rest of the paper is devoted to applications dealing with spin systems and complex 
fluids, all of them contained in ^ In §8.1l we start with the example of spin systems, since it 
illustrates the applicability of our theory in a very simple situation that exhibits, nevertheless, 
some of the key difficulties of more complicated fluid models. We then treat, in §8.2j - §8.7[ 
the examples of Yang-Mills and Hall magnetohydrodynamics for fluids and superfluids as 
well as the HVBK dynamics for ^He with vortices. In each case, we formulate in detail the 
Lagrangian and Hamiltonian reduction processes as well as the associated Poisson bracket 
and circulation theorems. In order to compare fluids and superfluids from a Hamiltonian 
point of view, we present in §8. 81 a summary of the models treated so far. In §8.9l we study the 
example of the Volovik-Dotsenko spin glasses and try to understand, from a Hamiltonian 
point of view, the passage from a given spin system to its hydrodynamic analogue. This 
process allows us to understand mathematically the link between the two approaches to spin 
glasses appearing in the current literature. All the Hamiltonian structures obtained so far 
by reduction of the canonical structure, coincide with the ones obtained previously and by 
various different methods in a series of papers by Holm and Kupershmidt. 

In the case of microfluids, treated in §8.101 we determine the Lagrangian and Hamiltonian 
structure of three models proposed by Eringen, namely, the micropolar, the microstretch, and 
the micromorphic models, corresponding to three groups associated to the internal structure 
of the fluid particles. Since our theory applies to any group, we can obtain new models 
of microlfuids such as the anisotropic micropolar or anisotropic microstretch models. In 
the case of microstretch and micropolar fluids, we show that the internal degree of freedom 
can be modeled by the group of invertible or unit quaternions, respectively. In §8.111 we 
determine and compare the Lagrangian and Hamiltonian structures of three models of liquid 
crystals dynamics. We also treat the case of Eringen's polymeric liquid crystals and give 
some information concerning the new model of anisotropic micropolar liquid crystals. 

We close this introduction by recalling some needed facts about Euler-Poincare and Lie- 
Poisson reduction for semidirect products. 

Notations for semidirect products. In the Euler-Poincare reduction for semidirect prod- 
ucts (see |21|]) one is given a Lie group G and a right representation p : G ^ Aut(l^) of G 
on the vector space V . As a set, the semidirect product 5* = G (8)1^ is the Cartesian product 
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S = G X V whose group multiphcation is given by 

(9i,vi)(g2,V2) = (gig2,V2 + Pg^ivi)). 

The Lie algebra of S is the semidirect product Lie algebra, s = q(§)V, whose bracket has 
the expression 

a%,«l)(6,f^2) = [{Cl,Vl),{^2,V2)] - ([6,61,^^16 -f^26), 

where denotes the induced action of g on V, that is, 

d 



dt 



t=0 



Prom the expression for the Lie bracket, it follows that for (^, e s and (/x, a) e 5* we have 

ad(^ (//, a) = (ad| // + ^; o a, , 
where G V* and f o a G 0* are given by 

d 



dt 



t=o 



and where (•, •)g : g* x g — >^ M and (-, ■)y : V* x V ^ M. are the duality parings. 
Lagrangian semidirect product theory. 

• Assume that we have a function L : TG x y* — > R which is right G-invariant. 

• In particular, if gq G V*, define the Lagrangian L^q : TG — >• M by LaQ{vg) := L{vg, ao). 
Then La^ is right invariant under the lift to TG of the right action of Ga^ on G, where 
GoQ is the isotropy group of ao- 

• Right G-invariance of L permits us to define / : g x — >^ M by 

KTgRg-^{vg),pl{ao)) = L(vg,ao). 

• For a curve g{t) G G, let ^{t) := TRg(^iyi{g{t)) and define the curve a{t) as the unique 
solution of the following linear differential equation with time dependent coefficients 

d{t) = -a{t)i{t), 

with initial condition a(0) = oq. The solution can be written as a{t) — p*(t)(ao)- 

Theorem 1.1 With the preceding notations, the following are equivalent: 
(i) With Qq held fixed, Hamilton's variational principle 

"12 



Lao{g{t),g{t))dt = 0, 
holds, for variations 5g{t) of g{t) vanishing at the endpoints. 



(ii) g{t) satisfies the Euler- Lagrange equations for La^ on G. 

(iii) The constrained variational principle 

S [ l{^{t),a{t))dt = 0, 
Jti 

holds on g X V* , upon using variations of the form 

df] 

6^= — -[^,ri], 6a = -ar], 

where r]{t) G g vanishes at the endpoints. 

(iv) The Euler- Poincare equations hold on q x V* : 



d_5l_ _ _ SI 
dt6^~~^^6^^Ta 



ad^ — + — oa. (1.2) 



It is worth noting that there is a remarkable symmetry breaking from L : TG x V* ^ 
to Lag : TG M. Whereas the function L is right G-invariant, that is, L{vhg, p*g-i{a)) = 
L{vh,a) for all Vh G TG, a G V*, g E G, the Lagrangian La^ is only Gag-invariant, that 
is, La(^{vhk) = Lag{vh) for all Vh G TG, k G Ga^- Note, however, that is a Lagrangian 
function on TG, whereas L is not since it is not defined on a tangent bundle. The relationship 
between these two different invariance properties are dealt with in the previous theorem and 
the evolution a(t) = Pg(t)Cto or, equivalently, the equation d{t) = —a{t)^{t), a(0) = ao, 
that needs to be added to (11. 2p in order to get a complete system of equations for all the 
unknowns, is exactly due to this symmetry breaking from G to Gag- We shall observe the 
same phenomenon throughout the paper when working out the affine reduction theorem on 
both the Lagrangian and Hamiltonian side. 

Hamiltonian semidirect product theory. Let S := G(§)V he the semidirect product 
defined at the beginning of this section. The lift of right translation of S" on T*S' induces a 
right action on T*GxV*. Consider a Hamiltonian function H : T*GxV* M right invariant 
under the S'-action on T*G x V*. In particular, the function Ha^ := -ff|r*Gx{ao} • T*G — > M 
is invariant under the induced action of the isotropy subgroup Gao := {g E G \ p*ao = ao} 
for any ao G V*. The following theorem is an easy consequence of the semidirect product 



reduction theorem (see [29]) and the reduction by stages method (see |26l |) 



Theorem 1.2 For a(t) G T*^^^G and fi(t) := T* Rg(t){,Ci{t)) G 0*, the following are equiva- 
lent: 

(i) a{t) satisfies Hamilton's equations for Ha^ on T*G. 

(ii) The Lie-Poisson equation holds on s* : 

d . , _ , , / _ 5h 5h\ 

— (/i, a) = - adfsh sh\[fJ', a) = - ad^^ fi -\- — o a, a— , a(0) = ao 

ot l37r''ss"J Y <5m da op J 



6 



where s is the semidirect product Lie algebra 5 = q(^V . The associated Poisson bracket 
is the Lie-Poisson bracket on the semidirect product Lie algebra s* , that is, 



' 6a 6a 5/i 



As on the Lagrangian side, the evolution of the advected quantities is given by a{t) = Pg(f)(cto)- 

For example, one can start with a Lagrangian as at the beginning of this section, 
suppose that the Legendre transformation WLag is invertible and form the corresponding 
Hamiltonian = o IFL~^^, where i^a;, is the energy of L^q. Then the function H : 
T*G X ^ M so defined is S'-invariant and one can apply this theorem. At the level of the 
reduced space, to a reduced Lagrangian / : g x V^* ^ M we associate the reduced Hamiltonian 
/i : 0* (S) V* ^ M given by 

61 

/i(/i,a) := -/(^,a), = 77. 

Since 

6h ^ 6h 61 

— =^ and — = - — , 
o/x 6a 6a 

we see that the Lie-Poisson equations for h on s* are equivalent to the Euler-Poincare equa- 
tions (11.21) for / together with the advection equation a + a^ = 0. 

Links with the reduction by stages. Consider the semidirect product Lie group S = 
G(§)V acting by right translation on its cotangent bundle T*S. An equivariant momentum 
map relative to the canonical symplectic form is given by 

Jij(a/, {u,a)) = T*L(^f^u){af, {u,^)) = {T*Lf{af) + uoa,a). 

Since is a closed normal subgroup of S, it also acts on T*S and has a momentum map 
Jv:T*S ^V* given by 

3v{c(f, (m, a)) = a. 

Reducing T*S by V at the value a we get the first reduced space {T*S)a = Jy^{a)/V. The 
isotropy subgroup Ga, consisting of elements of G that leave the point a fixed, acts freely 
and properly on (T*S)a and has an induced equivariant momentum map 3a '■ (T*S)a — >■ 0^, 
where Qa is the Lie algebra of Ga- Reducing {T*S)a at the point fia '■= /^ISa, we get the 
second reduced space {(T*S)a)^^ = Ja^(/^a)/(^<jJw- 

Using the Semidirect Product Reduction ([28|) or the Reduction by Stages Theorem 
([2^)5 the two-stage reduced space {{T*S)a)^^ is symplectically diffeomorphic to the reduced 
space (T*S')(^_a) = Ji?^(/^? ct)/G^(At,a) obtained by reducing T*S by the whole group S at the 
point (/i, a) G s*. 

The first symplectic reduced space {(T*S)a,fla) is symplectically diffeomorphic to the 

canonical symplectic manifold (T*G,Q) and the second reduced space (^{(T*S)a)^^ , (^a)/^aj 

is symplectically diffeomorphic to the coadjoint orbit (C(^,a)5 i^(^(,a)) together with its orbit 
symplectic form. Note also that we can consider the right G-invariant Hamiltonian H : 
T*G X ^ M as being the Poisson reduction of a S'-invariant Hamiltonian H : T*S — > M 
by the normal subgroup {e} x V since {T*S)/{{e} x V^) = T*G x V* . Here the G-action 
on T*G x V* is induced by the lift of right translation of S. Theorem IL2I is then a trivial 
consequence of these observations. 
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2 Affine Lagrangian Semidirect Product Theory 

Consider the right contragredient representation of G on the vector space V*. We 
can form an affine right representation 9g{a) — p*_i(a) + c{g), where c G J-'{G,V*) is a 
contragredient representation valued right group one-cocycle, that is, it verifies the property 
c{fg) = Pl-i{c{f)) + c{g) for all f,geG. This implies that c(e) = and c{g-^) = -p*g{c{g)). 
Note that 

t=o 

and 

(a^ + dc{^),v)v = (dc^(v) -voa, C)g, 
where dc : Q ^ V* is defined by dc(^) := Tec{$,), and dc^ : F ^ is defined by 

• Assume that we have a function L : TG x V* ^ M. which is right G-invariant under 
the affine action 

K,a) ^ {ThRg{vh),eg{a)) = {ThRg{vh), p;-i{a) + c{g)). 

• In particular, if ^ ^*) define the Lagrangian : TG ^ M by La^Xvg) '■= L{vg, ao). 
Then La^ is right invariant under the lift to TG of the right action of on G, where 
G^P is the isotropy group of ao with respect to the affine action 9. 

• Right G-invariance of L permits us to define I : g xV* ^'Rhy 

l{TgRg-i{vg),9g-i{ao)) = L{vg,ao). 

• For a curve g{t) G G, let ^(t) := TRg(t)-i{g{t)) and define the curve a{t) as the unique 
solution of the following affine differential equation with time dependent coefficients 

a(i) = -a(t)e(t)-dc(C(t)), 

with initial condition a(0) = ao- The solution can be written as a{t) — 6'g(t)-i (ao). 

Theorem 2.1 With the preceding notations, the following are equivalent: 
(i) With ao held fixed, Hamilton's variational principle 



S Lao{9{t),m)dt = 0, (2.1) 



holds, for variations Sg{t) of g{t) vanishing at the endpoints. 
(ii) g{t) satisfies the Euler- Lagrange equations for Lag on G. 
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(iii) The constrained variational principle 



t2 



6 mt),a{t))dt = 0, (2.2) 



holds on Q X V* , upon using variations of the form 

S^ = — - [^,r]], 6a = -ar] - dc{r]), 

where r](t) E q vanishes at the endpoints. 
(iv) The affine Euler-Poincare equations hold on g x V* : 

Proof. The equivalence of (i) and (ii) is true in general. 

Next we show the equivalence of (iii) and (iv). Indeed, using the definitions, integrating 
by parts, and taking into account that ri{ti) = ri{t2) = 0, we compute the variation of the 
integral to be 

5 j\m,ci{t))dt = 1'^ ((f ''^^) + (^"'^)) 

^, ?7 - adg r/^ - (^ar] + dc(r/), dt 



t2 



t2 
tl 



d 51 ,,5l\ I 51 , r / \ s , , 



'Ul d5l ^,51 51 ^t(^^\ W 

and so the result follows. 

Finally we show that (i) and (iii) are equivalent. First note that the G-invariance of 
L : TGxV* — >• M and the definition of a{t) = Og(t)~^ (cto) imply that the integrands in (12. ip and 
(12. 2p are equal. It is known that all variations 5g{t) G TG of g{t) with fixed endpoints induce 
and are induced by variations 5^{t) G of ^(t) of the form 5S, = t) — r]\ with rjit) vanishing 
at the endpoints; the relation between 5g(t) and rj(t) is given by r]{t) = TRg(t)-i{5g{t)). See 
[3] for details. 

Thus, if (i) holds, we define rj{t) = TRg^t)-^{5git)) for a variation 5g{t) with fixed end- 
points. Then if we let 5^{t) = TRg^iyi{g{t)), we have 5^ = rj — [C,,ri]. In addition, the 
variation of a(t) = 6g(iyi{aQ) is 5a{t) = —a{t)r]{t) — dc(r/(t)). Conversely, if 5C, = fj — [C,,ri] 
with ri(t) vanishing at the endpoints, we define 5g(t) = T Rg(^t){jl{'^)) ■ This 5g{t) is the general 
variation of g{t) vanishing at the endpoints. From 5a{t) = —a{t)ri{t) — dc{r]{t)) it follows 
that the variation of 6g(^t){ci'{t)) = cto vanishes, which is consistent with the dependence of 
^ao only on 5((t),5((t). ■ 
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3 Lagrangian Approach to Continuum Theories of Per- 
fect Complex Fluids 

Recall that in the case of the motion of a fluid on an orientable manifold the configuration 
space is the group G = Diff(T') of all diffeomorphisms of V. In the case of incompressible 
fluids, one chooses the subgroup Diffvoi('P) of all volume preserving diffeomorphisms, with 
respect to a fixed volume form on V. Besides the diffeomorphism group, the other basic 
object is the vector space V* of advected quantities on which G acts by representations. 
Typical advected quantities are for example the mass density, the specific entropy or the 
magnetic field. One can obtain the fluid equations by choosing the appropriate Lagrangian 
and Hamiltonian functions and by applying the semidirect Euler-Poincare or Lie-Poisson 



reduction processes ( Theorem 1 1 . 1 1 and 1 1 . 21) . see [28|| and [21 . 

The goal of this section is to extend these formulations to the case of complex fluids. At 
the reduced level, the Euler-Poincare equations for complex fluids are given in [l^ (equa- 
tions (3.23), (3.24), (3.32) and (3.33)). The two key observations we make regarding these 
equations are the following. First, the two equations (3.23) and (3.24) suggest that the 
configuration manifold Diff(P) has to be enlarged to a bigger group G in order to contains 
variables involving the Lie group O of order parameters. Second the two advection equa- 
tions (3.32) and (3.33) suggest that there is a new advected quantity on which the group 
G acts by affine representation. Making use of these two observations, we construct below 
the appropriate configuration space and the appropriate affine action for the dynamics of 
complex fluids. By using the general process of affine Euler-Poincare reduction developed 



before (Theorem 12. ip . we get (a generalization of) the equations given in [16 . 

Here and in all examples that follow, there are fields different from the velocity field for 
which we shall never specify the boundary conditions. We make the general assumption, 
valid throughout the paper, that all integrations by parts have vanishing boundary terms, or 
that the problem has periodic boundary conditions (in which case P is a boundaryless three 
dimensional manifold). Of course if one would try to get an analytically rigorous result, the 
boundary conditions for all fields need to be carefully specified. 

The configuration manifold. Consider a finite dimensional Lie group O. In applications 
O will be called the order parameter Lie group. Recall that in the case of the motion of a 
fluid on an orientable manifold V, the configuration space is the group G = Diff (I^) of all 
diffeomorphisms of V. In the case of complex fluids, the basic idea is to enlarge this group to 
the semidirect product of groups G = Diff (P) (S)JF(X), O). Here J-'{T>, O) denotes the group 
of all mappings x defined on T> with values in the Lie group O of order parameters. The 
diffeomorphism group acts on JF(D, O) via the right action 

(r/, x) e Diff (D) X J^{V, 0)^xov^ ^FiV, O). 

Therefore, the group multiplication is given by 

x){^, = {v° V, ix ° v)-^)- 

Recall that the tangent space to Diff(D) at rj is 

T^Difr(P) = {u, : I? ^ TP I u,(x) G r^(.)I?}, 
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the tangent space to J-{T>, O) at x is 

A direct computation shows that the tangent map of right translation is 

ri?(^,^)(u^, v^) = (u^ o ^, TR^(u^ oip)). 

For simphcity we fix a volume form /i on V. Therefore we can identify the cotangent space 
T* Diff (D) with a space of one- forms over 77, that is, 

T; Difr(P) = {m,:V^ T*V \ m,{x) G T;,)P}. 

The cotangent space of J^{T>, O) at x is naturally given by 

Using these identifications, the cotangent map of right translation is computed to be 

r*-R(<p,^)(m^, K^) = J{ip~^) (m^ o v9"\T*i?^o<p-i(K^ o ip-^)) , 

where J{ip^^) is the Jacobian determinant of the diffeomorphism (^~^, and the corresponding 
cotangent lift, defined by Rj^^p-j '■= T*R(^^^^yi, is given by 

^fJ,V)('^'7'%) = ° ip,T*R^-l{K^ Oip)). 

The Lie algebra q of the semidirect product group is 

= X(P)(DJ^(r>,o), 

and the Lie bracket is computed to be 

ad(u,z.)(v, C) = (adu V, ad^ ( + di^ • v - dC • u), 

where aduV = — [u, v], ad^^C ^ is given by ajd^({x) :— ad.^^x)C{^): and du • v E 

J-{'D, 0) is given by di/ ■ \{x) := dz/(a;)(v(a;)). 

Using the previous identification of cotangent spaces, the dual Lie algebra q* can be 
identified with 

Q\V)(i)J^{V,o*) 

through the pairing 

((m, «;), (u, I/)) = / (m • u + K • I/) /X. 
Jv 

The dual map to ad(u,i/) is 

ad(„ (m, = (i^uHi + (div u)m + K ■ dz/, ad* K + div(u/t)) . (3.1) 

This formula needs some explanation. The symbol k ■ dv E 0,^(1)) denotes the one- form 
defined by 

K ■ diy{vx) := K{x){du{v:^)) 
11 



and ad* k G J^{V, o*) denotes the o*-valued mapping defined by 

ad* := ad*(^)(K(x)). 

The expression u/t denotes the 1-contravariant tensor field with values in o* defined by 

UK^ax) := ax(}i{x))K{x) G o*. 

Since uk is a generalization of the notion of a vector field, we denote by X(T', o*) the space 
of all 1-contravariant tensor fields with values in o*. In (13.11) . div(u) denotes the divergence 
of the vector field u with respect to the fixed volume form /x. Recall that it is defined by the 
condition 

(div u)/i = -Cu/W- 

This operator can be naturally extended to the space X(T', o*) as follows. For w G X(T', o*) 
we write w = WaS"" where (e") is a basis of o* and Wa G X{V). We define div : X(V, o*) 
J^{V, 0*) by the equahty 

divw := {div Wa)£°'- 

Note that if w = uk, then 

div(uK) = dn ■ u + (div u)k. 

The space of advected quantities. In physical applications, the affine representation 
space V^* of G = DiS(V) (§)J^(V, O) is a direct product © ■, where V* are subspaces of 
the space of all tensor fields on T> (possibly with values in a vector space). Moreover: 

• Vi is only acted upon by the component Diff(T') of G. 

• The action of G on V^* is affine, with the restriction that the affine term only depends 
on the second component J^{V, O) of G. 



In this way, we obtain the affine representation 

(a, 7) G V* ® V; ^ {ar], 7(r/, x) + C(x)) G V* © V*, 



(3.2) 



where 7(77, x) denotes the representation of {rj, x) G G on 7 G V2, and C G JF(JF(T', O), V2) 
satisfies the identity 

G((xo^)^)=C(x)(v^,^) + G(^). (3.3) 

Note that this is equivalent to say that the representation p and the affine term c of the 
previous section have the particular form 

7) = («^, 7(^, X)) and c{r], x) = (0, C{x)). 
The infinitesimal action of (u, z/) G 5 on 7 G Vg* induced by the representation of G on 



^*is 



d 

~ di 



t=o 



7(exp(tu),exp(tz/)) = ^ 



7(exp(tu), e)(e, exp(tz/)) 



t=Q 



7(exp(tu),e) + ^ 



7(e, exp(tz/)) =: 7U + ju. 



t=o 
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Therefore, for {v,w) G © V2 we have 



(f , w) o (a, 7) = {v o a + w Oi 7, w 02 7), 

where Oi and 02 are associated to the induced representations of the first and second com- 
ponent of G on On the right hand side, the diamond operation o is associated to the 
representation of Diff(T') on V{. The space V{ is naturally the dual of some space Vi of 
tensor fields on V. For example the {p, q) tensor fields are naturally in duality with the 
tensor fields. For a G and w G Vi, the duality pairing is given by 

(a,w) = {a- v)fi, 
Jv 

where • denotes the contraction of tensor fields. 

Since the affine cocycle has the particular form c{rj, x) = (0, we obtain that 

dJ{v,w) = (0,dC^H). 

For a Lagrangian / = Z(u,z/,a,7) : [X{V) (i)J^{V, o)](i)[V^ © V^] M, the affine Euler- 
Poincare equations (12. 3p become 




^ 61 ,51 61 61 61 

-£uT (divu)- — ■du + — oa + — Oi'y 

ou ou ou da 07 

^* 5^ ^. f 6l\ 61 ^61 
- ad* ^ div u— + — 02 7 - dC^ 



'^61' \ 6v J 6^ \6^ 

and the advection equations are 

d + au = 

7 + 7U + 7Z/ + dC{u) = 0. 



(3.4) 



(3.5) 



Remark on the duality pairing. We remind the reader that, for simplicity, we have 
assumed that V is orientable and we have also fixed a volume form n G f2"(P). Thus, using 
the pairing, the dual space of the Lie algebra X{V) of vector fields is identified with the 
one-forms Therefore, the functional derivative 61 /6u is a one- form on V. 

The natural dual of X{V) is the space © Ven{V) of one-form densities on V. In 

this general approach, valid also for non-orientable manifolds, the duality pairing of a © G 
n^{V) © Ven{V) with u reads 

/ {a-\i)uj. 
Jv 

Then the functional derivative 61/6vl is interpreted as a one-form density on T). The same 
remark applies for the dual of vector valued tensors on P, in particular for the Lie algebra 
^(P,o). 

We chose to work on oriented manifolds and to avoid the use of Q}{T>) © T>en{V) be- 
cause this general case would considerably complicate the writing of various formulas and 
equations, without adding any essential new information. We emphasize that all results 
of this paper are also valid for non-orientable manifolds by appealing to the usual exterior 
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differential calculus and the Stokes theorem for twisted A;-forms on T> (see [2|], Supplement 
7.2A for details). ♦ 

Basic example. Take V2 '■= 0), the space of all one-forms on T> with values in 0. 

This space is naturally the dual of the space V2 = X{T>, 0*) of contravariant tensor fields 
with values in 0*, the duality pairing being given, for 7 e VL^{V, 0) and w G 0*), by 

Jv 

where 7 ■ w denotes the contraction of tensors. 

We consider for (13.21) the affine representation defined by 

(a, 7) ^^ (ar/, Ad^-i 77*7 + x^^Tx), (3.6) 

where Adi^,-i 77*7 + x~^Tx is the o-valued one-form given by 

(Ad^-i 77*7 + x'^Tx) {v.^) := Ad^(^)-i(r7*7(w^)) + xix^^T^xivx) , 

for Vx G TxT>. One can check that 7(?7,x) ■= -^^x~^ ^7*7 ^ right representation of G on 
V2 and that C(x) = X~^Tx verifies the condition (13.31) . In fact, (13.61) corresponds to the 
action of the automorphism group of the trivial principal bundle O x V, on the space of 
connections. 

For this example we have 

7U = jGu7, 7z/=— ad,y7 and dC^u) = du, 

where ad,y7 G Q^{M, 0) denotes the one-form given by 

(ad^7) (v.^) := ad,,(^)(7(t;^)) = [z/(x), 7(^7^.)]^ 

and du G Q^{T>, 0) is given by dulvx) '■= Txi'[vx) G 0. 
A direct computation shows that 

If Oi 7 = (div w) ■ 7 — w ■ i_ d7 G fi^ {V) , 
u;o2 7 = - Tr(ad;u;) eJ'iV.o*), 
dC'^iw) = -divw G J^(P,o*), 



where Tr denotes the trace of the 0*- valued (1,1) tensor 



ad! w : T*V x TV 



[ax.Vx] 



ad. 



'w{ax)). 



In coordinates we have Tr(ad* w) = ad*, w^. Making use of these computations, the affine 
Euler-Poincare equations (13. 4p become 



d 51 n 51 ^51 51 ^ 51 

— — = -£my~ ~ (divu)— - — ■dv+ — oa 
01 5\i ou ou 5v da 



+ ( div — 

07 



7 



— ■ i.d7 

07 



(3.7) 



d_5I_ 



u- 



5l_ 



Tr I ad; 



5l_ 

^7 
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and the advection equations are 



a + au = 

7 + - ad,^ 7 + dz/ = 0. 



These are, up to sign conventions, the equations for complex fluids as given in [16 



We now rewrite these equations using the covariant differentiation associated to a con- 
nection. A one-form 7 G 0) can be considered as a connection one-form on the trivial 
principal (9-bundle O x V ^ V, namely. 



K,eh) e T,v X no ^ Adh-i(7(x)K) + m-i(eh)) e o. 



(3i 



The covariant differential associated to this principal connection will be denoted by d"^. 
Therefore, for a function u G J-^{T>, 0), we have 



d^z/(v) :=dz/(v) + [7(v),z/]. 
The covariant divergence of w G 0*) is the function 

div''' w := divw — Tr(ad* w) G J-'(T>, 0*) 
defined as minus the adjoint of the covariant differential, that is. 



{d'^u ■ w) II 



[v ■ div'^ w) II 



(3.9) 



(3.10) 



(3.11) 



for all V G J-^{T>^ 0). Note that the Lie derivative of 7 G Q}{T>^ 0) can be written as 

£,7(v) = d(7(u))(v) + i,d7(v) 

= d^(7(u))(v) - [7(v),7(u)] +d7^(u,v) - [7(u),7(v)] 
= d^(7(u))(v) + i„i?(v). 



(3.12) 



where 



B := d'^7 = d7 + [7,7], 



is the curvature of the connection induced by 7. 

Note also that, using covariant differentiation, we have 

1(7 Oi 7 = (div w) ■ 7 — w • i_ d7 = (div'^ w) ■ — w ■ \_B . 

Therefore, in terms of d'*', div''', and B = (Pj, the equations read 

„ SI 



d 61 



'u— (divu) 



61 



51 51 , 51 

— ■ dv + — oa 

ou du da 

dW^ ^ ] ■ 7 — — • i_ i? 

07 / 07 

d 51 51 ( 5l\ ^ 51 

— — = - ad^ div u— + div^ — , 

ot5v 5v \ J ^1 



(3.13) 
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and 

a + au = 



7 + d^(7(u)) + iufi + d^z/ = 0. 

The ^-representation. In this example we want to reformulate ttie reduction process as 
well as the equations of motion fl3.13p in terms of another set of variables, namely, (u, z/, a, 5) , 
where B = d'^7 = d7 + [7, 7] G Q^{T>, 0) is the curvature of 7, instead of (u, z/, a, 7). As we 
shall see below, with this choice of variables the action f l3.6p becomes linear instead of affine. 
We shall also assume that the Lagrangian L, and hence also /, depend on 7 only through 
B. To do this we shall use the standard Euler-Poincare reduction for semidirect products 
(Theorem 11.11) . 

Indeed, if 7' = Ad^-i 77*7 + x~'^Tx then we have d''' 7' = Ad^-i ri*d'^'~f. Thus the repre- 
sentation of DiS{V)(s)J^{V, O) on ® VL^{V, 0) is given by 

(a, 5) {a7],Ad^~ir]*B). 

The associated infinitesimal action of (u, z/) G XiV) (S)JF(P, 0) is 

(a, -B)(u, z/) = (au, -B(u, z/)) = (au, £uB — ad,^ B). 

The space Q'^{T>, 0) is, in a natural way, dual to the space ^2{T^, 0*) of 2-contravariant skew 
symmetric tensor fields with values in 0*. The duality pairing is given by contraction and 
integration with respect to the fixed volume form /i. More generally, we can consider the 
space ^lk{1^, 0*) of fc-contravariant skew symmetric tensor fields with values in 0* and we can 
define the divergence operators, div, div''' : f2fc(P, 0*) Qk~i{1^, 0*), to be minus the adjoint 
of the exterior derivatives d and d'^, respectively. For example, div'^ is defined on fifc(I', 0*) 
by 



{d^a-uj)fi = - {a-dw^uj)fi, (3.14) 
V Jv 

where a G Vt^^^iV, 0) and uj G VtkiV, 0*) Note that we have used the notations VtiiV, 0*) = 
X(P, 0*) and Vto{V, 0*) = J^{V, 0*). 

Using the duality pairing defined before, for {v,b) G Vi © ^2(1^, 0*) and (a, B) G V{ © 
fi^(P, 0), the diamond operation is given by 

(f , b) o{a,B) = {voa + b Oi B, 6 02 B), 

where 

boiB = (divb) ■i_B - b-i_dB e n\V) 

and 

6 02 5 = - Tr(adB b) = - ad^^^ b'^ G J^{V, 0*). 

The Euler-Poincare equation (11. 2p are in this case 

d 61 ^61 , , ,61 61 61 

— — = -i-uT (divu)- — ■dz/+ — oa 

ar ou ou ou 6u 6a 



(3.15) 



mTu = -^<Tu^jB'^^-^'[^^^jB 
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and the advection equations are 

a + au = 

B + £^B -ad^B = 0. 

We now show that when B is the curvature of the connection 7, that is, B = 
then the affine Euler-Poincare equations (13.131) imply the standard Euler Poincare equations 
(I3.15p . Define the map 

^ ■.n\V,o) ^ n'^{V,o) by $(7) = d^7 

and suppose that the Lagrangians /^(u, u, a, 7) and /^(u, z/, a, B) are related by 

/^(u, a, $(7)) = /^(u, z/, a, 7). 

We have T^$(a) = da + [7, a] + [a, 7] = d^a and $(Ad^-i 77*7 + = Ad^-i r7*$(7). 

Therefore, by taking the derivative, the infinitesimal action verifies 

Bvl + Bu = d'''(7u + 7z/ + du), 

or, explicitly, 

£u5-ad^5 = d^(£u7-ad,.7 + dz/). (3.16) 

This proves that the advection equations for 7 implies the one for B. Using the preceding 
formula and the definition of the diamond operation we find that 

(div'' 6) Oi 7 = — 6 Oi 5 and (div''' 6) 02 7 + div(div'^ 6) = — 6 02 5, 

or, explicitly, 

div^ ((div^ 6)7) - div^ b ■ i_B = - divb ■ i_B + b ■ i_dB 

and 

div^(div^6) = Tr(ad^6). 



Using the equality 



div'^ ( o $ = 



\5B J 57 



we obtain 

7 = 777^1^ and ^027 + div — =— 02fi. (3.17) 
^7 oB d'j V^7/ 

This proves that the affine Euler-Poincare equations (I3.13P imply the standard Euler-Poincare 
equations (13.151) . 



4 Affine and Magnetic Lie-Poisson Reduction 

The goal of this section is to carry out a generalization of the standard process of Lie- 
Poisson reduction for Lie groups, which is motivated by the example of complex fiuids. The 
only modification lies in the fact that the Lie group G acts on its cotangent bundle by a 
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cotangent lift plus an affine term. The principal result of this section states that, under some 
conditions, reducing a canonical symplectic form relative to a cotangent lift with an affine 
term is equivalent to reduce a magnetic symplectic form relative to the right- cotangent lift. 
At the reduced level, we obtain affine Lie-Poisson brackets and affine coadjoint orbits, whose 
affine terms depend on the affine term in the action. 

Consider the cotangent lift of the right translation Rg on a Lie group G. Recall that 
is the right action of G on T*G given by 

Rf{af)=T*Rg-.{af). 

Consider the map : T*G T*G defined by 

mg{af):=Rl\af)+Gg{f), (4.1) 

where G : G xG ^ T*G is a smooth map such that Gg{f) G T^gG, for all f,gEG. The map 
is seen here as a modification of the cotangent lift by an affine term G. The following 
lemma gives the conditions guaranteeing that the map is a right action. 

Lemma 4.1 Consider the map "^g defined in (14. ip . The following are equivalent. 
(i) "^g is a right action. 

(ii) For all f,g,h& G, the affine term G verifies the property 

Gg,{f)=G,{fg) + Rr{Gg{f)). (4.2) 

(iii) There exists a one-form a G Q^{G) such that Gg{f) = a{fg) — R^ {c({f)). 

Proof. The equivalence between (i) and (ii) is a direct computation. Suppose that (ii) holds. 
By setting / = e in formula (14. 2p . we obtain that Ggh{e) = Gh{g) + R^* {Gg{e)). Therefore 
we can define the one-form a by a{g) := Cg(e), and we have Gh{g) = d^gh) — Rj^*{a{g)). 
Conversely, suppose that (iii) holds. Then a direct computation shows that ( 14. 2p holds. ■ 

We denote by C{G) the space of all maps G : G x G ^ T*G, {g, f) ^ Gg{f) G T/^G 
verifying the property (14.21) . Remark that given an affine term G G C{G), the one-form a 
in item iii, is only determined up to a right-invariant one-form. Denoting by f2^(G) the 
space of all right-invariant one-forms on G, we have an isomorphism between C{G) and 
Q-^ {G) / Q]^{G) . This space is clearly isomorphic to the space ^1{G) of all one-forms a on G 
such that a{e) = 0. From now on, when we say that the one-form a is associated to G, we 
shall always assume that a{e) = which then guarantees the uniqueness of this one- form. 

In order to carry out the symplectic reduction associated to an affine action of the 
form (14. ip . we make two crucial observations (see Theorem 14. 6p . 

• Let a G ^o{G) be the one-form associated to and consider the associated fiber 
translation ta on T*G defined by 

t^i(3j) ■.= (3f-aif). 

Then ta, viewed as a map from the canonical cotangent bundle {T*G, Qcan) to the 
magnetic cotangent bundle (T*G,flcan — ^r^da), is a symplectic map. Moreover, ta is 
equivariant with respect to the action \E'g on (T*G, Qca.n) and the cotangent lift on 
{T*G,n-n*ada). 
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• Suppose that da is G-invariant. Then the action '^g is symplectic relative to the 
canonical symplectic form Qcan- 

From these observations we conclude that, under some conditions to be specified later on, 
reducing the canonical cotangent bundle {T*G, fican) relative to an affine action of the form 
(14. ip is equivalent to reduce the magnetic cotangent bundle (T*G, flcan — vr^da) relative to 
the cotangent lift of right translations. It is therefore useful to recall below some facts about 
the reduction of magnetic cotangent bundle. 

Some facts about magnetic cotangent bundle reduction. We first recall the following 

result (Theorem 7.1.1 in ^]) about the existence of momentum maps associated to cotangent 
bundles with magnetic terms. 

Theorem 4.2 Let B be a closed two-form on a connected configuration manifold Q. Let 
^ : G X Q ^ Q be a free and proper right action leaving the form B invariant. Consider the 
cotangent lift $J* of the G-action $g to the symplectic manifold {T*Q, Qca,n — t^q^), where 
fican is the canonical symplectic form and ttq : T*Q Q is the cotangent bundle projection. 
Suppose that there is a smooth map (p '■ Q —>■ g* that satisfies 

ig^i3 = d(0,O 

for all C, € 0. Then the following hold: 

(i) The map J = Jcan~0°7rQ, where Jean is the standard momentum map for the G-action 
relative to the canonical symplectic form, is a momentum map for the cotangent lifted 
action of G on T*Q with symplectic form Qcan — t^q^- 

(ii) The momentum map is, in general, not equivariant. Its non-equivariance Q*-valued 
group one-cocycle a : G ^ g* is given by 

a((7) = -0(<|.,(g))+Ad;(0(g)), (4.3) 

with the right hand side independent of qEQ. The cocycle identity is 

a{gh) = Ml{a{g)) + a{h). 

Note that the G-invariance of B ensures that the action is symplectic relative to the 
symplectic form fican — '^*q^- The momentum map J : T*Q q* is equivariant relative to 
the right affine action of G on q* given by 

OgiX) ■.= Ml\ + a{g). 

We turn now to the particular case when the configuration manifold is a Lie group G 
and the action is the cotangent lift of right translation, that is, $ J* = . We first recall 
from [26] (Theorem 7.2.1) the magnetic Lie-Poisson reduction theorem. 

Theorem 4.3 Consider a closed G-invariant two-form B on G. The Poisson reduced space 
for the right cotangent lifted action of G on (T*G, Qcan — t^g^) ^■^ S* with Poisson bracket 
given by 

for f,genQ*). 
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5/ Sg_ 



i5(e) ( ^ ^ 



Note that for a Hamiltonian /i : g* — M, the corresponding Hamiltonian vector field is given 
by 

(Sh 

When the magnetic term B is absent, we recover the standard Lie-Poisson bracket on g*. 
In this particular case, the symplectic reduced spaces J^^(/i)/G^ are symplectically diffeo- 
morphic to the coadjoint orbits = {Ad*(/i) \ g E G}, where Jr is the standard right 
momentum map for the G-action relative to the canonical symplectic form and is the 
coadjoint isotropy group of /i G g*. The symplectic diffeomorphism is induced by the map 

if : J^^(/i) Of,, ifiag) = T*Rg{ag) = Ad*-i /i. 

When a magnetic term B is present, the existence of a momentum map is not guaranteed, 
and we shall use the result of Theorem 14.21 Recall that the infinitesimal generator of ^ for 
the right action is the left invariant extension of ^; that is, C,g = Suppose that there is a 
smooth map : G — g* that satisfies 

^LB = d{^,0 

for all ^ G g. Since is determined by this equation only up to a constant, we can always 
impose the condition 0(e) = 0. From Theorem 14.21 the map J = Jr — 0o7rG is a momentum 
map for the cotangent lifted action of G on T*G with symplectic form Qcan — ^h^- -'■^ ^^^^ 
case, using the relation fl4.3l) . the nonequivariance cocycle is simply given by cr = —0. We 
denote by G^ the isotropy group of /i relative to the affine action 6g{ii) = Ad*Jfi) + a{g). The 



;rop; 

following result (see 26|, Theorem 7.2.2) shows that the reduction of a magnetic cotangent 
bundle of a Lie group at a given point /i G g* is symplectically diffeomorphic with the affine 
coadjoint orbit (9^ passing through fi. 

Theorem 4.4 Consider a dosed G-invariant two-form B on G, suppose that there is a 
smooth map : G — g* that satisfies 

for all ^ E Q, and consider the momentum map J = Jean ~ ° ttq. Then for each /i G g*, the 
symplectic reduced space 3^^{^)/G'^ is symplectically diffeomorphic to 

0;^={eg{^^)=Ml^i + a{g)\geG], 

the affine orbit through fi. The tangent space at X = OgifJ^) G (9^ to is given by 

r,0;= {ad^ A -S(^,-) I ^Gg}, 

where S(^, ■) := —Tf.a{C,) = —B{e){^, ■). The symplectic structure on has the expression 

to+{\) (ad*A -S(e,-),ad;A-S(r7,-)) 
= (A, %ri])-m.v). 

which we call the magnetic orbit symplectic form. 



20 



The symplectic diffeomorphism between the symplectic reduced spaces and the affine 
coadjoint orbits is constructed as follows. Consider the smooth map : J~^(/i) — > 
defined for ag G J~^(/i) by 

(p{ag) := 6'g-i(/i) = Ad*-i /i + (r{g~^). 
Then ip is G^-invariant and induces a symplectic diffeomorphism 

Note that we have (p{ag) = T*Rg{ag). Indeed, since ag G J^^{fj,), we have = T*Lg{ag) — 
4>{g), therefore we obtain 

= T*Rg{ag)+Ad;-,{a{g)) + a{g-') 
= T*Rg{ag). 

The general theory of symplectic reduction implies that the affine coadjoint orbits (C^jCU^) 
are the symplectic leaves of the Poisson manifold (g*, { , where {,}b denotes the Poisson 
bracket in Theorem 14. 31 

The following proposition shows that when the magnetic term B is an exact two-form, 
then the magnetic cotangent bundle is symplectomorphic to the canonical symplectic cotan- 
gent bundle. Through this symplectomorphism, the cotangent lift is transformed into an 
affine action. In the particular case of a Lie group G acting on its cotangent bundle by right 
cotangent lift, this affine action is of the form (14.11) . 

Proposition 4.5 Assume that all the hypotheses of Theorem \4.S\ are satisfied and suppose 
that B = da. Let t^ '■ {T*Q, ^can) ~^ {T*Q, ^can — '^q^) ^he fiber translation defined by 

taiPg) ■= Pq " 

Then the following hold: 

(i) ta is a symplectic fiber translation. 

(ii) The symplectic action \I' on (T*Q, Qcan) induced by through the map ta, that is, 
^9 := t~^ o $ o for any g & G, is the affine action given by 

^g{(3g) = ^Yif3g) + Cgiq), whcrc ( g) i = « ( 5, ( ) - («( ^ ) ) . 
(iii) A momentum map relative to the G-action \I' on (T*Q, Qcan) is given by 

Ja = J ota = Jean ° - (p O TCq] 

its nonequivariance cocycle equals the nonequivariance cocycle of 3 . 
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Proof, (i) See Proposition 6.6.2 in 27 
(ii) We have 



= $f (/?,) + a($,(g))-<l>f(a(g)) 
= $f(/3,) + C,(g). 

(iii) This is a consequence of the fact that Jq, = J o t^- The nonequivariance one-cocycle of 
3a is 

= crig). m 

Note that da is assumed to be G- invariant. When the one- form a is also G-invariant, then 
the affine term C vanishes. Therefore, the most interesting case happens when a is not 
G-invariant but da is. 

AfRne Lie-Poisson reduction. We apply now the previous results concerning the reduc- 
tion with magnetic terms to our initial problem, that is, the Lie-Poisson reduction of the 
canonical cotangent bundle (T*G, fican) with respect to an affine action of the form (14. ip . 
We obtain below, as an easy consequence of the previous theorems, the main result of this 
section. 



Theorem 4.6 Consider the symplectic manifold {T*G, Qcan), o,nd the affine action 

where C G C{G). Let a G ^^o(^) one-form associated to ^ g. Then the following hold: 

(i) The fiber translation t^ : {T*G, fican) {T*G, ^7can ~ vr^da) is a symplectic map. The 
action induced by \E'c, on (T*G, flca.n — vr^da) through t^ is simply the cotangent lift 
Rf. 

(ii) Suppose that da is G-invariant. Then the action is symplectic relative to the 
canonical symplectic form flca.n- 

(iii) Suppose that there is a smooth map (p : G ^ q* that satisfies 

i^Lda = d(0, ^) 

for all^ G 0. Then the map Jq, = 3 Rota — (f>onG is a momentum map for the action 
relative to the canonical symplectic form. We can always choose such that 0(e) = 0. 
In this case, the nonequivariance one-cocycle of 3a is a = —0. 
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(iv) The symplectic reduced space /C^^VL^) is symplectically diffeomorphic to the 

affine coadjoint orbit {O^^uo^), the symplectic diffeomorphism being induced by the 
Gj^ -invariant smooth map 



an 



Proof, (i) That ta is a symplectic map follows from item (i) in Proposition 14. 51 From item 
ii in Proposition 14.51 we know that the action induced on {T*G, fican) by the right cotangent 
lift on {T*G, fican ~ ^gda) through the map ta is the affine action whose affine term is given 
by a. This is precisely the action "^g. Thus we conclude that induces the right cotangent 
lifted action on {T*G, Qcan — ttqcIq;). 

(ii) When da is G-invariant, we know that the right cotangent lift is symplectic relative 
to i7 — TTgda. Since t^ is symplectic, we conclude the result. 

(ill) This follows from item i and ii in Theorem 14.21 Since (p can be chosen modulo a 
constant term, we can impose the condition 0(e) = 0. From the relation (14. 3 p we obtain the 
equality = — cr. 

(iv) The symplectic diffeomorphism t^ induces a symplectic diffeomorphism between 
the reduced spaces. Therefore, by Theorem 14.41 we obtain that the point reduced space 
(J^^(/i)/G'^, fi^) is symplectically diffeomorphic to the affine coadjoint orbit (C^,ciJ^), the 
symplectic diffeomorphism being induced by the map ip := ip o t^ : Jq^(/u) — ^ O'p^. We have 

= ^{t.%)) = T*RgiPg - a{g)) = R^M^g) - R'^'Mig)) 



Ri:,%)+Cg-.{g) = ^g{P, 



9/1 



where in the fourth equality we used Cg{f) = a{fg) — Rg {a{f)). ■ 

The affine coadjoint orbits (^0'^,lj^) are symplectic leaves in the affine Lie-Poisson space 
(0*,{,}^), where 

'SI sg; 

Note that this bracket is the Lie-Poisson bracket on the Poisson submanifold g* x {1} c 
0* := g* X M, where g is the one-dimensional central extension of g defined by the cocycle 
-S. 



{f,g}t{fi) = {fi 



Assume that — S integrates to a group two-cocycle B : G x G 



that is. 



-s(e,r7) 



dtds 



{B{g{t),h{s))-B{h{s),gm 



t=s=0 



where t 



g{t) and s ^— h{s) are smooth curves through e E G with tangent vectors 



dg{t) 

dt 



t=0 



and rj 



dhjs) 
ds 



. Let G be the central extension of G defined by the two- 



cocycle B and recall that the Lie algebra of G equals g. Then the affine coadjoint orbit (9^ 

is obtained by usual Lie-Poisson reduction of T*G relative to the lift of right translation at 
(/i, 1). See |26|, §6.2 for more details. 
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5 Affine Hamiltonian Semidirect Product Theory 

This is the Hamiltonian version of Section O More precisely, we carry out the Poisson and 
symplectic reductions of a canonical cotangent bundle (T*S,flcan), where S = G (§)V is the 
semidirect product of a Lie group G and a vector space V and where S acts on its cotangent 
bundle by cotangent lift plus an affine term. We will see that this process is a particular 
case of the theory developed in the previous section. 

Consider the semidirect product Lie group S := G (§)V associated to a right representa- 
tion p : G ^ Aut(l^). The cotangent lift of the right translation is given by 

R{g,v)iaf, {u,a)) : = T(},,)(g„)%,^)-i(a/, {u,a)) 

We modify this cotangent lifted action by an affine term of the form 

C{g,v){f,u) := {Ofg,v + pg{u),c{g)), (5.1) 

for a group one- CO cycle c G J^{G, y*), that is, verifying the property c(/(?) = p*„i(c(/))+c((?), 
as in Section [21 The resulting affine right action on T*S is therefore given by 

^(g,v){af,{u,a)) : = (m, a)) + m) 

= {RY{af),v + Pg{u),p*g.,{a) + c{g)) (5.2) 

This action is clearly of the form (14.11) . We now check that property (14.21) holds. This will 
prove that "^{g^v) is a right action. Indeed, 

((/, u) {g, v)) + Rfl^) (C(<;,„) (/, u)) 

= C(^h,n,){fg, V + Pg{u)) + Rfh,nj)i^fg, V + Pg{u), c{g)) 

= (0/9/1, w + phiy + pgiu)), c{h)) + {Rl''{Ofg),w + ph{v + Pg{u)), pl-^{c{g))) 
= (0/g/,, w + phiv + Pgiu)), c{h) + pl^i{c{g))) 
= (Ofgh^w + Phiv) + pghiu),cigh)) 

= G{gh,W + Ph{v))if,U) = G(^g^y)(h,w)if,U). 

In the following lemmas, we compute the one-form a G ^^o(>S') associated to G and we 
show that it verifies the hypotheses of Theorem 14.61 Recall that a is defined by aig,v) : = 
C{g,v)ie,0). 

Lemma 5.1 The one-form a G ^liS) associated to the affine term (15. ip is given by 

aig,v)i^g,iv,u)) = {cig),u), (5.3) 

for i^g, iv,u)) G T(^g^y)S. Moreover B := da is S-invariant and its value at the identity is 
given by 

B(e,0)((e,M), ir^,w)) = (dc(0,w) - (dc(r/),M). (5.4) 
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Proof. For {^g, {v,u)) G T(^g,v)S we have a{g,v) := C(g_t,)(e,0) = {Og,v,c{g)). Therefore we 
obtain the equahty a{g,v){^g, {v,u)) = {c{g),u). We now prove the right-invariance of da. 
For {^,u) G s define the associated right-invariant vector field {^,u)^ G X{S) by 

{^,u)\g,v) :=Ti?(,,,)(e,^/) = {e{9),v, Pg{u)), 

where ^^(^f) := T^RgiO- Note that we have a (5-,^^) = {c{g),Pg{u)) = -{c{g-^),u). 

Therefore 

^(g,^)da(e,0)((C,M), {v,w)) = da{g,v){TR(^g^^^{^,u),TR^g^^){r],w)) 
= da{g, v) ((^, u)^{g, v), [t], w f{g, v)) 
= d{a{{v,wf)){g,v) {{r],w)\g,v)) 

- d (« ((r/, w)^)) {g, v) iiv, w)\g, v))-a{ [(e, uf, (r/, wf] ) {g, v) 

d 



dt 



a {{r], w)^) {exp{t^)g, v + pg{ut)) 



d 
' dt 



" u)^) (exp(t?7)fi', V + pgiwt)) 



t=o 



+ a {{[i,rj\,uri - wi)^) {g,v) 



d 
dt 



[c{g-'exp{tO-'),w) + 



t=o 



d 
dt 



[c{g ^exp{tr]) ^),u) 



t=o 



+ oiig, v) ([^, ri]^{g),pg{uri - w^)) 
{c{g~^)^ + dc(0, - {c{g'^)ri + dc(r/), m) - {c{g~^),ur] - 
(dc(0, w) - {dc{r]), u) = da{e, 0)((^, u), {t], w)), 



where in the fourth equahty we used the identity 
and for the sixth equahty the identity 



R 



d 
di 



[cig-'expitO-'),w)= 



{p:^^anicig-')) + ciexpitO-'),w) 



-{c{g-')^ + dc{0,w). 



This lemma shows that hypothesis (ii) of Theorem 14.61 is verified. This implies that the 
action "^(g^v) is symplectic with respect to the canonical symplectic form Qcan on T*S. We 
now check hypothesis (iii) of Theorem 14. 6[ This implies that the action "^(g^v) admits a 
momentum map relative to the canonical symplectic form. 

Lemma 5.2 The map (p : S ^ s* defined by 

(f>{9, = {dJ{v) -vo c{g), -c{g)), 

verifies the property 

i{C,«)ida = d{(j), 

where {^,u)^ G X{S) is the left-invariant vector field induced by {C,,u) G s. 
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Proof. We will use the following formulas: 



{<P{g, v), (e, u)) = -{c{g),u) + {v, c{g)i + dc(0), 

r%,,)-i {{i,u)'^{g,v)) = (Ad,^,p,-i(M + 0). 

By the right-invariance of da, we have 

kt,,u)Lda{g,v) {{r],w)^{g,v)) = da{g,v) {{^,u)^{g,v), {r],w)^{g,v)) 
= da{e, 0) ((Adg ^, p^-i {u + <)), (Ad^ r], pg-i {w + vrj))) 
= {dc{Adg^),Pg~i{w + vri)) - {dc{Adgr]), pg-i{u + v^)) 
= {Tgc{TLg{0), w + vf])- {Tgc{TLg{r])), u + vO. 



For the last equality, we use that 

pl-.{dc{AdgO) = P;-rj^ 



c{gexp{t^)g ^] 



t=Q 



d 
It 



(c(^exp(tO) - c{g)) 



= TgC{TLg{i)). 

On the other hand we have 

d{<P,{i,u)){g,v) {iv,wng,v))= 



{(j){g exp{tri),wt + Pexp(t^)(^^)), u)) 



t=0 



d_ 
It 



{c{gexp{tri)),u) + (wt + Pexp(t^)(f ), c(^exp(t?7))^ + dc(0) 



-{Tgc{TLg{r])), u) + {w + vr], c{g)^ + dc(0) + 



d 
dt 



{v,c{gexp{tr]j)^) 



t=o 



-{Tgc{TLgir])),u) + {w + vrj,TgciTLg{m - {vi,Tgc{TLg{r^))) 
{Tgc{TLg{i)),w + vr^)- {Tgc{TLg{r,)),u + <)• 



Thus we obtain that 



\^,u)Lda = d{(f), 



The momentum map. By item (iii) of Theorem 14.61 and using a G ^o{S) given by (I5.3p . 
we obtain that a momentum map for the right-action 



'^ig,v){af,{u,a)) = iaf),v + pg{u), p* ,{a) + c{g)) 



is given by 



JaiPf, {u, a)) = 3R{ta{pf, {u, a))) - u) 
= Jr(/?/, (m, a - c(/))) - u) 
= ^*^(/,n) {u, a - c(/))) - u) 

= {T*Lj{(3f) +uo{a- c{f)), a - c{f)) - {dJ{u) -uo c{f), -c(/)) 

= {T*Lf{Pf)+uoa-dc^{u),a), (5.5) 
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with nonequivariance one-cocycle 

u) = u) = {uo c(/) - dc'^H, c(/)) e 5*. (5.6) 

Poisson bracket and Hamiltonian vector fields. Using Theorem 14. 31 and the expression 
of B given in Lemma 15. ![ we obtain that the reduced Poisson bracket on s* is given by 

'^1 



+ B(e 0) ( ( — — ^ ^ 



+ ( dc 



- ( dc 



59\ Sf 



SI Sg_ 

Given a Hamiltonian function /i : 3* — M, the corresponding Hamiltonian vector field with 
respect to the bracket { , }b, is given by 



5h 
6a 



— adih, /i — — o a + dc 

6fj. " 



r6h 


6h' 






5a ^ 






, —a 


5h 




6fi 



6h 
6fi 



The symplectic reduced spaces. By item iv of Theorem 14.61 the reduced space 

(j^^(/i,a)/^(;^„),fi(^,a)) 



is symplectically diffeomorphic to the affine coadjoint orbit yO^^^^u^j. More precisely, we 
have 



Oi^,a) = {Ad^^^„)(/i,a) + a((7,n)| {g,u) G S} 

= {{Adln + uoOgia) -dJ{u),eg{a))\{g,u) G S} 



(5.7) 



where 6g{a) := p*_i(a) + c{g) is the affine action of G on V*. Note also that the bilinear 
form S appearing in the formula of the affine orbit symplectic form (see Theorem 14. 4p , is 
given by 



U),-) = -T(e,0)O-(^, U) = - 

= (dc^(«),-dc(0), 



d 

dt 



(tu o c(exp(t,^)) — dc'^(tu), c(exp(t^))) 



t=o 



where (^,m) G s. The tangent space to the affine coadjoint orbit O^^^^ at (A, 6) is equal to 

{ad^^,„)(A,6)-S((e,n),-)|(e,«)G5} 
{ (ad^ X + uob- dc^(M), + dc(0) | (^, m) G 5} . 



^(A,fe)C{'^,a) 
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The symplectic form on O^^^ is given by 

u+{\, h) |^(ad^ \ + uoh- dc^(M), + dc(0) , 

(ad* \ + w oh — d(F{w), hq + dc^rf))^ 
= ((A,6),[(e,w),(r/,«;)])-S((e,w),(r7,«;)) 

= (A, % r,]) + (6, ur, - <) + (dc(r/), - (dc(0, w). (5.8) 
Recall that the affine coadjoint orbits are the symplectic leaves of the Poisson manifold 

We now apply these results to the main result of this section, that is, the Hamiltonian 
counterpart of Theorem 12.11 

Consider a Hamiltonian function H : T*G x V^* ^ M right-invariant under the G-action 

K,a) ^ K),e,(a)) := K), p*-i (a) + c((?)). (5.9) 

This G-action on T*G x V* is induced by the ^-action (15.21) on T*S given by 

'^{g,v){ah, {u,a)) = R[g]^^{ah, {u,a)) + C(^g^y){h,u) 

= {R^* {ah),v + pg{u), pl-i{a) + c{g)) . 

Note also that we can think of this Hamiltonian H : T*G x V* ^ W as being the Poisson 
reduction of a ^-invariant Hamiltonian H : T*S ^ M by the normal subgroup {e} x V since 
{T*S)/{{e} X V)^T*G x V*. 

In particular, the function Ha^ := -f^|T*Gx{ao} • T*G — > M is invariant under the induced 
action of the isotropy subgroup G^^ of oq relative to the affine action 9, for any oq G V*. 
Recall that Og{a) := p*_i(a) + c{g) for any g E G and a G V*. The following theorem is a 
generalization of Theorem 11.21 and is also a consequence of the reduction by stages method 
for nonequivariant momentum maps, together with the results obtained in section H] and at 
the beginning of the present section. 

Theorem 5.3 For a(t) G T*fj_-^G and /i(t) := T*Rgi^t){,Ci{t)) G 0*, the following are equiva- 
lent: 

(i) a{t) satisfies Hamilton's equations for Ha^ on T*G. 

(ii) The following affine Lie-Poisson equation holds on 5* : 

d , , f ^ T f6h\ 6h , f6h 
— {fi,a)= - adsh fi- — oa + dc — , -a- de - 
ar \ *M da \da J o/i \djj 

a(0) = Oq, where s is the semidirect product Lie algebra s = 0(§)V. The associated 
Poisson bracket is the following affine Lie-Poisson bracket on the dual s* 



{f,9}B{f^,a) = (/i 



51 5g_ 



' 6a 6fi 5a 6p 



dc{^-l) ^-l\-/dc(^) ^ 
5n J ' 5a / \ \5n J ' 6a 
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As on the Lagrangian side, the evolution of the advected quantities is given by a{t) = 

Proof. Recall that the momentum map relative to the canonical symplectic form on T*S 
and to the action 

^ig,v){Pf,{u,a)) = {R^'{Pf),v + Pg{u),p;-^ia)+c{g)) 

is given by 

JaiPf, (u, 0.)) = iT*Lf{(3f) +UO a — dc^{u), a). 
The action '^{g,v) of S induces an action of V given by 

(/?/,(«, a)) ^ {(3f,{v + u,a)). 

Since is a closed normal subgroup of S", this action admits a momentum map given by 

Jv(/3/, (m, a)) = a. 

Since V is an Abelian group, the coadjoint isotropy group of a & V* is Va = V and the 
first reduced space (T*S)a^ = 3~^{a)/V is symplectically diffeomorphic to the canonical 
symplectic manifold (T*G,flca.n)- The action \E' of S" on T*S restricts to an action of 
Gl(S)V on Jy^(a), where 

G: := {geG\pl-,{a)+c{g)=a}. 

Passing to the quotient spaces, this action induces an action of G^ on (T*S)a, which is readily 
seen to be the cotangent lifted action of G^ on T*G. We denote by : {T*S)a {qD* the 
associated equivariant momentum map, where gl is the Lie algebra of G^. Reducing {T*S)a 
at the point fia '■= /^Ifla; we get the second reduced space {{T*S)a)^^ = Ja^{fJ^a)/{G^)^^, with 
symplectic form denoted by {Qa)fj.a- 

By the Reduction by Stages Theorem for nonequivariant momentum maps, the second 
reduced space is symplectically diffeomorphic to the reduced space 

obtained by reducing T*S by the whole group 5* at the point {p, a) G s*. Here S"^ ^-^ denotes 
the isotropy group of the afiine coadjoint action with cocycle a given in (15.61) . 

As shown at the beginning of this section, this reduced space is symplectically diffeo- 
morphic to the affine coadjoint orbit (^O^^ ^-^,00^^ endowed with the affine orbit symplectic 

symplectic form described in fl5.8l) . These afiine coadjoint orbits are the symplectic leaves of 
the Poisson manifold (5*, { , 

Note finally that we can consider the right-invariant Hamiltonian H : T*G x V^* ^ M, as 
coming from a S*- invariant Hamiltonian H : T* S ^ M.. 

The theorem is then a consequence of all these observations. ■ 

Reconstruction of dynamics. We give now some details concerning the passage from the 
reduced formulation (ii) to the canonical formulation (i). Let {fi{t),a{t)) ^ q* x V* be the 
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solution of the affine Lie-Poisson equations, with initial condition (/io,ao). Then the curve 
a(t) := T*Rg(^tyi{fi(t)), where g(t) G G satisfies the linear ordinary differential equation 
with time-dependent coefficients 



g{t)=TR,(t)(jl^y 9(0) 



is the solution of Hamilton's equations associated to ifao on T*G and with initial condition 
a(0) = fiQ. Note that for the curve g{t) define above, we have a{t) = ^g(t)~i(ao)- 

The preceding theorem is compatible with Theorem I2.3[ Indeed, we can start with a 
Lagrangian Lag : TO — > M as in Section[2|, that is, we have a function L : TG xV* —^M which 
is right G-invariant under the affine action {vh, a) (ThRg{vh), Og{a)) = {ThRg{vh), p*g-i{a) + 
c{g)), such that La^ivg) = L{vg,ao). Then is right invariant under the lift to TG of the 
right action of G'^^ on G. Suppose that the Legendre transformation FLag is invertible and 
form the corresponding Hamiltonian Ha^ = Ea^ o FL~^^ , where Ea^ is the energy of La^ . Then 
the function H : T*G x V* ^ so defined is S'-invariant and one can apply this theorem. 
At the level of the reduced space, to a reduced Lagrangian / : g x V* — M we associate the 
reduced Hamiltonian /i : g* (s) V* M given by 

/i(/i,a) := - = T7- 

Since 

5h ^ ^ 5h 51 

— =^ and — = - — , 

o/i da da 

we see that the affine Lie-Poisson equations for h on s* are equivalent to the affine Euler- 
Poincare equation (12. 3p for / together with the affine advection equation 

a + a^ + dc(0 = 0. 



6 Hamiltonian Approach to Continuum Theories of 
Perfect Complex Fluids 

This section is the Hamiltonian version of Section [31 Recall that in Section [3] we have applied 
Theorem 12.11 to the case of complex fluids. Here we apply the Hamiltonian analogue of this 
theorem, namely Theorem 15. 3[ 

Recall that for complex fluids we apply the abstract theory to the Lie group G = 
DiS{V) ®J^{V,0) and V* = V^* © V^*. The representations of G on V and V* are of 
the form 

P{v,x)i'^^'^) = i'^V,w{ri,x)), and p(^^^^)-i (a, 7) = (ar^, 7(7/, x)). 
This implies that the infinitesimal actions of g on V and V* are of the form 

{v,w){u,h') = {vu,wu + wu) and (a, 7)(u, z/) = (au, 7U + 7Z/). 

We therefore obtain the diamond operation 

(f , w) o (a, 7) = {v o a + w Oi 7, w 02 7). 
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Since the affine term has the particular form c(?7, x) = (0, C{x)): we obtain the equahties 

dc(u, z/) = (0, dC{iy)) and dJ{v, w) = (0, dC^{w)). 

We now compute some useful expressions in the particular case of complex fluids by using 
the general formulas of section [5l 

By Theorem 15.31 we obtain that the affine Lie-Poisson bracket for complex fluids is 

5f 5g 



m 



K ■ 



om om 

^9 , ^Sf 



V 



+ 
+ 
+ 



V 



7 • 



ad«/ . 

6a 6m 

5f 5g 



+ 



6k 6m 

6a 6m J 

5f6g 6g6f 



6k 6m 



dC 



6'~f 6m 6^ 6k 
Sf\ 6g 



-dC 





6'y 6m 6'y 6k 

6k J ^7 



/i. 



The first four terms give the Lie-Poisson bracket on the dual Lie algebra 

([X(P) ®T{V, o)] (s) [Vi ® V2])* = n\V) X J^(V, 0*) X V* X V*. 

The last term is due to the presence of the affine term C in the representation. Since C 
depends only on the group J-'(V,0), this term does not involve the functional derivatives 
with respect to m. 

The symplectic leaves of this bracket are the affine coadjoint orbits in the dual Lie 
algebra {[X(V)(§)J^{V, 0)] (s) [Vi © V2])*. The expression of the tangent spaces and of the 
affine orbit symplectic forms involves the bilinear form S which is defined in this case on 
[XiV)®T{V,o)](S)[Vi © V2] by 



S((Ui, Ui, Vi, Wi), (U2, Z/2, V2, W2)) = dC{u2) ■ Wi - dC (Ui) ■ W2 

For a Hamiltonian h = h{m,K,a,'j) : fi^(D) x J^(D, 0*) x V{ x - 



the affine 



Lie-Poisson equations of Theorem 15.31 become 



d 


- — £ Sh 
(5m 


m — div 


f 6h' 
\6m^ 


d 
dt 


— ad ah 

Sk 


K — div ^ 


' 6h 
6m 


d 


6h 






a^" = 


om 






d 

Ft^ = 


6h 
-7— 
om 


6h 
-7— - 
6k 


dC ( 



^6_h 

6k 



6h 

6'j 



- . 



6h 
6a 

dC^ 



6h 

^7 



6h\ 
6^) 



(6.1: 



As explained in the previous section, when the reduced Hamiltonian h is defined by a 
Lagrangian / through the Legendre transformation 

61_ 6l_ 
6u 6v , 



/i(m, K, a, 7) = ((m, k), (u, u)) - /(u, z/, a, 7), (m, k) 
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we have 

5h _ 5h _ 5h _ _5l Sh _ _5l 
Sm ^ 6 k, ^ 6a 6a 6'-f 6^ 

Using these equahties, we see directly that the affine Lie-Poisson equations fl6.ip are equiv- 
alent to the affine Euler-Poincare equations (13 ■4p . together with the advection equations 

Using formula ( 15. Sp . the momentum map of the affine right action of the semidirect 
product [Diff (P) (^J^{V, O)] © [Vi © V2] on its cotangent bundle is computed to be 

J„(m^, K^, (f , w), (a, 7)) = {T*r] o + T*x o + w o a + w Oi 7, 

T*L^ o K^ + wo2-f - dC'^{w), (a, 7)) . 

In this formula we need to elaborate on the meaning of the expression T*x o E Q}{T>). 
Thus, by definition, for any Ux G T^V we set ((T*x o -u^^) := (T^x , m^,.) = 

{Ky^{x) ^Txxiux)) ■ In this last expression, recall that : V ^ T*0 covers x'-T^ ^ ^■ 

By the general theory, the nonequivariance cocycle of Jq, is given by a = —0, where is 
computed (using Lemma [5^ to be 

: [Diff(P) (S)J^(P, O)] (S) [Vi © V2] n\V) X T(V, 0*) X V* X V^, 

(f){r], X, V, w) = (dc^(f , w) - {v, w) o c(?7, x), -c(?7, x)) 

= i~w oi C(x), dC^iw) -W02 C(x), 0, -C{x)). 

One can also compute the momentum map J(ao,7o) appearing in the proof of Theorem 15.31 at 
the second stage of reduction. It is associated to the cotangent lifted action of the isotropy 
group 

(Diff(I))(s)^(P,0))(,^^^^) = {(r^,x) I {aovMv,x) + C{x)) = (ao,7o)} 
on the canonical cotangent bundle T* (Diff (P) (§)J^{T>, O)). It is given by 

J(ao,7o) i^v^ f^x) = i'^*V o + T*x o K^, T*L^ o k^) . (6.2) 

Basic example. As in the example given in section [31 we consider the particular case when 
V2 = il^{T>, 0) and the affine representation is 

(a, 7) {at], Ad^-i 77*7 + x'^Tx)- (6.3) 

Recall that in this particular case we have 

7U = £u'y = d'^(7(u)) + iud'^7, 

7z/ = — ad,y7, dC(z/) = di^, and dC^{w) = — div(w). 
The diamond operations are given by 

w Oij = (div'' w) ■ 'J — w ■ i_d'^7 and w 02 7 = — Tr (ad* w) . 
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The affine Lie-Poisson equations (16.11) become 



d_ 

dt' 



-m 



5h 



5h\ 



5m 

5h 



.6h 6h 



£ shUi — div - — m — K ■ d- — o a 



5k 6a 



d_ 
dt' 
d 

d_ 

di 



-K 



5'j J 
— ad^ K, — div 
5h 



1 



5m 



7 



5h 
5m 



^7 
5h 
5m 



— div'*' 



5h 
5^ 



iif^d'''7 — d''' — 
5 k 



(6.4) 



So we recover, by a reduction from a canonical situation, the equations (3.44) of [16,], up 
to sign conventions, as well as their Hamiltonian structure. In matrix notation and with 
respect to local coordinates we have 



rhi 






Kbdi 


[Uo a)i 






dk^a 







a 












. 7" . 












- Cca7 • 






{5h/5mY 
{5h/5KY 

5h/5a 
{5h/5j)l 



(6.5) 



This matrix appears in [20| (formula (2.26a)) and in [16| (formula (3.46)) as the common 
Hamiltonian structure for various hydrodynamical systems. For another derivation of this 
Hamiltonian structure based on reduction see See [5|. 
The associated affine Lie-Poisson bracket is 



{f,9}{m, K,a,-f) 



m ■ 



V 



+ 
+ 
+ 



V 



V 



(5m' 5m 

ad*/ — + d— ■ d— • - — 

«t 5k 5k 5m 5k 5m 



(6.6) 



fx 



5_l5g__5_g_5l_ 
5a 5m 5a 5m 



fx 



V 



"4 

5k 



n \ ^9 



d^^ + £^7 

5k 



51 

^7 



fx. 



The momentum map is computed to be 



Ja(m^, K^, iv,w), 



(«,7)) 



{T*ri o m^ + T*x o k^ + v o a + (div'*' w) ■ •y — w ■ i_ (Pj, 
T*L^ o K-^ + div"' w, (a, 7)) . 



The ^-representation. As on the Lagrangian side, we consider the case where the Hamil- 
tonian is given in terms of the curvature B = d'^7. Recall that the affine representation of 
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Diff (X>) (§)J-'{V, O) on a connection one-form on the trivial principal O-bundle O x V ^ V 
induced by 7 G Q^(V, 0) is 

7(^, X) ■= Ad^-i r/*7 + x'^Tx- 
It induces the linear representation 

B{r],x) ■.= Adl-^7]*B 

on the curvature B = The Lie-Poisson reduction for semidirect products (see Theorem 
II. 2p gives the equations 



d_ 

dt 



m 



— £ shtn — div 



dlV-.l.fi 



6h 
dm 



Jh 6h 
m — K ■ d- — o a 

OK oa 



— • 1 dB 

5B - 



d 
d 

dt 



ad^ K - 

5h 
om 

-£ 6h B 
Sm 



6h 
6m 



div I —K I + Tr I ad^ — 



6h 



(6.7) 



d.dsh B 

6k 



Using formulas fl3.17p . we can prove, as on the Lagrangian side, that the Lie-Poisson equa- 
tions (16.71) are compatible with the affine Lie-Poisson equations (16. 4p . In matrix notation 
and with respect to local coordinates, the Lie-Poisson equations read 



rhi 




rrikdi + dkrrii 




(□oa)i 






" {5h/5mf ' 






dkHa 











{5h/5Kf 


a 




aOdk 













5h/5a 


















_ {5h/5B)f _ 



where 



-Bill + dkBii - diB. 



and 



ijk 



Btj,k + Bkjdi - Bl^dj. 



As before, we recover by a reduction from a canonical cotangent bundle, the Hamiltonian 
structures appearing in [i^] (formula (2.28)) and 16] (p. 152), and we have explained in 
which sense this matrix is Lie-Poisson, as already noted in these papers. More precisely, we 
have found the Lie group which corresponds to the Lie algebra underlying this Hamiltonian 
structure. This group is given by 

[Dm{v)®r{v,o)] (§) [V* © n\v,o)] , 

where Diff (P) (s)JF(X), O) acts on Q'^{V, 0) by the representation 

B ^ Ad^-i r]*B, 

and where the space V{ is only acted upon by the subgroup Diff(X>). 
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The Lie-Poisson bracket is computed to be 



{/,^}(m,K,a,5) 
+ 



+ 



V 



Using formulas (13.161) and 



we obtain that the map 



m 



V 



6m' Sm 



fx 



I ^ ^9 ^9 ^^9 \ 
K ■ adsi — + d— ■ d— ■ - — fi 

Ok ok dm ok dm 



a ■ 



6a 6m 6a 6m 



59 



51 
6B 



£ sf B — ad*/ B I • — [ £ sg B — ad^s B | • — 



5f . , Sf 
5^ = -^^" 



/i. 



(m, z/, a, 7) (m, a, d^7) 

is a Poisson map relative to the affine Lie-Poisson bracket { , } given in fl6.6p and the Lie- 
Poisson bracket associated to the ^-representation. 



7 The Circulation Theorems 

The Kelvin-Noether theorem is a version of the Noether theorem that holds for solutions of 
the Euler-Poincare equations. For example, an application of this theorem to the compress- 
ible adiabatic fluid gives the Kelvin circulation theorem 



d 
di 



vl^ = <h Tds, 

7t >^ It 



where 7t C P is a closed curve which moves with the fluid velocity u, T = de/ds is the 
temperature, and e, s denote respectively the specific internal energy and the specific entropy. 
The Kelvin-Noether theorem associated to Euler-Poincare reduction for semidirect product 



is presented in [21J. We now adapt this result to the case of affine Euler-Poincare reduction. 



Kelvin-Noether Theorem. We work under the hypotheses and the notations of ^ Let 
C be a manifold on which G acts on the left and suppose we have an equivariant map 
}C : C X V* ^ Q**, that is, for all g e G,a e V*, c eC, we have 

m9c,0g{a)),fi) = (/C(c, a). Ad* /i), 

where gc denotes the action of G on C, and 6g is the affine action of G on V*. 
Define the Kelvin-Noether quantity I : C x g x V* ^ R hy 

I 61 
/(c,^,a) := //C(c,a), — (^,a) 



The same proof as in [2l| yields the following result. 
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Theorem 7.1 Fixing cq G C, let ^{t),a{t) satisfy the affine Euler-Poincare equations (12.31) 
and define g{t) to he the solution of g(t) = TRg(^t)C,(t) and, say, g{0) = e. Let c(t) = g{t)cQ 
and I{t) := I{c(t),i(t),a{t)). Then 

|/w^(K(cw,.w),£..-dc-(| 

As we will see in the applications, some examples do not admit a Lagrangian formula- 
tion. Nevertheless, a Kelvin-Noether theorem is still valid for the Hamiltonian formulation. 
Keeping the same notations as before, the Kelvin-Noether quantity is now the mapping 
J : C X 0* X ^ M defined by 

</(c,/i,a) := (/C(c,a),/i), 

and we have the following result. 

Theorem 7.2 Fixing cq G C, let fi{t),a{t) satisfy the affine Lie-Poisson equations of Theo- 
rem \5.3\ and define g{t) to be the solution of 



git) = TRg^,) (^-^^ , <^(0) = e. 



Let c{t) = g{t)cQ and J{t) := J{c{t), fi(t),a{t)). Then 

This result follows from the reconstruction of dynamics as described in ^ Of course, 
when the Hamiltonian h comes from a Lagrangian by Legendre transformation, then this 
theorem is a corollary of Theorem 17. 1[ 

In the case of dynamics on the group G = Diff (P), the standard choice for the equivariant 
map /C is 

(K;(c,a),m) := /-m, (7.1) 

J c P 

where c G C = Emb(S'^, V), the manifold of all embeddings of the circle in D, m G 
and p is advected as {Jrj){porj). There is a generalization of this map in the case of the group 
G = Diff(r')(s)J^(P,0), see §7 in Therefore, Theorems O and O can be applied in 
the case of the affine Euler-Poincare and Lie-Poisson equations (13.41) and (16. ip . Nevertheless 
we shall not use this point of view here and we apply the Kelvin-Noether theorem only to 
the first component of the group G, namely the group Diff(P), and we obtain the following 
result. 

Theorem 7.3 Consider the affine Euler-Poincare equations for complex fluids (13.41) . Sup- 
pose that one of the linear advected variables, say p, is advected as {Jri){p o t]) . Then, using 
the map (17. ip . we have 

^ I ^ - I ^ ( du oa o 
dt J p 6u Jct P\ ^'^ ^7 ^ ^ 
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where Ct is a loop in V which moves with the fluid velocity u. 

Similarly, consider the afflne Lie-Poisson equations for complex fluids (16.11) . Suppose 
that one of the linear advected variables is advected as {Jrj){p o rj). Then, using the map 
(17. ip . we have 

d r 1 r 1 f Sh 5h \ 

— (p -m= (b - [-K - d- ^o«-^Oi7 , 

dt Jet P Jet P \ OK da d'j J 

where Ct is a loop in T> which moves with the fluid velocity u, defined he the equality 

5h 
(5 m 

7-circulation. The 7-circulation is associated to the equation 

d 

77:7 + = -dz^ + ad,, 7. 
at 

Let Tjt be the flow of the vector field u, let Cq be a loop in V and let Ct := rjt o cq. Then, by 
change of variables, we have 



^ 7 = ^ ^t7 = (7 + = Vti-^t^ + ad^7) = jf ad^7 G 0. 



8 Applications 
8.1 Spin Systems 

We thank D. Holm for challenging us with this example as a simple model for many applica- 
tions. It illustrates the applicability of Theorems 12. II and 15.31 in a very simple situation that 
exhibits, nevertheless, some of the key difficulties of more complicated fluid models. Let V 
be a manifold and O a Lie group thought of as the "order parameters" of some fluid model. 
Take G = J-'{T>, O) 3 x ^* = ^^{T>, 0) 3 7. Consider the affine G-representation on V* 
given by 

^^(7) :=Ad^-i7 + X"'Tx. 

This is simply the gauge transformation of the connection on the trivial principal (9-bundle 
O X V ^ V induced by 7 (see (13.81) ). The associated diamond operation is given by 

w o 7 = — Tr(ad* w). 

Since c{x) = x~^Tx we obtain, as in the case of complex fluids, 

dc{iy) = du and dc^(w) = — divw. 

As before we use the notations u G JF(D, 0) and w G X(T>,o*). The affine Euler-Poincare 
equations (12. 3p become 

— — = -ad^ — + div*^ — , (8.1) 
at ou ov 07 
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where (ad* k) (x) = ad*(^)(K(x)). The evolution equation for 7 is 

d 

7 + d^u = 0. 

ot 

Similarly, the affine Lie-Poisson equations of Theorem 15.31 become 

—K = — adsh K — div ' — 

dt s. 57 

at OK 

In the particular case D = and O = 50(3), this system of equations appears in the 
context of the macroscopic description of spin glasses, see equations (28) and (29) in and 



references therein. See also equations (3.9), (3.10) in [22|, system (1) in [23| and references 
therein for an application to magnetic media. In this context, the variable k is interpreted as 
the spin density, v is the infinitesimal spin rotation and the curvature d'^7 is the disclination 
density. 

Interesting choices for the Lagrangian are 

1 /• ,2 . , . 1 



/±(^,7) = / Ilh7]||>, hi^.l) = / IIM^,7)II>, 

and 



9 / H L 3 /J II r 3 II \ ' ' / 9 



Isg{^,i)='t; / IIHIV-^ / hllV (8.3) 



2 JD 2 j-p 

In the expression of the Lagrangian l^, [z^, 7] denotes the o-valued one-form on T> given by 
Vx ^ [i^{x)j'l{x){'Vx)] £ 0. The norm is associated to the metric (gk), on the vector bundle 
of o-valued one-forms on V, induced by a Riemannian metric g onV and by an Ad-invariant 
inner product k on 0. More precisely, the metric (gk) on the vector bundle of o-valued 
fc-forms on V is given in the following way. The Riemannian metric g induces a Riemannian 
metric 'g on the vector bundles A^V — X> of exterior /c-forms on V. For a^, Px ^ A^{V, o)x, 
we can write ax = a°'fa and Px = P°'fa, where {/„} is a basis of and a'*,/?" G (A'^M)^;. So 
we define 

{gk)x{ax,f3x) ■■= kabg{a'',l3^), 

where kab '■= k{fa, ft)- This construction is independent of the basis. If the manifold is taken 
to be M and if C = 50(3), then the Lagrangian l±_ reads 

This choice is reminiscent of the Skyrme model, a nonlinear topological model of pions in 



nuclear physics, see [31 



In the expression of the Lagrangian /y, fc(z^, 7) denotes the one-form on T> given by 
Vx I— > k{i'{x),'y{x){vx)). The norm is taken relative to the Riemannian metric g on V. 

Interestingly, when the Lagrangians l± or /y are used, the affine Euler-Poincare equations 
0|) simplify to 

d 61 51 

-— = div— , 8.4 

otdu d'j 
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where on the right hand side we have the usual divergence operator in M.^. 

In the expression of the Lagrangian Isg, "we used the symbols e and p for the constants of 
susceptibility and the rigidity. The norms are respectively associated to the inner product k 
and to the metric {gk). The associated Hamiltonian reads 



h 



SG 



and is used, with O = 50(3), in the context of the macroscopic description of spin glasses, 
see expression (26) in j^. In this case the affine Lie-Poisson equation (18.21) reads 




where 7* G X{V, 0*) and G J^{V, 0) are associated to 7 and k via the metrics. 
More general expressions, such as 



Jv ^^1 Jv ^ Jv ^ Jv ^ Jv 



are used in the theory of magnetic media, see e.g. equation (3) in 23 
Lagrangian reduction. Consider a Lagrangian 

L : TT{V,0) X 9}{V,o) ^ R, L = ^(2/^,7) 

such that, for all %Ij G J-^iV, O) we have 

L{TR^ o i/^, Ad^-i 7 + ij'^Tijj) = L{u^, 7). 

Let X be a curve in J^(V,0) and consider the curve u := TR^~i o x. For 70 G Q}{V,o) 
consider the solution 7 of the equation 



-7 + d-. 







with initial condition 70. This solution is given by 7 = Ad;^7o + X^X ^- Then, by Theorem 
12.11 X is a solution of the Euler-Lagrange equations associated to L^^ 



if and only if v is 

solution of (18. ip . Note that in the special case 70 = 0, the evolution of 7 is given by the 
important relation 

7 = XTX~^ 

and the disclination density vanisheSjthat is, d7 + [7,7] = 0. This hypothesis is usually 



assumed in the examples treated in [8| and [23[, and is referred to as the Maurer-Cartan 
constraint. Recall that the vanishing of the curvature is preserved by the flow of (18. 2p . In 
our approach the variable x can be interpreted as the Lagrangian evolution of the spin. 

Hamiltonian reduction. Consider a Hamiltonian 

H : T*J^{V, O) X 9}{V, 0) ^ R, H = H{k^, 7) 
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such that, for all ip G T{T>^ O) we have 

H{T*R^-i o /t^, Ad^-i 7 + i:-'Ti:) = H{k^, 7). 

Then, by Theorem 15. 3[ a curve G T*J^{V, O) is a solution of the Hamilton equations 
associated to if^^ if and only if the curve 

K := T*R^ o 

is a solution of the affine Lie-Poisson equation (18. 2p . The associated affine Lie-Poisson bracket 
is given by 

Thus, we have recovered the Poisson bracket of j^, by reduction of the canonical structure. 
The momentum map is 

Jc,(kx, (w, 7)) = iT*L^ oK,x + div"^ 7) 

and the reduced symplectic spaces are affine coadjoint orbits in the dual Lie algebra 
{J-'(T>, 0) (§)fl^{V, 0))*. Using formula (15.71) . we obtain that they are given by 

01^^) = {(Ad*«: + div^>^(^)^,^^(7)) I ix,w)eJ'iV,O)(S)XiV,0*)}. 
8.2 Yang-Mills Magnetohydrodynamics 

Magnetohydrodynamics models the motion of an electrically charged and perfectly conduct- 
ing fluid. In the balance of momentum law, one must add the Lorentz force of the magnetic 
field created by the fiuid in motion. In addition, the hypothesis of infinite conductivity 
leads one to the conclusion that magnetic lines are frozen in the fluid, i.e. that they are 
transported along the particle paths. This hypothesis leads to the equation 

d 

-B + £^B = 0. 

ot 

This model can be extend to incorporate nonabelian Yan g-M ills interactions and is known 



under the name of Yang-Mills magnetohydrodynamics; see [18| for a derivation of this model. 
Recall that for Yang-Mills theory, the field B is seen as the curvature of a connection A on 
a principal bundle. Clearly the connection A represents the variable 7 in the general theory 
developed previously, on which the automorphism group acts by affine transformations. This 
shows that the abstract formalism developed previously is very natural in the context of 
Yang-Mills theory. Note that there is a more general model of fiuid motion with Yang-Mills 
charged particles, namely the Euler- Yang-Mills equations. The Hamiltonian structure of 



these equations is given in [13j, see also [l2j for the associated Lagrangian and Hamiltonian 
reductions. 
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As remarked in [20(, at the reduced level, the Hamiltonian structure of Yang-Mills mag- 
netohydrodynamics is given by the matrix (16. 5p . In this paragraph we carry out the corre- 
sponding affine Lie-Poisson reduction. 

The group G is chosen to be the semidirect product of the diffeomorphism group with 
the group of C- valued function on P, that is, G is the group Diff (P) (s)^(P, O). The order 
parameter Lie group O represents here the symmetry group of the particles interaction. For 
example, O = corresponds to electromagnetism, O = SU{2) and O = SU (3) correspond 
to weak and strong interactions, respectively. The advected quantities are the mass density 
p, the specific entropy s, and the potential of the Yang-Mills field A. Therefore, we set 



The action of {rj, x) € G on (p, s) G V{ is the usual right representation of the fluid relabeling 
group on the mass density and entropy. It is given by 

The right affine action of (77, x) ^ G on A & Q^{M, 0) is given, as in the example (16.31) . by 

A ^ Ad^-i 7]*A + x'^Tx. 

Since the variable n G J-'(T>, 0*) is interpreted as the gauge-charge density, we use the nota- 
tions Q G TiV, 0*) = T*J^{V, O) and G T;^J^{V, O). 

The Hamiltonian i/(p,s,A) : T*(Diff (P) ®T{V, C)) ^ M, is given by 



where e denotes the specific internal energy, the norm in the first term is associated to a 
Riemannian metric g on "D, and the norm in the third term is associated to the metric 
{gk), on the vector bundle of o-valued fc-forms on P, induced by the metric g and by an 
Ad-invariant inner product k on 0. For details see, §8.11 

The metric {gk) can be used to identify fifc(P, 0*) and its dual VL^{T>, 0), by raising and 
lowering indices. Through this identification, the operators div and div"^ act also on Vt^{T>, 0). 
The Hamiltonian if(m^, Q^, p, s, A) is invariant under the right action of {ip, ip) given by 



(m^,Qx,p,s, A) 

(-^ (J(y9(m^ o y?), J(f(T*R^-i o Q^o Lf), Jip{p O Lp),s O Lp, Ad^-1 Lp* A + tp~^T^) . 

Therefore, the hypotheses of Theorem 15.31 are satisfied and the reduced Hamiltonian h : 
n^{V) X T{V, 0*) X V{ xV^ is given by 



a= {p,s) eV{ = J^{V) X J^(V) and A = -fen\V,o). 




(8.6) 



h{ni,Q,p,s 
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where the norms are respectively associated to the metrics g and (gk). We now compute the 
affine Lie-Poisson equations (16.41) associated to this Hamihonian. The functional derivatives 
are 

6h 1 u ^ 
u:= — = -m», z/:= — = 
om p oQ 

and 

5h 1 ,, ,,9 de 5h de 6h , a -.a . 

^ = --||uf + e + p^. - = p^. _ = -d,v-d-A 

where div"^ is defined by (13.141) . The advection equation are 

+ div(pu) = 0, + ds(u) = 0, + d^(A(u)) + i^B = 0. 

The equation for the gauge charge is 

^Q = - div(Qu) + div^ div^ d^A = - div(gu). 

Indeed, for all / G J-'(T>, o) we have 

[ fc(div^div^d^A,/)p= / {gk) {d^A,dM^f)p= [ (gk) (5, [E, /]) /x= 0, 
Jv Jv Jv 

where we used the equality we have d"^d"^/ = [B, /] for B = d'^A and the fact that, in an 
orthonormal frame with respect to g, we have 

(gk) {B, [B, /]) = Bt^[B, f%kat = Bt^[Bi„ f]'ka, 

= k{B,„ [B,^, /]) = -k{[B,^, 5,,], /) = 0. 

This proves that div^ div"^ d'^A = 0. 

Using the advection equation for p we obtain the equality 



+ + (div u)m j = P i + VuU + Vu^ ■ u ) . 



We also have 



f6h 6h A" P J,, ,,2 J J de\ de ^ 

\Jp ^ ~^ Js ^ J - pgrade - pgrad I p— I + p— grads 

p-Q^J - p^gradp 
dp, 



p( Vu ■ u) — grad I p — 



and 



— Oi A = div^ TT M - TT ■ 1 
5 A \ 5AJ 5A - 



+ (^A;)(div^ d^AXd^'A) 
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Therefore we obtain that the first hne of (16.41) becomes 

— u + VuU= -- (^gmdp+{gk) (i.5,div^5)' 
where B = d"^A and p = p^^. In summary, we have obtained the equations of Yang-Mills 



magnetohydrodynamics 
^ d 



^^u + VuU = -- (gradp + {gk) div^ 5)*) , B = d^A, 
-Q + div(gu) =0, -A + d^(A(u)) + i^B = 0, 



•7) 



d 

—p + div(pu) = 0, 



9 1 / X 

-. + d.(u) 



0. 



We now treat the particular case of magnetohydrodynamics, that is, the case O = S^. 
In order to recover the standard equations we suppose that V is three dimensional. In 
this case we can define the magnetic potential A := A" G X(T') and the magnetic field 
B := (^-B)" € X{V). Since the group is Abelian, covariant differentiation coincides with 
usual differentiation and the equality d^A = B reads curl A = B. Using the identities 
(divi?)'* = — curlB and (iu-B)" = B x u we obtain 

g{i_B,divBf = - (i(divB)«5)* = B X curlB. 

Suppose that all particles have mass m. The electric charge q is such that Q = p—, therefore 
the equation for Q in (18.71) becomes 



d_ 
di 



q + dg(u) = 0. 



If we suppose that at time t = all the particle have the same charge, then this charge 
remains constant for all time. By making use af these remarks and hypotheses, equations 
(18. 7p become 

^ d 1 

— u + VuU = — (gradp + B x curl B) , B = curl A, 
ot p 

d 

A + grad[^(A,u)] + B X u = 0, (8.8) 



dt 
d_ 
di 



p + div(pu) = 0, 



d_ 

dt 



s + ds(u) 



0. 



Thus, we have recovered the equations for magnetohydrodynamics. Turning back to the 
general case and using Theorem 15.31 we obtain the following result. 

Hamiltonian reduction for Yang-Mills magnetohydrodynamics. A smooth path 
(m^,Q^) G T*[DiS(V) (§)J^{V,0)] is a solution of Hamilton's equations associated to the 
Hamiltonian i^(po,so,Ao) given in (18. 6p if and only if the curve 

(pu^ Q) =: (m, Q) := Jir]'') (m, o r]-\ T*R^^^-. (Q^ o r^-^)) 
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is a solution of the system (18.71) with initial conditions (po, Sq, Aq). 
The evolution of the advected quantities is given by 



P = JiV ^)ipO°V S = SqOT] \ 

A = Ad^o^-i r],Ao + ix o V'')Tix'' ° V'') = V* (Ad^ + xTx'') ■ 

This theorem is interesting from two points of view. Firstly, it allows us to recover the 
non-canonical Hamiltonian structure given in 20| by a reduction from a canonical cotan- 
gent bundle. Secondly, it generalizes to the nonabelian case the Hamiltonian reduction for 
magnetohydrodynamics given in 28 . 

The associated affine Lie-Poisson bracket is that given in (16. 6p . where the third term 
takes the explicit form 



bp) 5m \6p) 5m/ \ \5s5m/ 



Since the Hamiltonian depends on A only through its curvature, the equations can be for- 
mulated using the 5-representation. One simply replaces the equation for A by those for its 
curvature B, namely, 

at 

In particular, we obtain that the force term ((jf/c) div"^ depends only on the curvature 
B and not on the connection one-form A. 

In the particular case of magnetohydrodynamics, the evolution equation of B reads 

d 

— B + curl(B X u) = 0. 

ot 

In the general case of Yang-Mills magnetohydrodynamics, the Kelvin- Noether theorem 
gives 

d 

and the 7-circulation gives 

dt 

where q is a loop which move with the fluid velocity u, that is, Ct = rjt o cq. 



({ W' = ({ Tds - I -{gk) (\_B,div^B) 

J Ct J Cf J Ct P 



i<i> A = 0, 



8.3 Hall Magnetohydrodynamics 

As we will see. Hall magnetohydrodynamics does not require the use of the affine Lie-Poisson 
reduction developed in this paper. However, in view of the next paragraph about superflu- 
ids, we quickly recall here from ij] the Hamiltonian formulation of these equations. We 
will obtain this Hamiltonian structure by a Lie-Poisson reduction for semidirect products, 
associated to the direct product group G := Diff(r') x DiSiV). The advected quantities are 

(p, s; n) G J^{V) x J^{V) x J^{V). 
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The variables p and s are, as before, the mass density and the specific entropy, on which 
only the first diffeomorphism group acts as 

(p, s) ^ {{Jr]){por]),sor]). 

The variable n is the electron charge density, on which only the second diffeomorphism group 
acts as 

By Lie-Poisson reduction, for a Hamiltonian h — h{m, p, s; n, n) defined on the dual Lie- 
algebra 

{[X{V)(§){J^{V) X X [X(P)(DJr(p)])* 

^ [Q\V) X J^{V) X J^{V)] X [n\V) X r{V)], 



we obtain the coupled Lie-Poisson equations 



( d 

— m — —£ 6hni — div 

at ^ 



6h 



6h 



6h 



r I ^ ^^^•^ 

dm J op OS 



d f 5h 



dt^ \Sm 
d / Sh 



P 



(8.9) 



and 





— n = —£6hn — div 
9i ^ 



-— n — — on 
on I On 



d 



—n = — div —n 



dt 



Sh 



Sn 



(8.10) 



Prom the second Lie-Poisson system we obtain that the evolution of n/n is given by 

(8.11) 

The Hamiltonian for Hall magnetohydrodynamics is 



^(n/n) = -£6h(n/n) -grad^. 
ot on 



/i(m, p, s; n, n) := 
where the one-from A, defined by 



/- 


m -A 


Jv P 


R 



'p + [ pe{p,s)p+l [ WdAfp, (8.12) 
Jv ^ Jv 



n 



n 



is the magnetic vector potential, the constants a, R are respectively the ion charge-to-mass 
ratio and the Hall scaling parameter, and the norms are taken with respect to a fixed Rie- 
mannian metric g on T>. The functional derivatives are computed to be 



u 1^ = im» - ^A\ 
dm p R 



Sh 1.. ..n a . de Sh de 
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6h 



V := 



ap 



R 



6h ap 



R 



— = -^u - -{divBf, p = . u + ■ (div5)« = --A ■ v, 



Sn n n ^ 5n ' ri^ ' R 

where B := dA. Recall from (13 .14^ that div is defined on VLk{T>) as the negative of the adjoint 
of the exterior differential d on Vt^{T>). However, the metric g on T> gives an identification 
of nfc(P) with n''{V) and hence we can regard div as defined also on Q^{T>). It follows that 
on Q^{V) we have div = —6, where 6 is the usual codifferential induced by g and d. 

The variable v is interpreted as the electron fluid velocity. The advection equations for 
p, s, and n are given by 



d_ 

di 



d 

p + div(pu) = 0, —s + ds(u 



d 

0, and —n + div(nv) 



0. 



Using the expression of v in terms of u we obtain that div(r2v) = — adiv(pu) (since B = dA 
implies that div ((divi?)^) = 0) which proves that 



0. 



Thus, if we assume that apo + hq = for the initial conditions, we have ap + n = for all 
time. 

Using the definition of A and (18. Ill) , we obtain that the equation for A is given by 



dt 



Using the equations for A and p, we obtain the following equation for u: 



d_ 

di 



U + VuU 



p 



gradp 



ap 

'r 



, P = P 



, de 
dp 



Suppose that the initial conditions po and uq verify the equality apo + uq = 0. As we have 
seen above, the equality remains valid for all times, and we obtain v = u + ^(div B)K Using 
this, the equations above simplify and we obtain the system 



d_ 

dt 
d_ 

di 
d_ 



U + VuU 



gradp- (i(divi3)»5)' 



P 



A 



div (pu) 



d_ 

di' 



—ds (u) 



(8.13) 



-LB 



R_^ 
ap 



(divB)tt 



B 



When V is three dimensional, we can define the magnetic potential A := A^ and the 
magnetic field B := (^-B)". In this case the previous equations read 



d 1 
u + VuU = — (gradp + B x curl B) 



dt 
d_ 

di 
d_ 

di' 



P 



p = — div (pu) 



R. 



d_ 

di' 



-ds (u) 



(8.14) 



u X B H B X curlB. 

ap 
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These are the classical equations of Hall magnetohydrodynamics. Note that we can pass 
from the equations for magnetohydrodynamics to those for Hall magnetohydrodynamics by 
simply replacing the advection law for A by the Ohm 's law. In terms of the magnetic field 
B, one simply replace the advection law 

where u is the fluid velocity, by the equation 

d 



dt 



B + £^B = 0, 



where v is the electron fluid velocity. 

In fact, the Hamiltonian h given in (18.121) . is the value at the identity of the right invariant 
Hamiltonian if(m^, p, s; n^, n) = if(p^s:n)(m^; n^), where 

H^p^s;n) : T*(Diff (P) X Diff (P)) ^ R 



is given by 



V 



ap 



o ^ ^ or] 



V 



no ^ ^ o r] 
n o 



2 




/U + 


i 



pe{pJr] \s)/i 



Hamiltonian reduction for Hall magnetohydrodynamics. Suppose that apo + no = 0. 
A curve (m^, n^) G T*(Diff(I') x Diff(I')) is a solution of Hamilton's equations associated 
to -f^(po,so;no) if and only if the curve 

(m, n) := ( J(r ')(m, o r]-'), J(r')K ° T')) 

is a solution of the equations (18.131) where A := ^n = — ^n, since ap + n = 0. Moreover 
the evolution of the advected quantities is given by 

p = J(7^-i)(po 77~^), s = soor7-\ n = J(^"^)(no o^-i). 

Let us assume from now on that the initial conditions po and no are related by apo + no = 0. 
We have seen that this implies that ap + n = 0. From the relations above we conclude the 
interesting result that the action of o r] fixes no, that is, J(^~^ o 77) (no o ° v) = ''^o- 
Conversely, given this relation and the condition apo+no = 0, it is easily seen that ap+n = 0. 

The Lie-Poisson bracket associated to these equations is clearly the sum of two Lie-Poisson 
brackets associated to the semidirect products Diff ("D) (s) [J^{V)xJ^{V)] and Diff ("D) (s) J'(I?). 

The Kelvin-Noether theorem associated to the variable m gives 
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which can be rewritten as 



where q is a loop which moves with the fluid velocity u, that is, Ct = rjt o cq. The Kelvin- 
Noether theorem associated to the variable n gives 

i<i A = 0, 



dt 

where dt is a loop which moves with the electron fluid velocity v, that is, dt = C,t° d^. 



8.4 Multivelocity Superfluids 

Superfluidity is a rare state of matter encountered in few fluids at extremely low tempera- 
tures. Such materials exhibit strange behavior such as the lack of viscosity, the ability to flow 
through very small channels that are impermeable to ordinary fluids, and the fact that it 
can form a layer whose thickness is that of one atom on the walls of the container in which it 
is placed. In addition, the rotational speed of a superfluid is quantized, that is, the fluid can 
rotate only at certain values of the speed. Superfluidity is considered to be a manifestation 
of quantum mechanical effects at macroscopic level. Typical examples of superfluids are ^He, 
whose atoms are fermions and the superfluid transition occurs by Cooper pairing between 
atoms rather than electrons, and '^He, whose atoms are bosons and the superfluidity is a 
consequence of Bose-Einstein condensation in an interacting system. 

For example at temperatures close to absolute zero a solution of ^He and ^He has its 
hydrodynamics described by three velocities: two superfluid velocities v], and one normal 
fluid velocity v„. If other kinds of superfluid are present, one needs to introduce additional 
superfluid velocities. For a history of the equations considered below and the Hamiltonian 
structure for multivelocity superfluids see [l9l]. 

The two-fluid model. We first treat the two-fluid model, that is, the case of one super- 
fluid velocity and one normal-fluid velocity v„. Remarkably, this Hamiltonian structure 
can be obtained by affine Lie-Poisson reduction, with order parameter Lie group O = S^. 
In this paragraph we also carry out the corresponding Lagrangian formulation, by apply- 
ing the general theory of affine Euler-Poincare reduction to the semidirect product group 

The linear advected quantity is the entropy density S on which a diffeomorphism t] acts 

as 

S 1-^ {Jri){S o rj). 

The affine advected quantity is the superfluid velocity v^, on which the element {r],x) ^ 
I)iS{V)®J^{V,S^) acts as 

This action is simply the affine representation (13. 6p for the Lie group O = S^. Here the 
advected quantity is a vector field and not a one-form, since it represents a velocity and 
hence the formula (13. 6p was changed accordingly. As will be seen, in this formalism, the 
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mass density does not appear as an advected quantity in the representation space V*; it is 
a momentum, that is, one of the variables in the dual Lie algebra q* = x JF(X)). 

The reduced Lagrangian for superfluids is 

/ : [X{V)(§)J^(V)] (S) [J^(V) © X(V)] — > R 

given by 

/(v„,z/, S',v^) := - / p||v„f/i+ / {p'^)fJ'- / e{p,S,\s -^n)fJ', (8.15) 
^ Jv Jv Jv 

where v„ is the velocity of the normal flow. The internal energy density e is seen here as 

a function of three variables e = e{p, S,y) : M x M x TV — > M. The norm in the first term 

is taken relative to a fixed Riemannian metric g on V. Note that for superfluids it is more 

convenient to work with the internal energy and entropy per unit volume and not per unit 

mass as in the preceding examples. We make the following definitions: 

/ichcm := -g^ip,S,Vs - v„) G J^iV), T := ^(P>'S',v^ - v„) G J^iV), 

de 

p:=— (p,5,v, -v„) en\V). 
or 

The interpretation of the quantities pchem^T, and p is obtained from the following thermo- 
dynamic derivative identity for the internal energy (superfluid first law): 

d{e{p, S,\s - v„)) = /ichemdp + TdS + p ■ V.(vs - v„) G ^^(V), 

where V denotes the Levi-Civita covariant derivative associated to the metric g. The function 
Pchem is the chemical potential, T is the temperature, and p is the relative momentum density. 
The mass density p is the function p = piy^, u, S, v^) defined implicitly by the condition 

de 1 

/Wchcm := i:-{p,S,Ws - V„) = -||v„f + (8.16) 

op 2 

Therefore, p is not a variable in this approach. By the implicit function theorem, the relation 
above defines a unique function p, provided the function e verifies the condition 

— (r, s, v^^) ^ 0, for all (r, s, i;^) G M x R x TV. (8.17) 

The afRne Euler-Poincare equations for the two-fluid model. Using the definition 
f l8.16p of the function p, we compute below the functional derivatives of /. We have 



51 ^ 1 2 ^P 



P^n + o l|Vn|| ^ \- J^T{ /ichcm^ h P = pV„ + p =: III 



51 1,, ..odp dp dp 

Y- = o Vn ^ + P + Pchem^ = P, 

5v 1 ov ov ov 

51 1 II 2 "^p dp dp 

Is ^ 2 ds^^'Is' ^'"^'^ds - ^ = 

51 1 .. ..n dp dp dp 
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Using the affine Euler-Poincare equations (13.131) . we obtain the following equations for p, S, 
and v^: 





—p + div(pv,„ + p«) = 0, —S + div(5v„) = 0, 

d_ 

di 



v^ + grad ^5'(vs,v„) + Pchem - ^||vnf^ + (ivndv^)" = 0. 



The last equation can be rewritten as 



d_ 

dt 



grad I^Pchcm - ^llv^ - v„f ^ + (iv.-v^dv^)" (8.1^ 



When T> is three dimensional, the last term reads 

(iv,-v„dv^) = (*dv^) X (v^ - v„) = curl x (v^ - v„) = (v„ - v^) x curl v^. 
Using the equality 

-pd ^/Xchem - ^||vn||^^ " SdT = -dp - p ■ V.v^ + m ■ V.v„, 

where p := —e{p, S, — v„) + pchemP + ST is the Euler pressure law, the equation for m is 
computed as follows: 

= --Cv„m - div(v„)m - pd ^pchcm - ^||v„f ^ - SdT 

-div(p«)v^ + dvt(_,p«) 
= -£^^ m - div(v„)m + m ■ V.v„ - dp - p ■ V.v^ 

-div(p«)v^ + dvt(_,p«) 
= -Vv„m - div(v„)m - dp - VptiV^ - div(p'')v^ 
= - Div T, 

where T is the (1, 1) superfluid stress tensor defined by 

T := v„(g)m + p«®v^+p5, 

5 is the Kronecker (1, 1) tensor, and Div is the divergence of a (1, 1) tensor, defined as the 
trace of the bilinear map 

{a,v)^W,T{a,_)&n\V). 

In coordinates, we have 

T} = v^nij- + pV,j- +p5i 

and 

(DivT), = (VT)}, = d{T] + Tjri, - T\T\^. 
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The equations for superfluid dynamics are therefore given by 
f d 



where 



m 



DivT, 



dt 
d 

—p + div(pv„ + p" 



0, + div(^v, 

ot 



— v^ + gradf 5((vs,v„) +pchem- 2l|vnf ) + (iv„dv^) = 0, 



^.19) 



T := v„ (g) m + p" (g) + and p := -£:(p, S", - v„) + /ichcmP + ST. 

Thus we have recovered the equations (la) - (Id) in p^], in the particular case of the two- 
fluids modeL By Legendre transformation of the reduced Lagrangian fIS.lSp . we obtain the 
reduced Hamiltonian 



/i(m,p, S", Vs) = / p||v„||V+ / (m-v„)/i+ / e(p, 5, - v„)p, 
^ Jv Jv Jv 

where v„ = v„(m, p, 5, v^) is the vector field defined by the implicit condition 

m - Pv^ = |^(P> \s - Vn) =: P- 
or 



^.20) 



(8.21) 



By the implicit function theorem, the above relation defines a unique vector field v„, provided 
the function e verifies the condition that 



(r, s, v^) ■ Ux - ru:, 



is a bijective linear map. If f l8.17p holds, then this condition is equivalent to saying that 
the Legendre transformation is invertible. Of course, the functions I and h are the values 
at the identity of the corresponding unreduced right invariant Lagrangian and Hamilto- 
nian L(v^,z/x,S',Vs) = L(5,v,)(v^,z/;^) and if(m^, p^, S*, v^) = H^s,ws)ij^v^ Px)^ where /^(5,v,) : 
T(Diff (P) (S) J^(I?, S^)) M and iy(5,v,) : T*(Diff (P) (S) J^(P, S^)) R. 
Lagrangian reduction for superfiuids. A curve (77, x) ^ Diff (X>) (^J^{T>^ S^) is a solution 
of the Euler-Lagrange equations associated to the Lagrangian Li^So,vso) if ^'^d only if the curve 

(v„, u) := (?) o r]-^ (TR^-ix) o r]'^) 



is a solution of the superfiuids equations fl8.19p with initial conditions (Scv^o). 
The evolution of the advected quantities is given by 

S = J{v-'){So0 7^-') and v,= (r^,(vto + dx-^))^ 
Note that the evolution of the superfluid vorticity is given by 
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therefore, the irrotationahty condition dv^ = (curlvg = for the three dimensional case) 
is preserved. 

Hamiltonian reduction for superfiuids. A curve (m^,p^) in the cotangent bundle 
T*[DiS{V) (§)J^(V, S^)] is a solution of Hamilton's equations associated to the superfluid 
Hamiltonian H(^Sq,vso) ^^"^ only if the curve 

(m,p) := J{ri'^) {m.^ o rj-^ o r]'^) 

is a solution of the system (18.191) with initial conditions (5*0, v^o). 
The associated Poisson bracket for superfiuids is 



{/,fi'}(m,p, S, \s 



m 



V 



5m' 6m 



+ / o • I d— • - d— ■ - 

ob dm ob dm 



P 



+ 



V 



'r 

5p 



dVg \ op ' 



p- 



(8.22) 



Multivelocity superfiuids. We now quickly explain how to generalize the preceding ap- 
proach to the case of superfiuids with m velocities. Consider the semidirect product 

G := Difr(P) (S) [J^(V, S^)x ...X J^{V, 



where the group on the right is the direct product of the groups J^(T>, S^). The semidirect 
product is associated to the right action of Diff ("D) given by 

The affine advected quantities are the m superfluid velocities (vj, v™) G X(V)"^ on which 
an element (77, x^, x™) acts as 



(v],...,vr)^ 



The reduced Lagrangian is defined on [X{V) (D J^(r')™] (s) [J^{V) © X(P)"'] and is given by 



/(v„,(j."),5, (v:)): = 



I I Phn\\p + Y. I P"^"- / ^((P"),5,(v:-v„))p, 



where p", the mass density of the condensate particles with flow v^, is the function p" = 
p^iyni i^""), S, (v")) defined by the condition 

/^chem := ^((P")' K - Vn)) = 2 llVnf + 

By the implicit function theorem, these conditions uniquely determine p", provided the 
matrix 



dp'^dpi^ 



al3 
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is invertible for all r*, s, f ^ G M x M x TV. The variable p is defined by p :— + ... + p^ 
and denotes the total mass density. Using the notations 

we obtain the functional derivatives 

^ = ,v„ + ^p, — = p, - = -T. — = -p 

a=l 

and the thermodynamic derivative identity for the internal energy 

m m 
d(£((p"), S, (V^ - V„))) = 5] //S,e.dp° + Td5 + 5] p'^ • V. (V,« - V„) e 

a=l a=l 

Using the affine Euler-Poincare equations associated to Diff (D) (§) [J^{T>, S^) x ... x J^{T>, S^)], 
we obtain the equations for multivelocity superfluids 

-m = -DwT, 

-p" + div (p"v„ + (p")«) =0, -S + diviSvn) = 0, (8.23) 
^v? + grad (^g{v^, vj + //".em - ^l|v„||') + (iv„d(v,")^)" - 0, 

where a = 1, m. The stress tensor T and the pressure p are given by 

m 

T : = v„ ® m + J2iP"f ® «f +P^^ 

m 

p : = -e{{p"), S, {v^ - v„)) + /^chemP" + ST. 
By Legendre transformation, we obtain the Hamiltonian 

/i(m,(p«),5,K)) = -^ / p||vnHV+ / (m-v„)/x+ / £((p"),5,(v«-v„))/., 

7l) 7x1 Jv 

where v„ = v„(m, (p"), 5", (v")) is the vector field defined by the implicit condition 

a=l 
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By the implicit function theorem, the above relation defines a unique function v„, provided 
the function e verifies the condition that the linear map 



E 

a,/3=l 



{r,s,Vcc) ■ Ucc - rux 



is bijective. 

Lagrangian and Hamiltonian reductions hold as in the two-fluids model. The evolutions 
of S and are given by 

S = J{v-'){So o r') and v," = r/, {{^r%Y + d(x")-^)" , 

and the irrotationality condition d(v°)^ = is preserved. 
The associated Poisson bracket is given by 



{/,5'}(m,p, S*, 



m 



(5m' 6m 



a=l 

-I' 

JV 

m 

a=l 



+ V / ° ■ f d— ■ — - d— • — ] 



5S 5m 5S 5m 



V 



5v? 



5r 



/X. 



(8.24) 



The 7-circulation gives 



d 
dt 



0, for all a = 1, ...,m, 



where Cf is a loop which moves with the normal fluid velocity v„. 



8.5 Superfluid Yang-Mills Magnetohydrodynamics 

In this paragraph we combine the Hamiltonian structures of Yang-Mills magnetohydrody- 
namics and superfluid dynamics, to obtain a new physical model for the theory of superfiuids 
Yang-Mills magnetohydrodynamics as well as the corresponding Hamiltonian structure. In 



the Abelian case we recover the theory and the Hamiltonian structure derived in [19|. We 
need a slight generalization of the geometric framework developed in §3] and §H1 namely we 
consider the group semidirect product 

DifT(P) (s) (^(P, O) X J^(V, S^)), 

where J-^{T>, O) x 5^) is a direct product of groups on which Diff (P) acts as 

(Xi,X2) ^ (Xi o^?,X2 0^?)• 
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The affine advected quantities are the potential of the Yang-Mills fluid A and the superfluid 
velocity v^, on which (?7,Xi;X2) acts as 

Ad^-iri*A + xi^Txi and v, ^^ (r^V^ + dx2)*. 

The reduced Hamiltonian is defined on the dual of the Lie algebra 

[X{V) (s) {J^{V, o) X J^{V))] (s) (J^(P) X n\V, o) X X{V)) 

and is given by 

h{m,Q,p, S,A,Vs) = 



o / PW^nW P 



V 



m ■ v„)/i 



where v„ is the normal fluid velocity defined as in fl8.2ip . This is simply the Hamiltonian 
(I8.20p plus the energy of the Yang-Mills field. The norms are respectively associated to the 
metrics g and {gk), where g is a. Riemannian metric on V and k is an Ad-invariant inner 
product on o. The affine Lie-Poisson equations associated to this Hamiltonian are computed 
to be 

—p + div(pv„ + p« 



— m = - Div T, —p + div(pv„ + p») = 

-Q + div(gv„) = 0, -S + div(5v„) = 0, 
ot at 

d_ 

di 



v^ + grad (^gi-Vs,^n) + Pchem - ^||vnf^ + (ivndv^)'' = 0, 



(8.25) 



+ d^(A(v„)) + = 0, B:= d^A. 

For superfiuid Yang-Mills magnetohydrodynamics the stress tensor is given by 

T := v„®m + p"®v^ + S-S + p5, 
where B-B is the (1,1) tensor field defined by 

(B-B)] 



and where the pressure is given by p := —e{p, S, — v„) + /ichcmP + ST — 

The corresponding Hamiltonian reduction and affine Lie-Poisson bracket can be found as 
before and the evolutions of the advected quantities are given by 

S = Jir]-')iSoor]-'), A = r]4Ad^Ao + XiTxl'), v, = (r/, (v^o + dx^ . 

As in the preceding example, it is possible to generalize this approach to multivelocity 
superfiuids. In this case, the 7-circulation gives 



di 



0, for all a 



m, and 



d 
dt 



A = 0, 



ct 



where q is a loop which moves with the normal fluid velocity v^. 
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8.6 Superfluid Hall Magnetohydrodynamics 

The Hamiltonian formulation of superfluid Hall magnetohydrodynamics is given in [iS?]. 
As one can guess, the Hamiltonian structure of these equations combines the Hamiltonian 
structures of Hall magnetohydrodynamics and of superfluids. This is still true at the group 
level and we will obtain the equations by affine Lie-Poisson reduction associated to the group 

G := [DiS{V)(i)J^{V,S^)] X Dm{V). 

In this expression, the symbol x denotes the direct product of the two groups. The advected 
quantities are 

{S,u;n) e J^{V) X X(I?) x J^{V). 

The variable S is the entropy density of the normal flow, the other variables will be inter- 
preted later. The action of (77, x; ^ ^ is given by 

(5, u; n) ^ (MS o r/), (r/*u^ + dx)h J^n o 0). 

The resulting affine Lie-Poisson equations consist of two systems, the affine Lie-Poisson 
equations associated to the variables (m, p, S, u) and the Lie-Poisson equations associated 
to the variables (n, n). 

The Hamiltonian of superfluid Hall magnetohydrodynamics is defined on the dual Lie 
algebra 

( [{X{V) ®J^{V)) 0) {T{V) © X{V))] X [X{V)®T{V)]y 
^ VL^{V) X T{V) X T{V) X Q}{V) X VL^{V) x T{V) 

and is given by 

^(m,p,^,u;n,n) := - \ j^P^^^W^ 1^ ^ (("^~ ' '^") ^ 



V 2 j-p 



where v„ is the velocity of the normal flow, := u — is the superfluid velocity, and e 
is the internal energy density. The one-form A is defined by 

n 

The norm in the first term is taken with respect to a fixed Riemannian metric g on V. The 
velocity v„ is the function v„ = v„(m, p, S", u; n, n) defined by the imphcit condition 

m - pv„ - —A = — (p, S, V, - v„) =: p. 

By the implicit function theorem, the above relation defines a unique function v„, provided 
the function e verifies the condition that the linear map 
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is bijective for all (r, s, v^, w^) eRxRxTV x TV. 
Using the notations 

Afchem := ^(p, S, Vs - v„) G T{V), T := — (p, S, - v„) e J^{V), 
de 

p:^—{p,S,Vs-Vn)en\V), 

the functional derivatives of h are computed to be 

5h 5h 1 ., ,,2 a 5h Sh 

dm dp 2 it dS du 

V := — = (apv„ + ap" + div By) , — = --A ■ v. 

on n ^ ' on R 

The vector field v is interpreted as the electron fluid velocity. The equations for p, 5, and n 

are given by 

d d d 

— p + div(pv„ + p«) = 0, — ,S + div(,5v„) = 0, and —n + div(nv) = 0. 

ot at at 

Using the expression of v in terms of v„ we obtain that div(nv) = — adiv(pv„ + p*) which 
proves that 

— {ap + n) = 0. 

Thus, if we assume that the initial conditions verify apo + uq — 0, then we have ap + n — 
for all time. The equations for A and u are computed to be 

w:^ = -iv„-B ipti-B i(divS)«-B 

OT p ^ ap ^ ' 

— u + grad U(vs,v„) + Pchem - 2l|vn|n + (iv„du'') = 0. 

Prom these two equations we obtain the evolution of the superfiuid velocity = u — ^A^ 
as 

= - grad g[Ys, Vn) + Pchem - 7,\Wn\ 



dt 

+ + ^i{divB)»^ - iv„dv^ 

Doing computations similar to those for superfluids we obtain that the equation for m + n 
is given by 

^(m + n) = -DivT, 

where the stress tensor T is given by 

T = v„ ® (pv^ + p) (^v^ + B-B + pS, p = ppchem + ST - E - ^\\Bf. 
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Thus, we have obtained the following equations 

d 



(m + n) = -DivT, 
p + div(pv„ + p») = 0, 



d_ 

di 



A 



dt 
d_ 

di 
d_ 

di 

= -grad I 5((Vs,V„) + /Uchcm - -\\^r. 



S + div(5v„) = 0, 



1 R 

-^v„B ipf-B i(divB)tt-B, 

p ap ^ ' 



(8.26) 



These are the equations for superfluid Hall magnetohydrodynamics as given in 19|] equations 
(35a)-(35e). When T> is three dimensional, the two last equations read 

= fv„ + -p« - — curlB^ X B and 

ot \ p ap J 



^v, = - grad ^5'(vs, v„) + /ichem - J II Vnf ) + v„ X curlv, 

+ - fcurlB - ^p" ) X B. 

p \ R 



Hamiltonian reduction for superfluid Hall magnetohydrodynamics. Consider the 
right-invariant Hamiltonian function if(m^, p^, S, u; n^, n) = if(5'_u;n)(m^, Px! induced by 
h and suppose that we have apo + uq = 0. A smooth curve 

(m^,p^;n5) G T* [{DiS{V)®J^{V,S')) x Diff(P)] 

is a solution of Hamilton's equations associated to H(^So,uo;no) ^i-nd with the initial condition 
Pq if and only if the curve 

(m, p; n) := (J(r')(m, o r^-^), J(r')(Px ° ^~'); '^(r')(n o T')) 

is a solution of the equations fl8.26l) . where = u — aA^/R = u — an/n. 

The Poisson bracket for superfluid Hall magnetohydrodynamics is the sum of the affine 
Lie-Poisson bracket associated to the variables (m, p, S, u) and the Lie- Poisson bracket as- 
sociated to the variables (n,n). 



The 7-circulation gives 



d 



0. 



Using the definition := u — -^A", we obtain 



dt 
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where q is a loop which moves with the normal fluid velocity v„. The Kelvin-Noether 
theorem associated to the variable n gives 




where dt is a loop which moves with the electron fluid velocity v. 



8.7 HVBK Dynamics for Superfluid ^He with Vortices 

The Hall-Vinen-Bekarevich-Khalatnikov (HVBK) equations describe superfluid Helium tur- 
bulence. We consider the version of HVBK equations, together with its Hamiltonian struc- 
ture, as given in 15|]. It turns out that this Hamiltonian structure is the same as that of 
superfluid Hall magnetohydrodynamics, that is, it is obtained by affine Lie-Poisson reduction 
associated to the group 

G := [DiS{V)(i)J^{V,S^)] X Diff(P). 

As before, the advected quantities are 

{S,n;n) G T{V) x X{V) x T{V), 

where S is the entropy density of the normal flow. The action of (r^, x; ^) G G is given by 

(5, u; n) ^ (MS o n), {r,*vt + dx)»; Ji{n o 0). 

The resulting affine Lie-Poisson equations consist of two systems, the affine Lie-Poisson 
equations associated to the variables (m, p, S, u) and the Lie-Poisson equations associated 
to the variables (n, n). 

For simplicity we assume that the manifold V is three dimensional. The Hamiltonian of 
HVBK dynamics is defined on the dual Lie algebra 

( [(X(P) (s)^(P)) (s) (^(P) © X \X{V)®J'{V)]y 

^ Q}{V) X T{V) X T{V) X 9}{V) X VL^{V) x T{V) 

and is given by 

/i(m, p, S", u; n, n) := j p||v„||V+ / (m - pA) ■ v„p + j £:(p, S", - v„, a;)p, 

where v„ is the velocity of the normal flow, := vl — A^ is the superfluid velocity, uj := curl 
is the superfluid vorticity, and e is the internal energy density. The one-form A is defined by 

n 

Note that the Hamiltonian of HVBK dynamics is similar to that of Hall magnetohydro- 
dynamics, for R = a = —1. The norm in the first term and the operator curl are taken 
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with respect to a fixed Riemannian metric g on V. The velocity v„ is the function v„ = 
v„(m, p, S, u; n, n) defined by the implicit condition 

de 

in - pvj, - = — (p, S, V, - v„, w) =: p. 

By the implicit function theorem, the above relation defines a unique function v„, provided 
the function s verifies the condition that the linear map 

is bijective, for all (r, s, v^, w^) G M x R x TP x TV. 
We make the following definitions 

ds Be 

p := — V, - v„,w) e A := — (p, v, - v„, e Q^(P). 

The interpretation of the quantities fichcm,T,p, and A is obtained from the following ther- 
modynamic derivative identity for the internal energy: 

d{e{p, S, V, - v„, u>)) = pchemdp + TdS + p • V. (v, - v„) + A • V.w e fl\V), 

where V denotes the Levi-Civita covariant derivative associated to the metric g. The func- 
tional derivatives of h are computed to be 

5h 5h 1 2 5h 5h 

7— =V„, = V„ - A ■ Vn + Pchem, Td^T, — = p CUrl A, 

dm dp 2 6S on 

Sh 1 . a , s\ . 

:= — = - (pv„ + p» + curl A) , — ^ A ■ vi. 
on n ^ 'on 

The vector field V; will be interpreted as the vortex line velocity. The equations for p, 5", and 

n are given by 

—p + div(pvn + p« + curl A) = 0, —S + div(5vn) = 0, and —n + div(nv/) = 0. 

Using the expression of in terms of v„ we obtain that div(nvi) = div(pv„ + + curl A) 
which proves that 

!<" -"'="■ 

Thus, if we assume that the initial conditions verify po = no, then we have p = n for all 
time. The equations for A and u are computed to be 

d 1 

—A^ = -(pvn + p" + curl A) X curl^l" = V; x curl A* 

+ grad ^^(Vs, v„) + Pchem - ^||Vn|P^ + curl U X v„ = 0. 
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As in superfluid dynamics, this equation preserves the condition curl u = 0, and we will 
suppose that it holds initially: curl uq = 0. In this case we have uj = — curl and the 
equations above read 



d_ 
di 



A« + V/ X Lj = and + grad (^^(v„ v„) + /iehem - ^l|v„lP^ = 0. (8.27) 



From these two equations we obtain that the evolution of the superfluid velocity = u — A" 
is given by 

^v, + u; X Vi = - grad |^fi'(v,,v„) + Hchem - ^||v„f^ , (8.28) 
which can be rewritten as 

+ Vv.v^ = - grad ^/Xchem - ^||vs - v„||^^ + f, f = (vi - Vs) X a;, 

(compare to (18.180 ). Doing computations similar to those for superfluids we obtain that the 
equation for m + n is given by 

^(m + n) = -DivT, 
where T is the HVBK stress tensor given by 

T = V„ ® (pV^ + p) + p« O - W ® + p = p/ichem + ST - E + LJ^ ■ \. 

Thus, we have obtained the following equations 

d 

-J = -D,vT, 

d d 

— p + div J = 0, —S + div(5v„) = 0, (8.29) 

+ Vv.Vg = - grad |^Pchem - ^||vs - v„f j + f, 

where f = (v; — v^) x a; and we have used the notation 

J := pv^ + p = m + n 

for the total momentum density. These are the equations of HBVK dynamics as given in 
equation (1) of [15] with R = 0. 

Hamiltonian reduction for HVBK dynamics. Consider the right-invariant Hamilto- 
nian if(m^, p^, S", u; n^, n) = /J(5u;n)(mr?, Px' induced by h. Suppose that po = Uq and 
curluo = and let {mr„p^;n^) be a smooth curve in T* [(DiS{V)®J^{V,S^)) x Diff(r')] 
solution of the Hamilton equations associated to -f^(5o,uo;no) and with the initial condition po. 
Then the curve 

(m, p; n) := (J(r/-^)(m, o r^-'), J{r'){Px ° r')\ AOi^ o C')) 
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is a solution of the HVBK equations (18.291) . where v>, = u — yl" = u + n/n. To obtain the 
converse of this assertion, it is not enough to assume that (m, p, S", v^, n) verifies fl8.29p . It 
is also required that A or u verify the corresponding equation in (18.27^ . 

Note that from the equalities = u — A" and curl u = we obtain that the variables 
u and are interpreted respectively as the potential and the rotational components of the 
superfluid velocity v^. 



As in ^8.61 we have 



where q and dt are loops which move with the normal fluid velocity v„ and the vortex line 
velocity v;, respectively. Using the equation (18.281) for the evolution of the superfluid velocity 
Vs, we obtain the vortex Kelvin theorem 

d 
dt 



vt = 0. 



dt 



By the Stokes theorem, this can be rewritten as 

I /■/■ (a;^n)ci5 = 0. 



St 



where n is the unit vector normal to the surface St whose boundary dSt is a loop which 
moves with the vortex line velocity v;. This is the conservation of the flux of superfluid 
vorticity through any surface whose boundary moves with the velocity v;. 



In 15| it is also supposed that |^ is coUinear to r, more precisely, that there exists a 



positive function ps such that 

P = Ps(Vs - V„). 

The function ps is interpreted as the superfluid mass density, and the density of the normal 
fluid is given hj pn '■= p — Ps- Using these notations, the total momentum J and the stress 
tensor T can be rewritten as 

J = PnVn + Ps^s 

and 



8.8 Classical Fluids Versus Superfluids 

We summarize below the examples that have been studied in the previous sections. 

Classical fluids. In order to compare the two theories, we use the notation v„ for the fluid 
velocity. We also express the dynamics in terms of the entropy density S and internal energy 
density e (and not in terms of the specific entropy s and internal energy e). We have 

6{p, S) = pe{p, S/p), 
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and the first law of tliermodynamics 

de = 4dp + rds, P^p'^ 

op 

reads 

ds = Pchemdp + TdS, p = ppchem + ST -E. 

(i) Basic hydrodynamics 

Symmetry group Diff(D), momentum m G Q^{T>) 
Advected quantities (p, S) G J^{V) x J^{V) 
Hamiltonian 



h{m,p,S) = l- / -||m||V+ / £{p,S)p, 
^ Jv P Jv 



m = pv^. 

Stress tensor formulation lii = — DivT, where 

T ^Vn0 pvl+pS, p = ppchem + -ST - £. 

(ii) Yang-Mills magnetohydrodynamics 

Symmetry group Diff (P) ®J^{V, O) 

Momenta (m, Q) G n\V) x J='{V, o*) 

Advected quantities {p,S,A) G J^{V) x T{V) x n^{V,o) 

Hamiltonian 

h{m,Q,p,S,A)^l f -||m||V+ / e{p,S)p+l f \\d^A\fp, 
^ Jv P Jv ^ Jv 

m = pv^. 

Stress tensor formulation rii = — DivT, where 

T = V„®pV^ + S-S+p5, p = pAtchem + '5r-£- 

(iii) Hall magnetohydrodynamics 

Symmetry group Diff(r') x Diff(D), momenta (m, n) G ^^(1)) x ^^(1)) 

Advected quantities (p, S\ n) G x J^{V) x 

Hamiltonian 

h{m,p,S;n,n) = ]- [ - m-^A p + [ e{p,S)p+l- [ \\dAfp, 
Jv P ^ Jv Jv 

where A :^ R-, 

t. 0,P A 

m = pv^ + -^A. 

Initial conditions 

apo + no = =^ {ap + n = and m + n = pv^} . 
Stress tensor formulation rii + h = — Div T, where 
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Superfluids. As before, we denote by v„ the velocity of the normal flow and by the 
superfluid velocity. For simplicity we treat the two fluid model in this summary. Gener- 
alization to multi fluids models follows as above in the examples. The superfluid first law 
reads 

d£ = Hchemdp + TdS + p • V. (Vs - v„) , 

where e is the internal energy density. For HBVK dynamics the term A • V_a; has to be 
added. 

(i) Basic superfluid hydrodynamics 

Symmetry group Diff (P) (s) S"^), momenta (m,p) e n^{V) x J='{V) 

Advected quantities {S,Vs) G J^iV) x X{V). 

Hamiltonian 

h{ni,p,S,Vs) ^ ~ / p||v„||V+ / (ni-v„)/x+ / e{p, S,v, - Vn)l^, 
^ Jv Jv Jv 

m = pv^ + p. 

Stress tensor formulation rri = — DivT, where 

T = V„ {pvl + p) + p" (g) p = p//chem + ST - S. 

(ii) Superfluid Yang-Mills magnetohydrodynamics 

Symmetry group Diff (P) © {^{V, O) x T{V, S^)) 
Momenta (m,Q,p) e Vt^{V) x J^{V, o*) x J^{V) 
Advected quantities (^,^,v,) e T{V) x Q}{V,o) x X{V) 
Hamiltonian 

/i(m, Q, p, S,A,Ws) 

Pl|vn||V+ / (m-v„)A(+ I e{p,S,Vs-Vn)l^+l I ||d^^||V, 
^ Jv Jv Jv ^ Jv 

in = Pv^ + P- 
Stress tensor formulation rii = — DivT, where 

T = V„ (pv^ + p) + p« ® + 5-5 p = p/iehem + ST - E - ^\\Bf. 

(iii) Superfluid Hall magnetohydrodynamics 

Symmetry group [BiS{V)(s)T{V, S^)] x Diff(P) 
Momenta (m,p,n) e n^{V) x n^{V) x J^('D) 
Advected quantities (5, u; n) e J^{V) x xfp) x T{V) 
Hamiltonian 

h{m,p,S,u;n,n) = ~ ^ ^ '^H'^"!!^'" + ^ (("^^ ^'^) ' '^") 

+ / e{p,S,\s -^n)p+7; / l|d-4||V, where 
Jv ^ Jv 

V, :=u--A», 
n R 
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m = pK + ^a + p. 



Initial conditions 



apo + ^0 = =^ [ap + n = and m + n = pv^ + p} , 
curl uo = =^ curl u = 0. 

Stress tensor formulation m + ri = — Div T, where 

T = v„ ® (pv^, + p) + p« ® + B-B + p6, p = p^chcm + ST-6- ^\\Bf. 

(iv) HBVK hydrodynamics 

Symmetry group [Diff (P) (s) JP'(P, S'^)] x Diff(r') 
Momenta (m,p,n) G n^{V) x ^^(D) x J^('D) 
Advected quantities u; n) G J^(I?) x xfp) x J^{V) 
Hamiltonian 

/i(m, p, S", u; n, n) = - ^ / p||vn||V+ / ((m - pA) ■ v„)/x 
+ / £(p, S',v^ - v„,a;)/i, 



A := , Vo := u — AK := curlv. 



m = pvl + pA + p. 



Initial conditions 



Po = {p = ''^ and m + n = pv^ + p =: j} , 

curl Uo = =^ curl u = 0. 

Stress tensor formulation rh + h = — Div T, where 

T = v„ O (pv^ + p) + p^ ® - a; (g) + p5, p = ppchem + ST - e + u;^ ■ X. 

The hypothesis p = Psiys — v„) implies the equalities J = p„v„ + p^v^ and 

T = v„ ® p„v^ + V, ® p,v^ + ® + 

8.9 Volovik-Dotsenko Theory of Spin Glasses 



In |32|, the authors use a Poisson bracket approach to derive the equations of nonplanar 
magnet with disclinations and of spin glass. These models are referred to as the Volovik- 
Dotsenko spin glasses. In 2^, the Hamiltonian structure of the Volovik-Dotsenko spin 



glasses is shown to be isomorphic to that of Yang-Mills magnetohydro dynamics. Thus, it 
can be obtained using the affine Lie-Poisson reduction developed in the present paper. In 



this section we also carry out the Lagrangian version of the approach given in (20| , to which 
we refer for additional comments on the physics of spin glasses. 
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The advected variables are 



peV* = T{V) and i ^V; = n\V, o). 

The variable p is the defect inertial-mass density and the curvature B = is interpreted 
as the disclination density. 

The reduced Lagrangian / : [X{D) ®J^{V, 0)] (S) [V^ © V^*] ^ M is given by 

/(u,z/,p,7) = i [ p||u||V + ^ / ||7(u) + z/||V-^ / PhllV (8.30) 

where the norms are associated to the metrics g, k, and (gk), respectively, and e is the con- 
stant of susceptibility. The unreduced Lagrangian L is the right-invariant function induced 
by I on the cotangent bundle. Since it has a complicated expression, we do not give the 
formula for L. We will justify the choice of this Lagrangian by showing that its Legendre 
transformation yields the Hamiltonian of the Volovik-Dotsenko spin glasses. 
This Lagrangian appears as a generalization of the Lagrangian 

/5g(^,7) = I ^ ll^^ll V - ^l^phfp, (8.31) 

associated to the macroscopic description of spin glasses in Q, where p is the constant of 
rigidity, see §8.1[ In order to understand mathematically the passage form the Lagrangian 
fl8.3ip to the Lagrangian fl8.30p . we consider the following general situation. 

General Case. Consider a Lagrangian Ip : JF(P, 0) ©^^(P, 0) ^ M, Ip = lp{v, 7) associated 
to a spin system and depending on a parameter p interpreted as the spin rigidity. Recall 
that the affine Euler-Poincare and advection equations are (see 



|f^ = -ad:f^ + div^f^, 

1:7 + = 0. 
ot 

To Ip we associate the Lagrangian / : [£{V)@T{V, 0)] (S) [TiV) © Vt^iV, 0)] M given by 

^(u, ^.P,l) Pll"llV + ^p(i^>7), 

where v := 7(u) + v and p is now a variable. Remark the analogy with the process of 
minimal coupling. In the case where Ip is given by (I8.3ip . we recover the Volovik-Dotsenko 
Lagrangian (18.301) . 

The functional derivatives of / are computed to be 

/€:= — = e T(V, 0* , m := — = pu^ + -^-7 G Vl\V), 
ov ov ou ov 

and 
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The advection equations are 

d d _ 

—p + div(pu) = and —7 + iu-B + d^u = 0. 

We now compute the affine Euler-Poincare equations. We have 

^ ?^ = - ad: - div r ?^uU div- f^A^+^A 



dt 5v ^ 5v \5v J \ 5v S'j 

- ad*-,_^(„^ ^ - div I I + div | 



-T.(ad;|u)..v.(| 
Using the equations for p, 7, the relation m = pu^ + ^ ■ 7, and the identity 



we obtain 



d ^ " 

— m + i'um + (div u)m 
ot 

= p + V.u + Vu^-u) - ad; ^-7 + div^ (f^) "7 - ^"iu^ 

-^■d^z^ + d^(u)-7+^-£u7 + (divu)^-7 
OU OV Of ou 

.J|„,V„u.VuO-ad;|.-, + d.v^a).-, 



|.d^Wu)-.) + d.v(|u).,, 



where in the last equality we used the identity fl3.12p and the definition of the curvature B. 
We also have 

51 V T 
— op = p Vu ■ u + grad — 

.^P J V h. 

and 

51 ( Sip 6lA f 5lp Sip 

— oi 7 = div^ -^u + 7^ -7 - T^u + ■i_B. 
07 \ou 0'-) J \ov 07 
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Thus the affine Euler-Poincare equation for the variable u is 



l-u + Vuu + Vu^-u^ - ad; + div 



Slf> 
^7 



■7 



+ ^-d^(7(u)-P) + div f^u 

Of \ Of 



■1 



sir 



6L 



-^■dz/ + p Vu' -u + grad-^ + div^ T^u + -7 



Sp 

Sip Sip . 
oz/ 07 



5L 



5v S'j 



which, after remarkable cancellations, reads: 



d „ Jlp 
u + VuU = grad -r^ — 



1 fSL-. 



-u + 



5L 



■i B. 



dt ' ' " ° P \Siy ' ^7 

Thus the affine Euler-Poincare equations associated to / are given by 



— u + VuU = grad 

at op p 



d Sip 6lp , ^ / Sip 



sip 6lp , 
7^u + 1 -kB 

ou 07 



d_ 

dt 



7 + iu5 + d^z/ = 0, 



(9 

—p + div(pu) = 0, 



B = d^7, 



z/ = 7(u) + z/. 



(8.33) 



This system of equations can be seen as a generalization of equations ( ]8.32p . If Ip is 
hyperregular, then / is also hyperregular, and the associated Hamiltonian is given by 



Jv 2p 



(8.34) 



where hp is the Legendre transformation of I p. 

The case of Volovik-Dotsenko Spin Glasses. We now particularize the previous dis- 
cussion to the case of the Volovik-Dotsenko spin glasses, that is, we consider the Lagrangian 
f l8.30p . In this case, the momenta are given by 



K = e(7(u) + z/) and m = pu^ + n-'y, 



and equations (18.331) read 



— u + VuU = -grad-||7| 

^K + div^(p7») = 0, 
d 

—p + div(pu) = 0, 



7" - -Ku 1 i B 



1-7 + iu5 + -d^«tt = 0, 
ot e 



B = d^7. 



(8.35) 
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In the case of the spin glass dynamics described by the equations above, the term — grad ^||7||^ 
is an analogue of the electrostatic force and the term ^7^ — ^^^uj i_B is an analogue of the 
Lorentz force. 

Lagrangian reduction for the Volovik-Dotsenko spin glasses. A curve (77, x) ^ 
Diff ("D) (S)JF(D, O) is a solution of the Euler-Lagrange equations associated to the Lagrangian 
if and only if the curve 

(u, u) := (?) o 7]-^, {TR^~ix) o 7]-^) 

is a solution of the spin glasses equations fl8.35p with initial conditions (po,7o), and where 
K, = e(7(u) + u). 

The evolution of the advected quantities is given by 

p = J(ry-i)(po or^-i) and 7 = r/, (Ad^ 70 + , 

and the evolution of the disclination density is 

B = 7], (Ad^d^°7o). 

The variable x is the orientation of Lagrangian particles in their reference configuration. 

By Legendre transforming i^(po,7o) and /, we obtain the right-invariant Hamiltonian if(p(, 
and the reduced Hamiltonian h, given by 

h(ni, fi:,p,7) = ^ j ^\\m - K--ff fx + ^ j j PhllV, 

see also fl8.34p . Note that this Hamiltonian differs by a sign in the first term from the 
Hamiltonian (2.26b) in 20|. This is due to our convention in the Hamiltonian structure 
which also differs from theirs. This justifies the choice we made for the Volovik-Dotsenko 
Lagrangian fl8.3Up . 

Hamiltonian reduction for the Volovik-Dotsenko spin glasses. A curve (m^, /t^) G 
T*[Diff (P) (S)JF(D, C)] is a solution of Hamilton's equations associated to the spin glasses 
Hamiltonian if(p(,^^y) if and only if the curve 

(pu^ + K-7, k) =: (m, k) := J(?7~^) (m^ o r/"^ T*i?^o^-i (k^ o r^"^)) 

is a solution of the system (18.351) with initial conditions (po;7o)- 

The associated Poisson bracket is identical to that of Yang-Mills magnetohydrodynamics, 
except the fact that the variable s is not present in this case. Notice that the variables have 
not the same physical meaning in the two theories; see 201 . 

The Kelvin-Noether theorem gives the complicated expression 

— (/" ( u'' H — 7 ) = ([ - (—Kdu + div'''(w)7 — w ■ i_B) , 
dt /cA P J Jet P 



where 

u := 

e 
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— — 7(u), w := Ku — pj^ 



which can be rewritten in the simpler form 



d C u /■ / u 1 



The 7-circulation gives 



Ct Ct 



8.10 Microfluids 

Microfluids are fluids whose material points are small deformable particles. Examples of 
microfluids include liquid crystals, blood, polymer melts, bubbly fluids, suspensions with de- 
formable particles, biological fluids. In this section we find the Hamiltonian structure of the 
equations governing the motion of non-dissipative microfluids in Eringen's formulation by 
showing that they appear b y E uler-Poincare and Lie-Poisson reduction. 



We quickly recall from [Hj some needed facts about microfluids. A material particle 
P in the fluid is characterized by its position X and by a vector S attached to P that 
denotes the orientation and intrinsic deformation of P. Both X and S have their own 
motions, X x = ri{X, t) and S ^ = S, t), called respectively the macromotion and 
micromotion. Since the material particles are considered to be of very small size, a linear 
approximation in S is permissible for the micromotion. Therefore, we can write 

e = x(^,t)s, 

where x{X,t) G GL(3)+ := {A G GL(3) | det(A) > 0}. The classical Eringen theory 
considers only three possible groups in the description of the micromotion of the particles: 
GL{?))^ D CSO{?>) D 5*0(3). These cases correspond to micromorphic, microstretch, and 
micropolar fluids, respectively, all of them discussed in detail below. The Lie group CS0{3) 
is a certain closed subgroup of GL{3)~^ that is associated to rotations and stretch. Of course, 
the general theory developed in this paper admits other groups describing the micromotion. 

Micromorphic fluids. A fluid in a domain V is called micromorphic if its macromo- 
tion and micromotion are described respectively by a diffeomorphism r] G Diff('D) and a 
function x ^ J-'{'D,GL{3)~^). As a consequence of this description, the configuration man- 
ifold of micromorphic fluids is isomorphic to the product of the two groups Diff(P) and 
J^{V,GL{3)^). We will show the remarkable facts that, as in the previous examples, the 
relevant group structure on the configuration manifold is given by the semidirect product 
DiS{V)(i)J^{V,GL{3)+). 

The general equations for micromorphic continua are given by 

tki,k + p(^fi- ^u,^ = 0, 

^klm,k + tml — Sml + p{llm — Im) = 0, (8.36) 

D , D 

—p + p dlV U = 0, —%kl - tkrJ^lr - Hrl^kr = 0, 

dt at 
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where D/dt denotes the material derivative. 

Since this is the first time this symbol D/dt appears, it is useful to make some comments 
and to place this discussion in the larger context of complex fluid motion on general Rieman- 
nian manifolds. The domain "D C is endowed with the usual Riemannian metric given 
by the inner product in Euclidean space. This Riemannian metric has an associated Levi- 
Civita connection V whose covariant derivative along an arbitrary vector field u is defined 
by Vu := u ■ V. Of course, if T) is replaced by a general Riemannian manifold, then Vu is 
the covariant derivative of the Levi-Civita connection associated with the given metric and 
does not have this simple expression. The material derivative of a vector field v along u is 
given, in general, by 

D d ^ 
-V = -V + V.v. 

On functions f G J-'(T>, V), where is a vector space, the material derivative is given by 

which is metric independent. 

In (18.361) the material derivative operator is applied to the vector field u and to the 
functions p and iki, even though i^i are the components of the symmetric tensor i. Thus, 
one should not interpret the last equation in fl8.36l) as the covariant derivative of the tensor 
i. This convention is in force throughout the rest of the paper and we shall comment on this 
in the relevant places. 



Equations fl836D are given in equations (2.2.10), (2.2.11), (2.2.31), and (2.2.32). The 
variable p is the mass density, Vki is the microgyration tensor, and the symmetric tensor iki 
is the microinertia. The vector field fk is the body force density and km is the body couple 
density. The spin inertia aki is given by 

_ • 

Ckl — Iml I '^^krn + l^knl^n 

The tensors tki,mkim, and Ski are respectively the stress tensor, the couple stress tensor, 
and the microstress tensor; see chapters 1, 2, and 17 in [ll] for details. The four equations 
in fl8.36p correspond respectively to the balance of momentum, the balance of momentum 
moments, the conservation of mass, and the conservation of microinertia. 

The constitutive equations for non-dissipative micromorphic fluids are given by 

tkl = 7^ — T^ki, mkim = 0, and Ski = -7; — r^z + 2p— — v/, (8.37) 
op~^ op~^ Olrk 

see (3.4.6) and (3.4.7) in [ll|. Here the function ^ = ^(p"\ z) : M x Sym{3) M is the free 
energy and Sym{3) denotes the space of symmetric 3x3 matrices. It is usually assumed 
that 

^{p-\ A-HA) = ^{p-\i), for all A e 0(3). 

This condition is imposed by the axiom of objectivity. Since i is symmetric, this implies that 
the free energy depends on i only through the quantities Tr(i),Tr(z^) and Tr(i^). 
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Using the constitutive equations (18.371) and assuming that fk 
(KWf read 

D . 

OlT{ 7, CTZm — —27- irl, 



0,lki = 0, equations 



D , D 
— p + pdivu = 0, —iki 
at at 



dirm 



(8.38) 



0. 



These are the equations for non-dissipative micromorphic fluids. 

Microstretch fluids. A microstretch fluid is a micromorphic fluid whose micromotion x 
takes values in the four dimensional real Lie group CS0{3) defined by 

CS0{3) = {Ae GL{3y I there exists A G M such that AA^ = A/3} 

and called the conformal special orthogonal group. The Lie algebra of (75*0(3) is 

C3o(3) = {u E qI{3) I there exists p G M such that u + = /i/3} . 

Note that in the relations above we necessarily have A^ = det(A)^ and 3/i = 2Tr(z/). Note 
also that each u G cso(3) decomposes uniquely as 



(8.39) 

The equality (I8.39p can be rewritten 



where z/q G M, and (t'l, 1^2, t^s) =: is a vector in \ 
as 

V = z/0/3 + 9, 

where u G so (3) is the matrix whose entries are given by Dij = —Sijkfk- 

The material particles of microstretch fluids have seven degrees of freedom: three for 

translations, three for rotations, and one for stretch. 

The general equations for microstretch continua are given by 



tkl,k + PUl - 



0, 



'mk,k + t - s + p{l - a) = 0, mki,k + eimntmn + p{h - = 0, 

^p + pdivu = 0, 

at 

D D 
-rjo - 2jot'o = 0, —Jki - 2uojki 



(8.40) 



dt 



dr 



These equations are given pjj], equations (2.2.38) to (2.2.41). Following the conventions 
introduced in the discussion following (]8.36p . the material derivative operator D/dt acts on 
the vector field u and on the functions uo, vi, p, jo, and jki. As before, the variable p is 
the mass density. The variables and u are respectively the microstretch rate and the 
microrotation rate, constructed from Uki as in (18.391) . The microstretch microinertia jo aiid 
the microinertia jki are constructed from iki as follows: 

jo := 2zfcA: and jki := ^jo4i - ki- 
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The microstretch spin inertia a and the spin inertia are given by 

D , 1 2 



JklT^k^l- 



The microstretch vector m^, the couple stress tensor rriM, the microstretch force density I, 
and the couple density I are defined from rrikim and Iki by the decompositions 



1 



mklr. 



Im 



rSimrrrikr and / 



kl 



kl 



Z^klr^r- 



We also used the notations t := t^k and s := s^k- Using these definitions and the fact that the 
fiuid has the microstretch property, one can obtain the equations f l8.40p from the equations 



The constitutive equations for non-dissipative microstretch fiuids are given by 
(9^ 



tkl 



dp' 



-bkh ruki = 0, mk = 0, and s - t = 2p [ ^—jki + 77-jo 

\d3ki djo J 



where \1/ = \E'(p ^,j,jo) is the free energy, see (3.4.16) in 

Micropolar fiuids. A micropolar fluid is a micromorphic fiuid whose micromotion x takes 
values in the Lie group 50 (3). The material particles of micropolar fiuids have six degrees 
of freedom: three for translations and three for rotations. Micropolar fiuids are, therefore, a 
particular case of microstretch fiuids. 

The general equations for micropolar continua are given by 



tkl,k + P I // 



D 
di 



0, 



nT'kl,k + £lmntmn + p{k — (^l) — 0, 

D 

Jt' 
D 

dt- 



iP + p div u = 0, 
I 

:jkl + i^kprjlp + £lprjkp)l^r = 0, 



(8.41) 



where the variables and tensors have the same meaning as in the previous examples. These 
equations are given [nl, equations (2.2.43) to (2.2.46). They are derived from those of mi- 
crostretch continua, using that the microgyration tensor u takes values in the Lie subalgebra 
50 (3) of cso(3), that is, = 0. In particular, the spin inertia au is given by 

. D , D.. ^ 

Cfe 3kl-^1-'l + ^klmjmn^l^n — -^[.Jkll^lj- 



The constitutive equations for non-dissipative micropolar fiuids are given by 

tkl 



—6ki and = 0, 



dp 



where \1/ = \E'(p is the free energy, see (3.4.27) in 11 
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Lagrangian and Hamiltonian formulation for micromorphic fluids. We now show 
that the equations for non-dissipative micromorphic fluids can be obtained by Euler-Poincare 
reduction associated to the semidirect product Diff(T') (§)J^(V, GL{3)~^). The advected quan- 
tities are the mass density p G and the microinertia tensor i G J-'{V, Sym{3)). The 
symmetry group DifT(P) (S)JF(P, G'L(3)^) acts linearly on the advected variables {p,i) by 

(p, i) ^^ {Jri{p o r]), o v)x) ■ 

The choice of this group representation is dictated by the form of the advection equations 
for the mass density and the microinertia in equations ( 18. 38^ . The infinitesimal actions 
p(u, z/) = pu and i(u, i/) of (u, i/) G X(V) (§) (V , qI{3)) on p and i are given, respectively, 
by 

pu = div(pu) and i(u, u) = di(u) + u^i + iu. 

Given two matrices a,b ^ 0^(3) we denote by ab their product and by a-b := Tr(a^b) the 
contraction. By identifying the dual of J-'{T>, Sym{3)) with itself through the pairing 



{m,i) = / m{x)-i{x)p, 
Jv 

we obtain the diamond operations 

mOii = —m-di and m02i = —2im. 

In order to obtain the equations for micromorphic fiuids we consider the Lagrangian / : 
[XiV)®J^{V, gl(3))] (S) [TiV) © T{V, Sym{3))] ^ M given by 

/(u, z/,p,i) := ^ / p||ufp+^ / p{iy-y)p- I p^{p~^,i)p, (8.42) 
Jv ^ Jv JV 

where the function \& represents the jree energy^ and where the norm in the first term is 
taken relative to a fixed Riemannian metric g on D. The functional derivatives are 

bl ^ bl . 5/ 1,, ,,2 1 . ^, _i \ , ., 
ou dz/ dp 2 2 pop ^ 

and - = -puu -p— (p ,.). 

A computation, involving remarkable cancellations, shows that the Euler-Poincare equations 
(fT:2|) associated to the Lagrangian (18:421) and to the group G = DiS{V) ® T{V, GL{3)+) 
acting linearly on V* = J^(V) © J^(T>, Sym{3)) are given by 

f d 1 
— u + VuU = - grad 



dt p dp~^ ' 

f d \ 
' [di'' + '^''^''^ -^^j= 2^^, (8.43) 

d d 

—p + div(pu) = 0, —i + dz(u) + u^i + iu = Q. 

C/6 kJL 



74 



Thus we have recovered the equations (18.361) by Euler-Poincare reduction, up to a change of 
variables, replacing v by —v^. 

Consider the right-invariant Lagrangian i^(po,jo) on T [Diff ("D) (s) JF(I?, GL(3)+)] induced 
by the Lagrangian (ElSl). A curve (r/, x) e Diff (P) (s) J^(r', GL(3)+) is a solution of the 
Euler-Lagrange equations associated to i^(po,io) if and only if the curve 

(u, v) := {fj o r]-\ XX'' o V~') ^ X(P) ©^(D, GL(3)+) 

is a solution of the micromorphic fluid equations (18.431) with initial conditions {po,io)- Note 
that, in our approach, the relation u = (xx~^) ° V^^ between the microgyration tensor 
and the micro and macromotions is due to the passage from the Lagrangian to the spatial 
representation. This relation coincides with that given in [llj], up to conventions. It is 
considered there as a definition of u. 

The evolution of the mass density p and the microinertia i is given by 

P = J{v~^){po°V~^) and i = (^(x^) ^ iox"^^ ^ ■ 
Remark that the evolution of the determinant of i is 

Therefore, if the initial microinertia iq is invertible, then i is invertible for all time. Under 
this hypothesis we can take the Legendre transformation of the Lagrangian / and we obtain 
the Hamiltonian 

h{m,K,p,i) = ]- / -||m||V + ^ [ -{r^K-K)p+ ! pm{p~^,i)p, 
^ Jv P Jv P Jv 

where {m,K,p,i) G Q^{V) x T{V,gi{3)) x JF(T') x J^{V,Sym{?,)), which consists of the 
sum of the kinetic energy due to macromotion and micromotion and the free energy. By 
Lie-Poisson reduction, we obtain that the equations (18.431) are Hamiltonian with respect to 
the Lie-Poisson bracket 



{/,^}(m,K,p,i) = / m 
Jv 



(5m' 6m 



P 



+ n - ( adsi — + d-— • d— • -- — p 

OK OK dm OK dm 



V 



+ /P|d|^)#^-d(^)l^l^ (8.44) 



V 



+ i\ div 



6p J 6m \Sp J 
6f6g\ 6g6f 6ff6g^^ 



V 



6i 6m J 6k 6i 6i \ 6k 



_ div {^-i^] + + ^ p 

6i 6m 1 6k 6i 6i \6k ) 
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The Kelvin- Noether circulation theorem applied to micromorphic fluids yields the simple 
relation ^ 

dtfct^ jet di 

Lagrangian and Hamiltonian formulation for microstretch fluids. The symmetry 
group of microstretch fluid dynamics is DiS(V) (§)J^{V,CS0{3)). As before, the advected 
quantities are the mass density p G J-'(V) and the microinertia tensor i G J-'(V, Sym{3)), 
on which the symmetry group acts linearly as in the micromorphic case. 

The Lagrangian of the microstretch fluid has the same expression as that of the micro- 
morphic fluid, namely, 

/ : [X{V)(§)J^{V, cso(3))] (S) [J^{V) © J^{V, Sym{3))] M 

is given by 

/(u,z/,p,i) := - / p\\ufp+- / p{iv-v)p- I p'^{p-\i)p. 
Jv Jv Jv 

Using the change of variables i ^ j := Tr(i)/3 — i, this Lagrangian reads 

l(u,jy,p,j) = ^ j^p\\ufp+^ j^p(^^jou^+ju-v'^ p- j^p^{p-\j)p, (8.45) 

where jo := Tr(j). The functions uq G J-'(V) and u = {ui, 1^2, 1^3) G J-'iV, M.^) are deflned as 
in the equation ( 18.39^ . The expression ( 18.45^ for the Lagrangian has the advantage of giving 
the expressions of the energy due to microstretch and microrotation separately. 

The associated Euler-Poincare equations are the same as in (18.431) . We now rewrite these 
equations, using the fact that, for microstretch fluids, the microgyration tensor u has values 
in the Lie subalgebra cso(3). We will make use of the following two lemmas. 

Lemma 8.1 Suppose that the free energy \E' verifies the axiom of objectivity, that is, 

^{p-\A-HA) = ^{p-\i), for allAe 0(3). 

Then the matrix is symmetric. 

Proof. Consider a curve A{t) G ^0(3) such that A{0) = h and i(0) = ^ G so (3). 
Differentiating the equality {A{t)~^iA{t)) = \l/(z) at t = 0, we obtain the condition 

D^(i)(i^ - ^t) = 0, for all ^ G so (3). 

Using the equalities D\l/(i)(2^ ~ ^i) = Tr {^{i^ — ^i)) = Tr {{^i — i^) ^), we obtain that 
the matrix ^i — i^ is symmetric. Since this matrix is clearly also antisymmetric, we obtain 
^t-i^=0. ■ 

oi ai 

Lemma 8.2 Fori G J^{V, Sym{3)) define 

3:=i:T{{)h-ieJ'{V,Sym{3)). 
For V G TiT)., cso(3)) define Uq G J^{T>) and v G JF(X),M^) hy the condition 

V = z/0/3 + 9. 

Then, 
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the equation 



is equivalent to the system 



di 



(8.46) 



Jo fD 



2 \dt 

D , , 

j—v - 2uoju - [ju) X 1/ = 0, 
at 



djo 



dj J 



]A7) 



the equation 

is equivalent to the equation 



D 

di 



i + v i + iu = 



D 

dt 



J + 2uoj + [j, u] = 0. 



Proof. The results follow by direct computations. The two equations in (18.470 are obtained 
by taking respectively the trace and the antisymmetric part of the equation (18.461) . For the 
computation of the trace, we use the equalities 

1 1 
Tr(iz/) = z/oTr(z) = -z/Qjo, Ti{iuu) = -i'q - {ju)-u, and 



Tr I i 



9^ 

di 



(9* 5^ 



For the computation of the antisymmetric part, we use the equalities 



IV — \%V) = JU 



IVV — [IVV] 



2vQ[iv 



[iv] ) + [iuu 



[ivu] 



and the fact that the matrix i-^ is symmetric, by the preceding lemma. ■ 

Using this lemma, we obtain that the Euler-Poincare equations ( 18.431) become 

( d _ 1 3^ 
— u + VuU = - grad - — -, 
at p op~^ 



jo fD 



J-^ty - 2i/ojV - (jV) xu = 0, 



(8A8) 



d_ 
di 



p + div(pu) = 0, 



D 
~dt 



J + 2z/oj + [j, P] = 0. 



The third equation can be rewritten as 



D 
It 



{ju) = 0. 
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From the last equation, we deduce that the conservation law for jo reads 

D 

-nJo + St'ojo = 0. 
at 

Thus we have recovered the equations fl8.40l) by Euler-Poincare reduction, up to a change of 
variables, replacing u by —u'^. 

The Lagrangian reduction for microstretch fluids follows from that of micromorphic fluids 
and so we do not repeat it. Making use of the fact that the micromotion x takes values in 
CS0{3), we obtain that the evolution of the microinertia i with initial value io is 



O 7] ^. 



It follows that the evolution of the variable j has the same expression and that the evolution 
of the microstretch inertia jo, with initial value (jo)o5 is given by 

As in the micromorphic case, if the initial microinertia tensor io is invertible, we can obtain 
the equations and its associated Poisson bracket by Lie-Poisson reduction. 

Lagrangian and Hamiltonian formulation for micropolar fluids. The symmetry 
group of micropolar fluid dynamics is DiS (T>) (§) J-" (T>, SO {3)). In terms of the variable j, 
the Lagrangian reads 



l{u,u,p,j) = l [ p||u||V + ^ / p{jiy-iy)p- I p^(p \ 
Jv ^ Jv Jv 



j)p- 



Using that z/q = for micropolar fluids, we obtain that the Euler-Poincare equations are 
given by 

d ^ 1,9* D ^ ^ 

-u + V.u = - grad j-u- (ju) x . = 0, 

d D (8-49) 

—p + div(pu) = 0, -J + [j, ly] = 0. 

Note the analogy between the second equation and the equation for the rigid body. 

The Lagrangian reduction for micropolar fluids follows from that of micromorphic fluids 
so we shall not repeat it. Making use of the fact that the micromotion x takes values in 
5*0(3), we obtain that the evolution of the microinertia j with initial value jo is 

Thus, the evolution of its determinant is given by 

det(j) = det(jo) o r^"^ 

which shows that j is invertible if and only if jo is invertible. Thus, if the initial microin- 
ertia tensor jo is invertible, the equations of motion (I8.49P can be obtained by Lie-Poisson 
reduction. Using the equalities 

5/ , 51 . 

m := — = pu and «;:= — = pji/, 
ou ou 
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the Legendre transformation yields the Hamiltonian 



Jv P 



|m|| /i + 



1 

"2 Jv P 



representing the total energy of the system. In the particular case of micropolar fluids, and 
in terms of the variable j, the Poisson bracket fl8.44p becomes 



V 



{/,^}(m,K,p, j 
+ 

+ 

+ / J- 

V 



m 



6m' (5m 



ad*/ — + d— ■ d— P 

*t OK OK dm OK dm 



5.50) 



P 



V 



d 

div 



5f\ 5g 



6p J 6m 
6j 6m 



d 

+ 



Sl\ 61 
6p J 6m 
SI 6_g_ 

6j ' 6k 



P 
— div 



5l6l 
6j 6m 



5i 6_f 

6j ' 6k 



P, 



where the brackets in the last term denote the usual commutator bracket of matrices. 

A quaternionic point of view on microstretch and micropolar fluids. We now show 
that, in the case of microstretch fluids, we can use the Lie group of invertible quaternions 
to describe the micromotion of the particles. To see this, recall that there is a 2 to 1 surjective 
group homomorphism 

n:S^^ SU{2) — y S0{?>) 



given by 



a + j/3 



a —[3 
P a 



-2Re{ap) -2Im(a/3) 



2 Re (a/5) 
-2Im(a/3) 



where the universal covering group 5"^ = G H | |Q;p + |/5p = 1} C denotes 

the Lie group of unit quaternions. The tangent map of vr at the identity is the Lie algebra 
isomorphism 

p = k(j/i + iu2 +31^3) e TiS^ I — >2i>e so(3), 

where u = {ui, 1/2, 1^3) and the Lie algebra TiS^ consists of pure quaternions. 
Remarkably, the map vr extends to a 2 to 1 surjective group homomorphism 



TT : 



CS0{3) 



given by the same expression. Thus the group of invertible quaternions can be seen as 
the universal covering group of the conformal special orthogonal rotations. The Lie algebra 
isomorphism reads 

p = z/Q + k(i/i + ii^2 +31^3) e H I — ^2 (z/q/s + 0) e cso(3). 

These observations show that the micromotion of a microstretch fluid can be described by a 
map X : D — i> with values in the group of invertible quaternions. 
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The same remark applies to micropolar fluids where one can use the group of unit quater- 
nions instead of SO (3). In both cases, the only difference in the equations of motion is 
that the variable u takes values in EI or in the Lie algebra of pure quaternions, respectively. 

Remarks on the use of other groups and the anisotropic cases. As is clear from 
the general theory in [|3] and ^ the Lagrangian and Hamiltonian theory of microfiuids can 
be formulated for any order parameter subgroup O C GL{3)^ to model different kinds of 
micromotion. We have chosen two specific subgroups relevant to anisotropy. The first one, 
SO{K), was suggested to us by D. Holm (see [1^], Vol I, §1.10.2) and the second, CSO{K), 
is its conformal version. 

For K an invertible and symmetric 3x3 matrix, we define the group 

SO{K) = {Ae GL{3y I A^KA = K] . 

This is a three dimensional Lie group, called the special K -orthogonal group with Lie algebra 
given by 

so{K) = {z/ G 0[(3) I p'^K + fsTz/ = 0} . 

Note that SO{K) is the group of orthogonal linear maps relative to the (possibly indefinite) 
inner product u^Kv. 

Recall than when K = the linear isomorphism given by the hat map^ : —>■ so(3) 
pulls back the commutator bracket on so(3) to the cross product on M^. We shall search 
for the analogue of this map for 5o{K). For an arbitrary invertible symmetric matrix L we 
define : _ so{K) by 

:= R-^Lu e so{K). 

Working in this generality allows us later to make specific choices of L that simplify certain 
formulas. The adjoint and coadjoint actions of SO{K) on are computed to be 

Ad^ u = L'^ALu and Ad;^-i v = LA^^L'^v. 

By differentiating the adjoint action we find the expression of the Lie bracket on associated 
to the group SO{K). This Lie bracket generalizes the cross product u x v associated to 
SO (3), thus we shall use the notation u v and we call it the cross-product associated to 
K. We have 



u Xk V : = — 

dt 



d 
it 

L-^K-\Lu X Lv) = det{L)L'^K-^L-\u x v 



Ad /,'-K\ V = L ^u^ Lv = L ^LuLv 

expl tu 1 

t=0 



where we used the relation 

Lu X Lv = det{L)L~^{u x v), 

valid for any invertible and symmetric matrix L. To show this formula, we first note that it 
holds when L is a diagonal matrix. Then it suffices to write L = B^DB, where B G SO (3) 
and D is the diagonal matrix of eigenvalues. Note that, by construction, we have the formula 

= u ,v 
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The infinitesimal coadjoint action is computed to be 



ad^w = det{L){L-^K-^L-^w) x u. 



Note that so{K) = [5o{K),5o{K)], as is easily seen using the formula for x^- Thus, 
since dimso(ir) = 3, the Lie algebra so{K) is simple (see 2^, example at the end of ^2, 
Chapter I). According to the classification of all three dimensional real Lie algebras (see [24| . 
problems 28-35 in §15, Chapter I), it follows that if K is definite, then so{K) = so{3) and 
if K is indefinite, then so{K) = s[(2,M) (see also the discussion at the end of §14.6 in (j?!). 

A Casimir function for the Lie-Poisson bracket on (M^, Xk) is given by 



C{w) 



since we have 



ad|c/5^^ w = 0. 



This shows that the coadjoint orbits of SO{K) are ellipsoids (when K is definite) or hyper- 
boloids (when K is indefinite). 

Note that if we choose L = det{K)K~^ , the cross-product and the Casimir function take 
the simple expressions 

u Xk V = K{u X v) 



and 



C{w) = -w^Kw. 



Other relevant choices for L are L = I3 and L = K. 

When K is positive definite, the group O = SO{K) can be used for the description of the 
anisotropic version of the theory of micropolar fluids. In such micropolar fluids, the particles 
have preferred axes of deformations. Note that if -u G M'^ is given, then the 5'0(-ft')-orbit 
{Au I A e SO{K)} through u in coincides with the ellipsoid £c ■= {v E IRj^ \ v^Kv = c} 
containing u. Typical cases are given by 



K 



1 
1 

e 



or K 



e 
£ 
1 



, < £ < 1, 



corresponding to extremely oblate and prolate elhpsoids. 

As in the micropolar case, the Lagrangian of the anisotropic micropolar fluid is deduced 
from that of the micromorphic case. For the microrotation energy, we have identity 

Ti{{iu)^i^) = Jki^- u, 

where = v and where 

3k ■■= Tr {iK'^) - LR-HR-^L G Sym{3). 
If we choose L = R, then we obtain 

jK := Tr {tR'^) R^ -IE Sym{3), 
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which is the anisotropic generahzation of the relation 

j = Tr(z)/3 - i 

that was used in the micropolar case. Thus the Lagrangian of the anisotropic micropolar 
fluid can be written as 

^ Jv ^ Jv Jv 

In terms of the variables (u, z/, p, i) the associated motion equations are given by fl8.43p . 
except that now v takes values in so{K). One can write these equations in terms of the 
variables (u, v^p^j^) by simply replacing z/ by and i by 

^ Tr(j^L-2)ir2 - KL~^jkL'^K. 

in the micromorphic fluid's equations (18.431) . 

We can also consider an anisotropic version of the microstretch theory, by choosing as 
internal symmetry group the Lie group 

CSO{K) = {Ae GL+{3) I there exists A G M such that A^KA = XK} 

called the conformal special K-orthogonal group. The Lie algebra of CSO{K) is 

cso{K) = {z/ G qI{3) I there exists p G M such that u^K + Ku = pK}. 

In such anisotropic microstetch fluids, in addition to have preferred axes of deformations, 
the fluid particles can also stretch. 

As in the microstretch case, we have, for all z/ G cso{K), the decomposition 

z/ = + Z/0/3, 

for a unique 1/ G and z/q G M. The adjoint action and Lie bracket read 

AdA{i>^ + i^oh) = ^^^(]^y/3 [L^^ALu] ^ + z/0/3 and [z/, p] = i/x]^^^ 

As in the microstretch case, the Lagrangian of the anisotropic microstretch fluid is deduced 
from that of the micromorphic case. In the anisotropic case the microrotation energy takes 
the complicate expression 

Tr((2z/)^z/) = JkJ^ u- 2z/oTr {{KL-'jkL-'K)u'') 

If we choose L = K, then this equality reads 

Triizufu) = jKU-u - 2Z/0 Tr(jxi>'') + z^o'^ MjkK-') Tt{K') - ul Tr(jx), 
which generalizes the formula 

Tr((iz/)'^z/) = + ^JqUq 

valid when K = h. 
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8.11 Liquid Crystals 



The liquid crystal state is a distinct phase of matter observed between the crystalline (solid) 
and isotropic (liquid) states. There are three main types of liquid crystal states, depending 
upon the amount of order in the material. The nematic liquid crystal phase is characterized 
by rod-like molecules that have no positional order but tend to point in the same direction. 
In cholesteric (or chiral nematic), molecules resemble helical springs, which may have op- 
posite chiralities. As for nematics, the molecules exhibit a privileged direction, which is the 
axis of the helices. Smectic liquid crystals are essentially different form both nematics and 
cholesterics, in that they have one more degree of orientational order. Smectics generally 
form layers within which there is a loss of positional order, while orientational order is still 
preserved. See for example 0], [3], and 33] for more information. 

There are various approaches to the dynamics of liquid crystals. In this section we carry 
out the Lagrangian an Hamiltonian formulation for three of them: 

- the director theory due to Oseen, Frank, Zocher, Ericksen and Leslie, 

- the micropolar and microstretch theories, due to Eringen, which take into account the 
microinertia of the particles and which is applicable, for example, to liquid crystal polymers, 

- the ordered micropolar approach, due to Lhuillier and Rey, which combines the director 
theory with the micropolar models. 

For simplicity we suppose that the fluid container P is a domain in and all boundary 
conditions are ignored. This means that in all integration by parts we assume that the 
boundary terms vanish. 



8.11.1 Director theory 

In this theory it is assumed that only the direction and not the sense of the molecules 
matter in the description of the physical phenomena. Thus, the preferred orientation of the 
molecules around a point is described by a unit vector n : V ^ S^, called the director, and 
n and — n are assumed to be equivalent. 

This description is convenient for nematics and cholesterics. We will consider the director 
as a map n : V —>■ M'^, and we will show that the condition ||n|| = 1 is preserved by the 
Euler-Poincare dynamics. 

We shall obtain the Ericksen-Leslie equation by Euler-Poincare and Lie-Poisson reduction. 
For nematic and cholesteric liquid crystals, the order parameter Lie group is O = 5*0(3). In 
this paragraph we always identify the Lie algebra so (3) with M.^, that is, we have ad^ v = uxv. 
Identifying the dual through the canonical inner product on R^, we have ad* if = w x u. 

The symmetry group is the semidirect product DiS{T>) (§)!F{V, S0{3)). In the original 
Ericksen-Leslie approach, the liquid crystal flow is supposed to be incompressible. In this case 
the subgroup DiS(V)^o\®^{1^, S0{3)) should be used. Here we treat general compressible 
flows. An element (?], x) in the group DiS(V) (§)J^(V, S0{3)) acts linearly on the advected 
quantities (p, n) G J^iV) x J^{V,R^), by 

(p, n) t-^ {Jri{p o T]), x"^(n o r])) . 

The associated infinitesimal action and diamond operations are given by 

nu = Vnu, nv = nxiy, m Oi n = — Vn^m and m02n = nxm. 
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where t*, m, n G T{T>, M^). We use here and in the rest of the paper Vn instead of the usual 
derivative dn of the director n G JF(X),M^), since this notation is standard in the hquid 
crystals literature. Thus Vn is a 3 x 3 matrix whose rows are the vectors Vni, Vn2, Vna and 
hence its columns are 9in, 92n, Ssn; Vn^ denotes the transpose of Vn. The Euler-Poincare 
equations (O) associated to the group (Diff(P) (8)^^(1), 50(3))) (S) {J^iV) x J^iV.W)) are 



d 51 ^51 51 51 ^ 51 r 

— — = -i^u- divu- — du + pd- Vn • — 

at ou ou ou 5u dp \ on 



d 51 51 (51 \ 51 

— = 1/ X div — u + n X — . 

dt5v> 5v> \5h> I on' 



(8.51) 



Note that the first equation is in Vt^{T>), the dual of the Lie algebra X(P). The advection 
equations are 

d 

—p + div(pu) = 0, 

d (8-52) 
— n + (Vn)u + n X = 0. 

The reduced Lagrangian for nematic and cholesteric liquid crystals is of the form 

l{u,u,p,n) :=]- I p||uf/i + i / pJ\\u\\^p- [ pF(p"\ n, Vn)p, (8.53) 



2 



V ^ Jv Jv 



where the constant J is the microinertia constant and F is the free energy. The axiom of 
objectivity requires that 

F{p-\ A-^n, A-^WnA) = F{p-\ n, Vn), 

for all A G 0(3) for nematics, or for all A G 5*0(3) for cholesterics. 

A standard choice for F is the Oseen-Zocher- Frank free energy given by 

pF{p~^, n, Vn) = K2{n ■ curln) (divn)^ +^^22 (n ■ curln)^ 




twist 



^.54) 



bend 



where K2 7^ for cholesterics and K2 = for nematics. The free energy can also contain 
additional terms due to external electromagnetic fields. The constants i^'n, ^^22, -^^33 are 
respectively associated to the three principal distinct director axis deformations in nematic 
liquid crystals, namely, splay, twist, and bend. In general, these constants are different, but 
the expression of the resulting equations of motion can be greatly simplified if we make the 
one-constant approximation Ku = K22 = -^"33 = K. In this case the free energy becomes, 
up to the addition of a divergence, 

X 

pF(p-\ n,Vn) = -K\\Vnf. 
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The functional derivatives of the Lagrangian fl8.53p are computed to be 

51 , 51 
m := — = pu , K := — = pJu, 
ou 5u 



51 1, „2 1 „2 ^ 1 , 51 dF ^ f dF 



and 



5p 2 2 p dp^^ 5n (9n \ dn 

The vector field h is referred to as the molecular field. In the case of the free energy (18.541) 
for nematics {K2 = 0), the vector h is given by 

h = Kii graddivn — _ft'22(^curln + curl(^n)) + -^'33(6 x curln + curl(n x B)), 

where A := n-curl n and B := n x curln. In the case of the one-constant approximation we 
have h = —KAn. 

Using the Lagrangian (I8.53p . the Euler-Poincare equations (I8.5ip become 



( fd „ \ , dF ^ f dF ^ 

< 



pJ^i^ = h X n, 

at 



1.55) 



where we have used the standard notation ^ := ^ + u-V for the material derivative acting 
on every component of u. The advection equations are 

d 

—p + div(pu) = 0, 

D (8-56) 
— n = v X n. 

dt 

The evolution of the advected quantities is given by 

p = J{ri'^){pQ o 77-^) and n = (^Hq) o rj'^ . 

We now show that under some conditions, equations (I8.55P and (18.561) are equivalent to 
the Ericksen-Leslie equations for liquid crystals. This will use the following lemma. 

Lemma 8.3 Let v and n he solutions of the Euler-Poincare equations (I8.55P and (18.561) . 

Then 

(i) II no II = 1 implies ||n|| = 1 for all time. 

(ii) — (n-i/) = 0. Therefore, no-i^o = implies n u = for all time. 
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(iii) Suppose that no-i^o = and ||no|| = 1. Then the second equation of fl8.56p reads 



D 

u = n X — n 

dt 



and the second equation of f l8.55p reads 



pj— n - 2gn + h = 0, 



|2 



w/iere 2g := n ■ h — pJ II = n ■ h — pJ ||i/|| . 

Proof, (i) This is clear from the evolution n = (xng) o r^^^, since x ^ SO (3). 
(ii) Using the second equations in (I8.55P and (I8.56p . we have 

D , ^ f D \ [D \ , , 1 

n w) = I — n •!/ + n- 1 —ly = [v x nj-v + n — -(h x n) = 0. 



dt \dt J \dt J pJ 

(iii) Using the first two results, we obtain 

n X — n = n X (i/ X n) = (n n)!^ — iri u)ri = u 

dt 

which proves the first relation. 

To prove the second, we take the material time-derivative of the relation above to rewrite 

the second equation of fl8.55p as pJ {ri x = h x n. Taking the cross product with n 

on the left we obtain the equation 

Defining 2q := n-h + pJ (^n-^nj , this equation reads 

pJ— — n — 2qn + h = 

dt'^ 

which is the second relation since n-£^n = — ||-^||^ = — II^^IP- The last two equalities 
are obtained by taking two consecutive material time-derivatives of ||n|p = 1 and using the 
identity i/ = n x ^n, respectively. ■ 

Thus, we obtain the following theorem. 

Theorem 8.4 Let (u, i/,p, n) be a solution of the equations f l8.55p and fl8.56p . with the 
initial conditions no and v>q verifying 

||no|| = 1 and no-^'o = 0. 
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Then (u, p, n) is a solution of the Ericks en- Leslie equations 



grad 



dF 



dj p 



dF 



■Vn 



2gn + h = 0, 



(8.57) 



d_ 
dt 



p + div(pu) = 0, 



where 2g := n ■ h — pj H'^lT- 

Conversely, let (u, p, n) be a solution of the Ericks en- Leslie equations fl8.57p such that 
||n|| = 1, and define 



V :- 



D 

n X — n. 

dt 



Then (u, i/jp, n) is a solution of the equations fl8.55p and f l8.56p . 

Remark. One should think of the function q in the Ericksen-Leshe equation the way one 
regards the pressure in ideal incompressible homogeneous fluid dynamics, namely, the q is 
an unknown function determined by the imposed constraint ||n|| = 1 in the following way. 
The dot product of the second equation of (18.571) with n yields the formula of q given in the 
statement of the theorem by imposing ||n|| = 1. 

This q does not appear in the Euler-Poincare formulation relative to the variables (u, p, n) 
since in this case, the constraint ||n|| = 1 is automatically satisfied by Lemma [8.3( i). 

As a consequence of Theorem 18. 4[ we obtain the Ericksen- Leslie equations for liquid 
crystals by Lagrangian reduction. Consider the right-invariant Lagrangian 

: T [Diff (P) ®T{V, SO{m ^ K 

induced by the Lagrangian (18.531) . and assume that ||no|| = 1 and i/q'^o = 0. A curve 
(r/, x) G Diff ("D) (s) JF(P, 50(3)) is a solution of the Euler-Lagrange equations associated to 
-Z^(po,no)5 with initial condition if and only if the curve 

(u, v) := (?) o ?7~\ XX"^ o ^7"^) 

is a solution of the Ericksen-Leslie equations (I8.57p . where 

p = J(r^-i)(po or^"^) and n = (xno) o 7^"^ 

As in the case of microfluids, the curve t] G Diff(D) describes the Lagrangian motion of 
the fluid or macromotion, and the curve x ^ J-'{T>, S0{3)) describes the local molecular 
orientation relative to a fixed reference frame or micromotion. A standard choice for the 
initial value of the director is 



no(x) := (0,0,1), for all x eV. 



In this case we obtain that 



n 



Xl3 
X23 
X33 
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This relation is usually taken as a definition of the director, when the 3-axis is chosen as the 
reference axis of symmetry. 



By the Legendre transformation, the Hamiltonian for liquid crystals is given by 
/i(m,K,p,n) := i / i||m|p;U + [ -\\K,\\'^fi+ [ pF(p"^, n, Vn)/i. 



The Poisson bracket for liquid crystals is given by 



V 



{/,fi'}(m,p, K,n) = / m 
Jv 



(5m' 6m 



+ 



V 



6f 6f\6g 
OK, dm / on 



The Kelvin circulation theorem for liquid crystals reads 



d 
dt 



-Vn^h, 

ct P 



where h is the molecular field defined by 



, dF ( dF 
on \ on , 



+ K ■ [ — X - — h a— • d— ■ - — p 

Jd \ok ok ok om ok om 



^Sg 6g\ Sf 

OK dm / on 



p. 



8.11.2 Micropolar theory of liquid crystals 



This approach is based on the equations for micropolar continua given in (18.411) . The dif- 
ference from the micropolar fluid treated before is that for liquid crystals the free energy \1/ 
depends also on a new variable 7 = (7^ ) G n\V,so{3)) called the wryness tensor. This 
variable is denoted by 7 = (7") when it is seen as a form with values in M^. 
The constitutive equations in the non-dissipative case are given by 



tkl 



dp~ 



T^ki - pTT-^lt and rriki = p 



according to (12.5.18) in [Uj, where the function \l/ = \l>(p ^, j, 7) : M x Sy'm{3) x ^[(3) R 
denotes the free energy. The axiom of objectivity requires that 

*(p-\A-^jAA-V) = ^(p"\j,7), 

for all A G 0(3) (for nematics and nonchiral smectics), or for all A G 5*0(3) (for cholesterics 
and chiral smectics). See paragraphs 12.6, 12.8 and 12.9 in [ll| for the choice of the free 
energy for nematics, cholesterics and smectics respectively. 
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Assuming that fi = and Ik = 0, the equation fl8.4ip for micropolar continua become 

— p + p dlV U = 0, — jfci + {Skprjlp + SlprjkpjT^r = 0. 

These are the equation for non-dissipative micropolar liquid crystals as given in section 12 



of 11 



To these equations one needs to add the evolution of 7 given by 

fpt = dii^a + yaa\ - 7"5/U. (8.60) 
which is equation (12.3.13) in pjj]. Here ^ acts on the one-form 7 G fi^(T',]R^) as £7 : = 

|7+ At 

Lagrangian and Hamiltonian formulation of micropolar liquid crystals. We now 

show that the system of equations (18.591) and (18.601) can be obtained by affine Euler-Poincare 
and affine Lie-Poisson reduction. As in the case of micropolar fluids, the symmetry group is 
the semidirect product 0111(1?) (^T{T), 50(3)). The advected quantities are the mass density 
p, the microinertia tensor j G J-'iV, Sym{3)), and the wryness tensor 7 G fl^{V,so{3)). The 
action of the symmetry group on the variables p and j is linear and is the same as for 
micropolar fluids. The action of (77, x) on the wryness tensor 7 is affine and is given by 

The Lagrangian is the same as that for micropolar fluids, except for the fact that the free 
energy depends also on the wryness tensor. We thus have 

l(u,u,p,j,-f) = ^ / p||u||V + ^ / p{ju-u)fi- / p^(p-\j,7)p. (8.61) 
^ Jv ^ Jv Jv 

The computation of the affine Euler-Poincare equations is similar to that for micropolar 
fluids, except for the equation associated to the variable u for which we give some details 
below. 

Using the Lagrangian fl8.6ip we find that the evolution for v is 



where the index A denotes the antisymmetric part of the matrix. According to our con- 
ventions, ^ji^ on the right hand side means that ^ is applied to every entry of the matrix 
valued function : P — > so(3). Using the equality 

= 2 

97, (97i ' 
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we obtain ^ 

where ~ denotes the inverse of ^. We now use the axiom of objectivity to simphfy this 
expression. 

Lemma 8.5 Suppose that the free energy \I/ verifies the axiom of objectivity, that is, 
'^{p-\A-HA,A-^-rA) = ^{p'\i,-f), for all A e 50(3) (or 0(3)). 

Then 



(8.62) 



Proof. Consider a curve A{t) G 50(3) such that A{0) = h and A{0) = ^ G so (3). 
Differentiating the equahty '^{A{t)'^jA{t),A{ty^jA{t)) = ^(j, 7) at t = 0, we obtain the 
condition 

D^(j, 7)(je - - ^7) = 0, for all ^ G so(3). 

Using the equalities 

Dvi>(j, 7)(je - 0-,7e - ^7) = Tr (^^(j^ - O)) + Tr f (^^^ (7^ - ^7) 



Tr 



9^ 



and the identity 2 (^j^j = ~ fjj; "we obtain the result. 
Note that equation fl8.62p . can be rewritten in as 



Using these results for the equation for u, the affine Euler-Poincare equations associated to 
the Lagrangian (18.611) read 



,-.-(j.)x. = --div^p-J+7 x^, 



— p + div(pu) = 0, 

1-7 + £„7 + d^i> = 0. 

at 



(8.63) 
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Recall that (Pu = du + [7, u]. Thus, we have recovered equations (18.591) for non-dissipative 
micropolar liquid crystals, together with the equation (18.601) . up to a change of variables 
7 I — > —7. 

Consider the right-invariant Lagrangian -Z^(po,jo,7o) induced on the tangent bundle 
r[Diff(r') (8)^^(1), 50(3))] by the Lagrangian(IMI]). A curve (r^, x) e DiS{V) ®J^{V, SO (3)) 
is a solution of the Euler-Lagrange equations associated to jg^,,) if and only if the curve 

(u, i>) := {fj o r]-\ xx'' o ^ ^ X{V)®J^{V, S0{3)) 

is a solution of the equations (I8.63P with initial conditions (po; JojTo)- The evolution of the 
mass density p, the microinertia j, and the wryness tensor 7 is given by 

p = J(r^-i)(po o r/"^), j = {xjoX'^) °V~^, and 1 = V* {xioX~^ + xTx~^) ■ 
If the initial value 70 is zero, then the evolution of 7 is given by 

1 = V* {xTx~^) ■ 

This relation is usually taken as a definition of 7 when using equation (18.591) . We consider 
7 as an independent variable and therefore the system (18.631) contains an evolution equation 
for 7. 

The Legendre transformation and the Hamiltonian formulation can be carried out as 
in the case of micropolar fluids. The affine Lie-Poisson bracket consists of the sum of the 
Lie-Poisson bracket (I8.50p with the term 



V 



OK Sm J d'J \ OK J ()7 



due to the presence of the variable 7. 

The Kelvin-Noether circulation theorem applied to micropolar liquid crystals yields the 
relation 

The 7-circulation yields the relation 

d_ 
dt 

in R^. 



7 = (p U X 'J 

ct J ct 



Remark. According to equation (12.9.1) in HI], a liquid crystal flow is called smectic if the 
constraint 



i=l 



is satisfied. However, note that this constraint is not preserved by the evolution 

7 = V* ixioX'^ + xTx~^) 

in general. This is consistent with the fact that the last equation in (18.631) . which can be 
written equivalently in vectorial form as 

-f + £u7 + di^ + 7X i/ = 0, 

at 

does not imply that if the initial condition for 7 has trace zero then Tr^ = for all time. 
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8.11.3 Microstretch theory of polymeric liquid crystals 



This approach is based on the equations for microstretch continua given in (18.40p . The 
difference from the microstretch fluid treated before is that for polymeric liquid crystals the 
free energy \1/ depends also the wryness tensor 7 G Q^{T>,so{3)) and on the microstrain 
e e fl^{V). The constitutive equations in the non-dissipative case are given by 



(9^ 



rrik 



dp~ 
dm 

dck 



jSki - p 



(9^ 



(9^ 



rriki = p 



and 



t = 2p 



(9^ . . 

-77— Jkl + 77— Jo 
OJkl OJo 



see equations (16.3.13) and (16.3.13) in [ll(], where 

^ = ^(p-i, j,7,e) : M x Sym{3) x g[(3) x ^ M 

denotes the free energy. The axiom of objectivity requires that 

mip-\ A-'jA, A-'jA, A-'e) = mip~\j, 7, e), 

for all A G 0(3) (for nematics and nonchiral smectics), or for all A G 5*0(3) (for cholesterics 
and chiral smectics). See paragraphs 16.4, 16.6 and 16.7 in [llj] for the choice of the free 
energy for nematic, smectic, and cholesteric polymers respectively. 

Assuming that fi = and Ik = 0, the equation fl8.40p for microstretch continua become 



D ^ 



dm 



dm 



per = d; 



D 



dm\ 

dk P 



2p 



dYk 

dm . dm . 

TT—jkl + TT—jo 

djkl djo 



dt 



p + p div u = 



^ImnP 

D 



dm „ dm 



(8.64) 



dt 



jkl ~l~ i^^kprjlp ~l~ ^lprjkp)^r 



These are the equations for non-dissipative microstretch polymeric liquid crystals as studied 
in section 16 of 

To these equations one adds the evolution of 7 G r2^(P,]R'^) and e G Q^{V) given by 

dt 



and 



D 

-n&k — l^O,k + CjUj 

dt 



which is equation (16.3.8) in llj. Like in equation (18.601) . D/dt acts on the one-forms 7 and 
e as D/dt = d/dt + 

Lagrangian and Hamiltonian formulation of polymeric liquid crystals. We now 

show that the equation (I8.64p can be obtained by afline Euler-Poincare and afline Lie-Poisson 
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reduction. As in the case of microstretch fluids, the symmetry group is the semidirect product 
DiS(V) (§)J^(V,CS0{3)). The advected quantities are the mass density p, the microinertia 
tensor i G J^{V, Sym{3)), the wryness tensor 7 G Q^(V,so{3)), and the microstrain e G 
fl^{T>). The action of the symmetry group on the variables p and i is linear and is the same 
as for microstretch fluids. The action of {ri,x) G Diff (P) (s)^(P, CS'0(3)) on the wryness 
tensor 7 is affine and is given by 

7 ^ X'\v*l)x + X^^Tx, (8.65) 
where x G J^iV, 50(3)) is defined by the equahty 

X = det(x)'/'x, 

and is the microrotation part of x G J-^(T>, CS0{3)). The determinant det{xY^^ can be seen 
as the microstretch part of x- The action of {rj, x) on the microstrain e is also affine and is 
given by 

e ^r/*e + -det(x)"M(det(x)). (8.66) 
3 

We now explain why these affine actions are natural. The variables e and 7 can be seen as 
the symmetric and antisymmetric part of a strain tensor ( G c5o(3)). More precisely, 

we have 

C = e/3 + 7, 

where 3e := Tr(^) and 7 := (a- The affine action of Diff(T') (s)jF(P, 050(3)) is the nat- 
ural action of the automorphism group onto the connections of the trivial principal bundle 
050(3) X V ^ V, as defined in the general theory (see fl3.6l) ). that is, 

C ^ x''^vXx + X'^Tx- 

By taking the trace of this action we find the affine action (18.661) . whereas the antisymmetric 
part gives the affine action (18.651) . 

We now give the associated right infinitesimal actions and diamond operations. We have 

7U = £u7, 7z/=[7,i>], w Oi 7 = (divw)-7 - t(7-i_d7, w 02 = [ji, Wi] 
eu = i^uC, 61^ = 0, / Oi e = (div /) -e — /■i_de, / 02 e = 0. 

Relative to the two group one-cocycles 

Ci{x) = X-'Tx and 02(x) = ^ det(x)-M det(x) 

we have 

dOi(z/) = dr), dOf(w) = -dww 
d02(z/) = dz/o, d02^(/) = -i(div/)/3. 

The Lagrangian is the same as that for microstretch fiuids, except for the fact that the free 
energy depends also on the wryness tensor 7 and on the microstrain e. We thus have 

/(u,z/,p,i,7,e) = ^ / p||u||V+^ / p{iv-y)p- [ p^(p"\ 7, e)p, (8.67) 

where iu ■ 1/ = Tr((iz/)"^z/). The computation of the associated Euler-Poincare equations 
involves the following generalization of Lemma 18.51 
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Lemma 8.6 Suppose that the free energy \1/ verifies the axiom of objectivity, that is, 
'^{p-\A-HA,A-^-fA,A-^e) = ^{p-\i,-f,e), for all A e SO (3) (or 0(3)). 

Then 



2 J 



dj 



\dj J 



de 



(8.68) 



Using this lemma, a direct (but long) verification shows that the affine Euler-Poincare 
equations associated to the Lagrangian fl8.67p are 



Jo fD 



dt 



dt 



9^ 



77 IT'^o - '^o + (J*^)-^ = 2 joT^ + - -div p— 



div I p 



<9j / P 



9e / 



9^ 
^7 



+ 7 X 7^ Hex — . 

0*7" ae 



(8.69) 



at 
d_ 
di 



p + div(pu) = 0, 



D 



J + 2uoj + [j, 9] = 0, 



7 + £ul + [l, u]+du = 0, 



d_ 
di 



e + £-ae + dz/Q = 0. 



Thus, we have recovered equations (18.641) for non-dissipative polymeric liquid crystals to- 
gether with the conservation laws for 7 and e, up to a change of variables v ^ —v^ 
and 7 I—*- —7. Recall that the microstretch rate z/q G J-'{V) and the microrotation rate 
G J-'{T>,M.^) are constructed from the variable u G ^^(1), cso(3)) through the decomposi- 
tion 

As before, the Lagrangian reduction can be carried out, by starting with the right- 
invariant Lagrangian -^^(po,io,7o,eo) • T [DiS(V) (§)J^(V, CS0{3))] — > M induced by the La- 
grangian (18.671) . The evolution of the linear advected quantities p and i is the same as in 
the case of microstretch fluids. The evolution of the affine advected quantities is given by 

7 = V* (x7oX"^ + xTx~^) and e = ^Cq + ^ det(x)d (det(x)"^)^ • 
If the initial values are zero, then the evolution of 7 and e is given by 
1 = V* {xTx~^) and e = ^r]^(^ det(x)d (det(x)"^) 



ll| takes these relations as definitions of 7 and e. 

The Legendre transformation and the Hamiltonian formulation can be carried out as in 
the case of microstretch fiuids. The affine Lie-Poisson bracket consists of the sum of the 
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Lie-Poisson bracket (I8.44p with the two terms 



+ 



V 



SI 
6k 



OK 



(5m 



(5 m 



be 



Sk 



£sa.e I — 



(57 
5f 



due to the presence of the variables 7 and e. 



The Kelvin-Noether circulation theorem applied to polymeric liquid crystals yields the 
relation 



d 



(9* 



di 



(97 



(9^ 



-—■di + •i_d7 div p—- -7 + — — -ide div p^r— e. 



The 7-circulation yields the relations 



in M^, and 



d 

dt^^ 



d 
dt 



(97 



1/ X 7 



9e 



(9* 



(9e 



e = 0. 



Ct 



8.11.4 Ordered micropolar theory 

This approach is developed in |25|. It is based on the micropolar theory and use the Oseen- 
Zocher- Frank free energy. As we will see, it gives a direct generalization of the Ericksen- Leslie 
equations. 

Lhuillier and Rey consider the general equations for micropolar continua fl8.4ip . together 
with the constitutive relations in the non-dissipative case, 

dF ^ OF „ ^ OF 



see equations (10) in |25|. Note that here we adapted these relations to the compressible 
case. Thus, equations f l8.4ip read 



D 



dF 



OF 



dp 



OF 

pai = dk I SlabnapTTT- 



dF 
<9n" 



-n 



(8.70) 



— p + p div U = 0, —jkl + [Skprjlp + eiprjkpj^r = 0. 



To these equations one needs to add the evolution for n 

D 



dt 



-n = 1/ X n. 
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which is equation (9) in [25|]. Recall that a = ^{ju). 

Lagrangian and Hamiltonian formulation of ordered micropolar theory of liquid 
crystals. We now show that the equations fl8.70p can be obtained by Euler-Poincare and 
Lie-Poisson reduction. As in the case of micropolar fluids, the symmetry group is the semidi- 
rect product DiS(V) (§)J^{V, S0{3)). The advected quantities are the mass density p, the 
microinertia tensor j G J^{T>^ Sym{3)), and the director n G JF(D,R^). The representation 
of the symmetry group on the variables p and j is the same as for micropolar fluids. The 
representation on the director is the same as that for the director theory. 

The Lagrangian is the same as that for micropolar fluids, except the fact that it involves 
the elastic free energy F(p~^,n, Vn), which is usually taken to be the Oseen-Zocher- Frank 
free energy. We thus obtain 

^(u,i^,p, j,n) = i [ p\\ufp + ^ I p{jiy-iy)^i- [ pF(p"\ n, Vn)/i. 
Jv Jv Jv 



The associated Euler-Poincare equations are 

d 



P 



dt 



-u 



The advection equations are 



VuU 



h X n. 



dF 



dF 
'dn_i 



■Vn 



(8.71) 



d_ 
dt 
D 

di 
D 
dt 



p + div(pu) = 0, 
J + [j, u] = 0, 



(8.72) 



v X n. 



The evolution of the advected quantities is given by 

p = J(r^-i)(po 077"^), j = (xjoX^^) <=> V'^ and n = (xno) o 77"^ 

Thus we have recovered equations fl8.70p together with the evolution of the director. One 
just needs to prove that the equations for the variable u are equivalent in fl8.7ip and fl8.70p . 
As shown in the following lemma, this is a consequence of the axiom of objectivity. 

Lemma 8.7 Suppose that the free energy F verifies the property 

F{p-\A-^n,A-WnA) = F{p-\n,Vn), for all A e SO (3) (or 0(3)). 

Then the matrix 



dF 

avn 



(Vn) - (Vn) 



dF 
dVn 



n 



dF 
dn 



is symmetric. 
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Note that in this expression the last term denotes the multiphcation of a column vector 



therefore, the equations associated to the variable v in f l8.70p and (18.711) are equivalent. 

The Hamiltonian and Lie-Poisson bracket can be computed as in the preceding examples. 
The Kelvin-Noether circulation theorem has the same form as that of the Ericksen- Leslie 
equations, namely 



8.11.5 Comparison of the three theories 

In this subsection we summarize the known relationships between the three theories for liquid 
crystals presented in this paper. We shall prove that the director theory of Ericksen-Leslie 
is a particular case of the ordered micropolar theory of Lhuillier-Rey. Therefore, one needs 
to compare the latter with the micropolar theory of Eringen. As will be shown, these two 
theories, while close, do not seem to be equivalent. 

Theorem 8.8 The Ericksen-Leslie equations are a particular case of the equations given by 
the ordered micropolar theory. 

Proof. As we have seen in Theorem 18.41 if ||n|| = 1 and ly = n x ^n, then the Ericksen- 




Leslie equations (18.571) are equivalent to the equations (18.551) and (18.561) . It turns out that 
these equations are a particular case of the equations (18.711) and (I8.72p given by the ordered 
micropolar theory. To see this, it suffices to assume that the microinertia j is constant and 
given by j = JI3, where J is the microinertia constant appearing in the director theory, and 
J3 is the identity 3x3 matrix. ■ 

Thus, the ordered micropolar theory is a generalization of the director theory, which takes 
into account the variation of microinertia. We now compare the ordered director theory 
with the micropolar theory of Eringen for a particular choice of the initial condition for the 
microinertia j. This choice imposes the condition that j is a moment of inertia. This will 
use some technical lemmas. 

Lemma 8.9 Let j and n be solutions of the equations (I8.7ip and (I8.72p . Let jo and no be 

the initial values and suppose that 



with a row vector and the result of this operation is a 3 x 3 matrix whose {i, j) entry is n. 
Using this lemma, we obtain the equality 



dF 
■« dni 






Jo 



J(/3 - no ® no) 



where J is a scalar constant. Then 



J 



J(J3 - n ® n) 



for all time. 
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Proof. From the Lagrangian formulation, we know that the evolution of j and n is 

n = (xno) o and j = {xjoX'^) ° V~^- 
Since x(n ® n)^^-'^ = (xn) ® (xn) for all x ^ S0{3), we obtain 

j = {xjoX'^) o V'^ = J{h- (xno) ® (xno)) = J(/3 - n ® n). ■ 



Lemma 8.10 Let u, j, and n he solutions of the equations (18.711) and (18.721) . Define 

7 := (Vn) X n G 9}{V,W^), 
that is, 7(v) = [(Vn)v] x n G T{T>^M?) for every v G X(P). T/ien 7 verifies the equation 



d j 

77:7 + "Cu7 + 7 X + ^di/ = 0. 



Proof. Using the equation 
and ||n|| = 1, we compute 



D 
It 



-n = X n. 



D „ [D \ ^ D 

—7 = V — n X n + Vn x — n 

dt \dt J dt 

= {du X n) X n + (i^ X Vn) x n + Vn x (1/ x n) 
= (ndi>')n — di/ + (ni/)Vn — {Vni/)n 
= u X (Vn X n) — {dv — (n-di/)n) 
= ^'X7 — (13 — n® n)di>'. ■ 

Theorem 8.11 Suppose (see the following lemma) that the free energy F can be written in 
terms of 

j := J{I'i — n ® n) and 7 := Vn x n, 
that is, there exists a function \1' : M x Sym{3) x 0[(3) — M such that 

F{p-\n, Vn) = ^(p-\j,7). 



Let (u, I/, p, j, n) be a solution of the equations (18.711) and (18.721) . and suppose that the initial 
conditions verify Jq = J{Iz — no ® ng). Then (u, p,j, 7) is a solution of the system 

x. = --cUv(^,--j+T, X—, ^^^^^ 



d 

—p + div(pu) = 0, 

d j 

77:7 + -^u7 + 7 X + ^di/ = 0. 

at J 



This system is very close but distinct from the system (I8.63P studied by Eringen. The differ- 
ence is due to the presence of the factor j in the second and last equations. 
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Using that the relation 7 = Vn x n reads 7^ = el^n\n'^ in coordinates, we obtain 



Proof. Let {u,u, p, j,n) be a solution of the equations fl8.7ip and (18 ■721) . The equations 
associated to the conservation of mass and microinertia in (18.741) are clearly verified since 
they are identical. Using that jo = J{h — no ® no), we obtain j = J{I^ — n n), by lemma 
18.91 Thus by Lemma 18. 101 the last equation in fl8.74p is verified. 

dF _ d-^ dYi _ 
dn\ dn^, d'j' 

This shows that 

Thus, the first equation is verified. We also have 

= -(/a - n (g) n)ia- — = 



a 





9^ 










-( 


Jd^J 



Ik 



Using this equality and the relation (18.731) . the second equation of (I8.7ip reads 

a dF\ dF ^ 



This can be written as 



which is exactly the second equation in (18.740 . ■ 

No Lie-Poisson or Euler-Poincare interpretation of the system (18.740 is known. This 
system would coincide with Eringen's system (18.630 . if j/J could be taken equal to one. 
This, however, cannot be achieved since j is an advected variable; J is a constant. We 
interpret Theorem 18.111 as saying that the the Lhuillier-Rey and Eringen formulations are 
very close but not equivalent. 

We now show that the hypothesis made on the free energy F is verified in the case of the 
Oseen-Zocher-Frank free energy. 

Lemma 8.12 Let n a unit vector. Define 

7 := Vn X n and j = /s — n ® n. 
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Then 



Tr(7) = n-curln, 
Tr(7^7) = Tr(Vn^Vn) = (divn)^ + (n ■ curln)^ + ||n x curln||^ 

+ div((n ■ V)n — n div n), 

Tr(7'^7(n ® n)) = ||n x curl n||^. 



Thus, 



Tr(7'^7j) = (divn)^ + (n ■ curln)^ + div((n ■ V)n — ndivn). 
Proof. Note that the relation 7 = Vn x n reads 

li 



dill X n 



n dm 



Thus, we obtain 



n ■ curln = n^{d2n'^ — d^n^) + n^i^d^n^ — din^) + n^i^div? — d2n^) 

= {n^din'^ - n^din^) + {n^d2n^ - n^d2n^) + {n^dsn^ - 71^8371^) 
= 7}+72 + 73 = Tr(7). 

This shows the first equality. 

In order to check the second assertion, we compute 7"^7 in terms of n. We have 

7^(u) ■ V = u ■ 7(v) = u ■ (Vn(v) x n) = Vn(v) ■ (n x u) = v ■ Vn^(n x u), 

therefore, we have 7"^(u) = Vn-^(n x u) and we can compute 

7'^7(u) = Vn^[n x (7(u))] = Vn^[n x (Vn(u) x n)] 

= Vn^[(n ■ n)Vn(u)) - (n ■ Vn(u))n] = Vn'^Vn(u). 

We have obtained the equality 7"^7 = Vn-^Vn. The formula 

Tr(Vn'^Vn) = (divn)^ + (n ■ curln)^ + ||n x curln||^ + div((n ■ V)n — ndivn) 



can be checked directly, see for example Lemma 3.3 in j33|. This proves the second equality. 
The third assertion follows from the equalities 

Tr (7^7(n (g) n)) = Tr ((7n)^7n) = ||7np = ||(Vn)n x n||^ 



||(Vn)n| 



n X rot n 



where the last one follows from the identity (Vn)n = — nxrot n; see Lemma 3.3 in 33 



As a consequence of this lemma, the expressions associated to chirality, twist, splay, and 
bend appearing in the Oseen-Zocher-Frank free energy can be expressed in terms of the 
variables j and 7. We have 



n-curln = Tr(7), 



fn- curln) 



Tr(7)^ 



(divn)^ = Tr(7"^7j) — Tr(7)^, modulo a divergence, 
|n X curln||^ = Tr(7"^7) — Tr(7"^7j). 
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Thus, in terms of j and 7, the Oseen-Zocher-Frank free energy reads 

7) =i^2Tr(7) + ^Ku (Tr(7^7j) - Tr(7)') + ^K.^Tiiyf 

+ ^Kss (Tr(7^7)-Tr(7^7j)), 
up to the addition of a divergence. This functional clearly satisfies the axiom of objectivity. 

8.11.6 Remark on the use of other groups 

From a mathematical point of view, the previous approaches generalize to any order param- 
eter Lie group O C GL{3)~^. 

7-theory. In the case of Eringen's theory, it suffices to consider the group Diff (D) (s)jF(X), O) 
acting on the advected quantities 

(p, i, C) e J^{V) X Sym{3) x n\V, 0) 

as 

If O = 5*0(3) then = 7 is the wryness tensor, and we recover the theory of micropolar 
liquid crystals (see §8.11.21) . If O = CS0{3) then ( = (7,e), where 7 is the wryness tensor 
and e is the microstrain and we recover the microstretch theory of polymeric liquid crystals 
(see §8.11.31) . In general, we obtain a theory of ^^O-liquid crystals''^ whose Lagrangian given 
by 

/(u,z/,p,z,C) = - / p||ufp + 7j / p{ii'-i')f^- / p^'(p"\i,C)/^, 
Jv Jv Jv 

and the variable C, is interpreted as a connection on the trivial (9-principal bundle over T>. 

The associated affine Euler-Poincare equations (13.131) are 



' p(|u + V.u 



5^ / (9^ 



d d'^ 
P(,(g^. + d.(u)-..-2- 



0, 



d_ 

di 
d_ 

di 



p + div(pu) = 0, 



d_ 

di 



1.75) 



i + dz(u) + z/ i + iu = 0, 



C + i:uC+[C,'^]+di/ = 0, 



where P : 0t(3) — denotes the orthogonal projection onto the Lie algebra 0, associated to 
the inner product 

a-b = Tr(a'^6), a,b E 0. 

n-theory. Recall that in the case of the Ericksen-Leslie and Lhuillier-Rey theories the 
director is a map n : V ^ M^, on which Diff (P) (S)JF(T', 5*0(3)) acts linearly as 



in^ X (nor/). 



(8.76) 



101 



This representation can be generalized to other groups in two ways. First the action (18.761) 
clearly still makes sense for any matrix Lie group O. In this case, the advection equation for 
n reads 

D 

— n = un, 

dt 

where u G J^{T>, o), and the variable n evolves by 

n = (xno) oT]-'^. 

Note that, in general, the norm ||n(a;)|| is not constant. For example, if O = SO{K), for 
K positive definite, then n describes an ellipsoid. As in the case of microfluids, the typical 
cases to consider are 



1 
1 

e 



or Ko 



€ 
e 
1 



< e < 1. 



The group SO{Ki) should be useful for the description of smectic liquid crystals. If O = 
CS0{3), then the norm ||n|| evolves as 

||n(r/(x))|| =det(x(x))^/='||no(x)||. 

Therefore, the norm depends on the "stretch part" det(x)^/'^ of x ^ CS0{3). More generally, 
we can use the group CSO{K) to model the theory of anisotropic microstretch liquid crystals. 
This constitutes a first way of generalizing the representation f l8.76p to other groups. 

A second way to generalize the representation (18.761) consists in replacing the director n : 
D — > by a Lie algebra valued variable n : — > o on which the group Diff (P) (S)JF(D, O) 
acts as 



n 



Ad^-i(r;, o r/). 



(8.77) 



In this case, the advection equation reads 



D 



—n = adiyU, 
dt 

and the variable n evolves as 

n = (Ad^ no) o r]~^. 

Remark that when O = 5*0(3), both approaches coincide because the birth representa- 
tion and the adjoint representation of S0{3) on are identical. When O = CS0{3), 
(microstretch case) then we can write n = n + m/3 and the action (18.771) reads 



where 



(n,m)i-^(x {n o rj) , m o rj) , 
1 



X 



det(x)i/3 



X--V-^S0{3), 



and the evolutions are given by 



n 



(xno)o?^ ^ and m = niQ o r]~ 
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Therefore n can be seen as a director, since ||n(a;)|| is constant in time. Since the evolution 
of the variable m does not depend on the micromotion and the evolution of n depend only 
on the S0{3) part of the micromotion x £ CSO{3), this approach cannot be used for the 
description of microstretch liquid crystals. This shows that the first generalization seems 
physically more interesting. 
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